Nektar::FilterCellHistoryPoints Class Reference

#include <FilterCellHistoryPoints.h>

Inheritance diagram for Nektar::FilterCellHistoryPoints:

Collaboration diagram for Nektar::FilterCellHistoryPoints:

Public Member Functions
	FilterCellHistoryPoints (const LibUtilities::SessionReaderSharedPtr &pSession, const ParamMap &pParams)
	~FilterCellHistoryPoints ()
void	SetCellModel (CellModelSharedPtr &pCellModel)
Public Member Functions inherited from Nektar::SolverUtils::FilterHistoryPoints
SOLVER_UTILS_EXPORT	FilterHistoryPoints (const LibUtilities::SessionReaderSharedPtr &pSession, const ParamMap &pParams)
SOLVER_UTILS_EXPORT	~FilterHistoryPoints ()
Public Member Functions inherited from Nektar::SolverUtils::Filter
SOLVER_UTILS_EXPORT	Filter (const LibUtilities::SessionReaderSharedPtr &pSession)
virtual SOLVER_UTILS_EXPORT	~Filter ()
SOLVER_UTILS_EXPORT void	Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
SOLVER_UTILS_EXPORT void	Update (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
SOLVER_UTILS_EXPORT void	Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
SOLVER_UTILS_EXPORT bool	IsTimeDependent ()

Static Public Member Functions
static SolverUtils::FilterSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const ParamMap &pParams)
	Creates an instance of this class. More...
Static Public Member Functions inherited from Nektar::SolverUtils::FilterHistoryPoints
static FilterSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const std::map< std::string, std::string > &pParams)
	Creates an instance of this class. More...

Static Public Attributes
static std::string	className
	Name of the class. More...
Static Public Attributes inherited from Nektar::SolverUtils::FilterHistoryPoints
static std::string	className = GetFilterFactory().RegisterCreatorFunction("HistoryPoints", FilterHistoryPoints::create)
	Name of the class. More...

Protected Member Functions
virtual void	v_Update (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
Protected Member Functions inherited from Nektar::SolverUtils::FilterHistoryPoints
virtual SOLVER_UTILS_EXPORT void	v_Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
virtual SOLVER_UTILS_EXPORT void	v_Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)
virtual SOLVER_UTILS_EXPORT bool	v_IsTimeDependent ()

Protected Attributes
CellModelSharedPtr	m_cell
Protected Attributes inherited from Nektar::SolverUtils::FilterHistoryPoints
SpatialDomains::PointGeomVector	m_historyPoints
unsigned int	m_index
unsigned int	m_outputFrequency
unsigned int	m_outputPlane
	plane to take history point from if using a homogeneous1D expansion More...
bool	m_isHomogeneous1D
std::string	m_outputFile
std::ofstream	m_outputStream
std::stringstream	m_historyPointStream
std::list< std::pair < SpatialDomains::PointGeomSharedPtr, Array< OneD, NekDouble > > >	m_historyList
std::map< int, int >	m_historyLocalPointMap
Protected Attributes inherited from Nektar::SolverUtils::Filter
LibUtilities::SessionReaderSharedPtr	m_session

Friends
class	MemoryManager< FilterCellHistoryPoints >

Additional Inherited Members
Public Types inherited from Nektar::SolverUtils::Filter
typedef std::map< std::string, std::string >	ParamMap

Detailed Description

Definition at line 45 of file FilterCellHistoryPoints.h.

Constructor & Destructor Documentation

Nektar::FilterCellHistoryPoints::FilterCellHistoryPoints	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const ParamMap &	pParams
	)

Definition at line 52 of file FilterCellHistoryPoints.cpp.

                             :
    FilterHistoryPoints(pSession, pParams)
{
}

Nektar::FilterCellHistoryPoints::~FilterCellHistoryPoints

(

)

Definition at line 63 of file FilterCellHistoryPoints.cpp.

{

}

Member Function Documentation

static SolverUtils::FilterSharedPtr Nektar::FilterCellHistoryPoints::create ( const LibUtilities::SessionReaderSharedPtr &  pSession, const ParamMap &  pParams  )

inlinestatic

Creates an instance of this class.

Definition at line 51 of file FilterCellHistoryPoints.h.

References Nektar::MemoryManager< DataType >::AllocateSharedPtr().

                                     {
            SolverUtils::FilterSharedPtr p
                            = MemoryManager<FilterCellHistoryPoints>
                                    ::AllocateSharedPtr(pSession, pParams);
            return p;
        }

void Nektar::FilterCellHistoryPoints::SetCellModel ( CellModelSharedPtr &  pCellModel )

inline

Definition at line 68 of file FilterCellHistoryPoints.h.

References m_cell.

Referenced by Nektar::Monodomain::v_InitObject().

        {
            m_cell = pCellModel;
        }

void Nektar::FilterCellHistoryPoints::v_Update ( const Array< OneD, const MultiRegions::ExpListSharedPtr > &  pFields, const NekDouble &  time  )

protectedvirtual

Reimplemented from Nektar::SolverUtils::FilterHistoryPoints.

Definition at line 72 of file FilterCellHistoryPoints.cpp.

References Nektar::iterator, m_cell, Nektar::SolverUtils::FilterHistoryPoints::m_historyList, Nektar::SolverUtils::FilterHistoryPoints::m_historyLocalPointMap, Nektar::SolverUtils::FilterHistoryPoints::m_historyPoints, Nektar::SolverUtils::FilterHistoryPoints::m_index, Nektar::SolverUtils::FilterHistoryPoints::m_isHomogeneous1D, Nektar::SolverUtils::FilterHistoryPoints::m_outputFrequency, Nektar::SolverUtils::FilterHistoryPoints::m_outputPlane, Nektar::SolverUtils::FilterHistoryPoints::m_outputStream, and Nektar::LibUtilities::ReduceSum.

{
    // Only output every m_outputFrequency.
    if ((m_index++) % m_outputFrequency)
    {
        return;
    }

    int j         = 0;
    int k         = 0;
    int numPoints = m_historyPoints.size();
    int numFields = m_cell->GetNumCellVariables();
    LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
    Array<OneD, NekDouble> data(numPoints*numFields, 0.0);
    Array<OneD, NekDouble> gloCoord(3, 0.0);
    std::list<std::pair<SpatialDomains::PointGeomSharedPtr,
                        Array<OneD, NekDouble> > >::iterator x;

    Array<OneD, NekDouble> physvals;
    Array<OneD, NekDouble> locCoord;
    int expId;
    int nppp = 0; // Number of points per plane

    // Pull out data values field by field
    for (j = 0; j < numFields; ++j)
    {
        if(m_isHomogeneous1D)
        {
            for (k = 0, x = m_historyList.begin(); x != m_historyList.end();
                 ++x, ++k)
            {
                locCoord = (*x).second;
                expId    = (*x).first->GetVid();
                nppp     = pFields[0]->GetPlane(0)->GetTotPoints();

                physvals = m_cell->GetCellSolution(j) + m_outputPlane*nppp
                            + pFields[j]->GetPhys_Offset(expId);

                // interpolate point can do with zero plane methods
                data[m_historyLocalPointMap[k]*numFields+j]
                    = pFields[0]->GetExp(expId)->StdPhysEvaluate(
                                                    locCoord,physvals);

            }
        }
        else
        {
            for (k = 0, x = m_historyList.begin();
                 x != m_historyList.end();
                 ++x, ++k)
            {
                locCoord = (*x).second;
                expId    = (*x).first->GetVid();

                physvals = m_cell->GetCellSolution(j)
                            + pFields[0]->GetPhys_Offset(expId);

                // interpolate point
                data[m_historyLocalPointMap[k]*numFields+j]
                    = pFields[0]->GetExp(expId)->StdPhysEvaluate(
                                                    locCoord,physvals);
            }
        }
    }

    // Exchange history data
    // This could be improved to reduce communication but works for now
    vComm->AllReduce(data, LibUtilities::ReduceSum);

    // Only the root process writes out history data
    if (vComm->GetRank() == 0)
    {

        // Write data values point by point
        for (k = 0; k < m_historyPoints.size(); ++k)
        {
            m_outputStream.width(8);
            m_outputStream << setprecision(6) << time;
            for (int j = 0; j < numFields; ++j)
            {
                m_outputStream.width(25);
                m_outputStream << setprecision(16) << data[k*numFields+j];
            }
            m_outputStream << endl;
        }
    }
}

Friends And Related Function Documentation

friend class MemoryManager< FilterCellHistoryPoints >

friend

Definition at line 48 of file FilterCellHistoryPoints.h.

Member Data Documentation

std::string Nektar::FilterCellHistoryPoints::className

static

Initial value:

= SolverUtils::GetFilterFactory().RegisterCreatorFunction(
            "CellHistoryPoints", FilterCellHistoryPoints::create)

Name of the class.

Definition at line 61 of file FilterCellHistoryPoints.h.

CellModelSharedPtr Nektar::FilterCellHistoryPoints::m_cell

protected

Definition at line 78 of file FilterCellHistoryPoints.h.

Referenced by SetCellModel(), and v_Update().