40 #include <boost/format.hpp>
53 FilterHistoryPoints::FilterHistoryPoints(
58 ParamMap::const_iterator it;
61 it = pParams.find(
"OutputFile");
62 if (it == pParams.end())
68 ASSERTL0(it->second.length() > 0,
"Missing parameter 'OutputFile'.");
78 it = pParams.find(
"OutputFrequency");
79 if (it == pParams.end())
93 it = pParams.find(
"OutputPlane");
94 if (it == pParams.end())
106 it = pParams.find(
"Points");
107 ASSERTL0(it != pParams.end(),
"Missing parameter 'Points'.");
130 "No history points in stream.");
137 int dim = pFields[0]->GetGraph()->GetSpaceDimension();
146 NekDouble Z = (pFields[0]->GetHomogeneousBasis()
150 cout <<
"Reseting History point from " << gloCoord[2]
151 <<
" to " << Z << endl;
161 gloCoord[1], gloCoord[2]);
172 int vRank = vComm->GetRank();
178 std::vector<Array<OneD, NekDouble> > LocCoords;
183 for (i = 0; i < vHP; ++i)
191 idList[i] = pFields[0]->GetExpIndex(gloCoord,locCoords,
195 LocCoords.push_back(locCoords);
203 pFields[0]->GetExp(idList[i])->GetGeom();
207 for (
int j = 0; j < g->GetCoordim(); ++j)
209 e->BwdTrans(g->GetCoeffs(j), coordVals);
210 NekDouble x = e->PhysEvaluate(locCoords, coordVals)
225 for (i = 0; i < vHP; ++i)
227 if (dist_loc[i] == dist[i])
246 = pFields[0]->GetZIDs();
247 for(j = 0; j < IDs.num_elements(); ++j)
255 if(j != IDs.num_elements())
275 for (i = 0; i < vHP; ++i)
279 if (procList[i] != vRank)
299 if (vComm->GetRank() == 0)
301 for (i = 0; i < vHP; ++i)
310 + boost::lexical_cast<std::string>(gloCoord[0])
312 + boost::lexical_cast<std::string>(gloCoord[1])
314 + boost::lexical_cast<std::string>(gloCoord[2])
315 +
" cannot be found in the mesh.");
321 cout <<
"Warning: History point " << i <<
" at ("
322 << gloCoord[0] <<
"," << gloCoord[1] <<
","
323 << gloCoord[2] <<
") lies a distance of "
324 << sqrt(dist[i]) <<
" from the manifold." << endl;
332 for (i = 0; i < pFields.num_elements(); ++i)
346 for (i = 0; i < vHP; ++i)
382 int numFields = pFields.num_elements();
386 std::list<std::pair<SpatialDomains::PointGeomSharedPtr, Array<OneD, NekDouble> > >
::iterator x;
393 for (j = 0; j < numFields; ++j)
400 locCoord = (*x).second;
401 expId = (*x).first->GetVid();
403 physvals = pFields[j]->GetPlane(
m_outputPlane)->UpdatePhys() + pFields[j]->GetPhys_Offset(expId);
406 if(pFields[j]->GetPhysState() ==
false)
408 pFields[j]->GetPlane(
m_outputPlane)->GetExp(expId)->BwdTrans(pFields[j]->GetPlane(
m_outputPlane)->GetCoeffs() + pFields[j]->GetCoeff_Offset(expId),physvals);
412 data[
m_historyLocalPointMap[k]*numFields+j] = pFields[j]->GetExp(expId)->StdPhysEvaluate(locCoord,physvals);
420 locCoord = (*x).second;
421 expId = (*x).first->GetVid();
423 physvals = pFields[j]->UpdatePhys() + pFields[j]->GetPhys_Offset(expId);
426 if(pFields[j]->GetPhysState() ==
false)
428 pFields[j]->GetExp(expId)->BwdTrans(pFields[j]->GetCoeffs() + pFields[j]->GetCoeff_Offset(expId),physvals);
432 data[
m_historyLocalPointMap[k]*numFields+j] = pFields[j]->GetExp(expId)->StdPhysEvaluate(locCoord,physvals);
442 if (vComm->GetRank() == 0)
449 for (
int j = 0; j < numFields; ++j)
451 m_outputStream <<
" " << boost::format(
"%25e") % data[k*numFields+j];
464 if (pFields[0]->GetComm()->GetRank() == 0)
std::ofstream m_outputStream
#define ASSERTL0(condition, msg)
std::map< int, int > m_historyLocalPointMap
static boost::shared_ptr< DataType > AllocateSharedPtr()
Allocate a shared pointer from the memory pool.
static const NekDouble kVertexTheSameDouble
boost::shared_ptr< SessionReader > SessionReaderSharedPtr
std::stringstream m_historyPointStream
boost::shared_ptr< Comm > CommSharedPtr
Pointer to a Communicator object.
unsigned int m_outputPlane
plane to take history point from if using a homogeneous1D expansion
NekDouble Evaluate() const
SpatialDomains::PointGeomVector m_historyPoints
unsigned int m_outputFrequency
std::map< std::string, std::string > ParamMap
LibUtilities::SessionReaderSharedPtr m_session
SOLVER_UTILS_EXPORT ~FilterHistoryPoints()
virtual SOLVER_UTILS_EXPORT void v_Finalise(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
virtual SOLVER_UTILS_EXPORT void v_Initialise(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
StandardMatrixTag boost::call_traits< LhsDataType >::const_reference rhs typedef NekMatrix< LhsDataType, StandardMatrixTag >::iterator iterator
virtual SOLVER_UTILS_EXPORT bool v_IsTimeDependent()
static const NekDouble kGeomFactorsTol
FilterFactory & GetFilterFactory()
boost::shared_ptr< StdExpansion > StdExpansionSharedPtr
virtual SOLVER_UTILS_EXPORT void v_Update(const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
void Vcopy(int n, const T *x, const int incx, T *y, const int incy)
boost::shared_ptr< Geometry > GeometrySharedPtr
std::list< std::pair< SpatialDomains::PointGeomSharedPtr, Array< OneD, NekDouble > > > m_historyList
boost::shared_ptr< PointGeom > PointGeomSharedPtr
tKey RegisterCreatorFunction(tKey idKey, CreatorFunction classCreator, tDescription pDesc="")
Register a class with the factory.