#include <FilterHistoryPoints.h>

Inheritance diagram for Nektar::SolverUtils::FilterHistoryPoints:

Public Member Functions
SOLVER_UTILS_EXPORT	FilterHistoryPoints (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, const ParamMap &pParams)

SOLVER_UTILS_EXPORT	~FilterHistoryPoints () override

Public Member Functions inherited from Nektar::SolverUtils::Filter
SOLVER_UTILS_EXPORT	Filter (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation)

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 FilterSharedPtr	create (const LibUtilities::SessionReaderSharedPtr &pSession, const std::weak_ptr< EquationSystem > &pEquation, 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...

Protected Member Functions
SOLVER_UTILS_EXPORT void	v_Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time) override

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

SOLVER_UTILS_EXPORT void	v_Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time) override

SOLVER_UTILS_EXPORT bool	v_IsTimeDependent () override

bool	GetPoint (Array< OneD, NekDouble > gloCoord, int I)

SOLVER_UTILS_EXPORT void	v_WriteData (const int &rank, const Array< OneD, NekDouble > &data, const int &numFields, const NekDouble &time)

virtual void	v_Initialise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)=0

virtual void	v_Update (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)=0

virtual void	v_Finalise (const Array< OneD, const MultiRegions::ExpListSharedPtr > &pFields, const NekDouble &time)=0

virtual bool	v_IsTimeDependent ()=0

Protected Attributes
Array< OneD, Array< OneD, const NekDouble > >	m_historyPoints

size_t	m_historyPointsSize = 0

unsigned int	m_index = 0

unsigned int	m_outputFrequency

int	m_outputPlane
	plane to take history point from if using a homogeneous1D expansion More...

std::vector< int >	m_planeIDs

bool	m_isHomogeneous1D

bool	m_waveSpace

std::string	m_outputFile

std::ofstream	m_outputStream

std::stringstream	m_historyPointStream

std::list< std::tuple< Array< OneD, const NekDouble >, Array< OneD, const NekDouble >, int, int > >	m_historyList

std::map< LibUtilities::PtsType, Array< OneD, NekDouble > >	m_pointDatMap

std::map< LibUtilities::PtsType, Array< OneD, int > >	m_pointNumMap

unsigned int	m_outputIndex = 0

bool	m_outputOneFile

bool	m_adaptive

Protected Attributes inherited from Nektar::SolverUtils::Filter
LibUtilities::SessionReaderSharedPtr	m_session

const std::weak_ptr< EquationSystem >	m_equ

Friends
class	MemoryManager< FilterHistoryPoints >

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 FilterHistoryPoints.h.

Constructor & Destructor Documentation

◆ FilterHistoryPoints()

Nektar::SolverUtils::FilterHistoryPoints::FilterHistoryPoints	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const std::weak_ptr< EquationSystem > &	pEquation,
		const ParamMap &	pParams
	)

Definition at line 53 of file FilterHistoryPoints.cpp.

    : Filter(pSession, pEquation)
{
    // OutputFile
    auto it = pParams.find("OutputFile");
    if (it == pParams.end())
    {
        m_outputFile = m_session->GetSessionName();
    }
    else
    {
        ASSERTL0(it->second.length() > 0, "Missing parameter 'OutputFile'.");
        m_outputFile = it->second;
    }
    if (m_outputFile.length() >= 4 &&
        m_outputFile.substr(m_outputFile.length() - 4) == ".his")
    {
        m_outputFile = m_outputFile.substr(0, m_outputFile.length() - 4);
    }
 
    // OutputFrequency
    it = pParams.find("OutputFrequency");
    if (it == pParams.end())
    {
        m_outputFrequency = 1;
    }
    else
    {
        LibUtilities::Equation equ(m_session->GetInterpreter(), it->second);
        m_outputFrequency = round(equ.Evaluate());
    }
 
    // output all data into one file
    it = pParams.find("OutputOneFile");
    if (it == pParams.end())
    {
        m_outputOneFile = true;
    }
    else
    {
        std::string sOption = it->second.c_str();
        m_outputOneFile     = (boost::iequals(sOption, "true")) ||
                          (boost::iequals(sOption, "yes"));
    }
 
    // OutputPlane
    m_session->MatchSolverInfo("Homogeneous", "1D", m_isHomogeneous1D, false);
    if (m_isHomogeneous1D)
    {
        it = pParams.find("OutputPlane");
        if (it == pParams.end())
        {
            m_outputPlane = -1;
        }
        else
        {
            LibUtilities::Equation equ(m_session->GetInterpreter(), it->second);
            m_outputPlane = round(equ.Evaluate());
        }
 
        it = pParams.find("WaveSpace");
        if (it == pParams.end())
        {
            m_waveSpace = false;
        }
        else
        {
            std::string sOption = it->second.c_str();
            m_waveSpace         = (boost::iequals(sOption, "true")) ||
                          (boost::iequals(sOption, "yes"));
        }
    }
 
    // Points
    if (pParams.end() != (it = pParams.find("Points")))
    {
        m_pointNumMap[LibUtilities::ePtsFile] = Array<OneD, int>(1, 0);
        m_historyPointStream.str(it->second);
    }
    else if (pParams.end() != (it = pParams.find("line")))
    {
        vector<NekDouble> values;
        ASSERTL0(ParseUtils::GenerateVector(it->second, values),
                 "Failed to interpret line string");
 
        ASSERTL0(values.size() > 2, "line string should contain 2*Dim+1 values "
                                    "N,x0,y0,z0,x1,y1,z1");
 
        double tmp;
        ASSERTL0(std::modf(values[0], &tmp) == 0.0, "N is not an integer");
        ASSERTL0(values[0] > 1, "N is not a valid number");
 
        int dim  = (values.size() - 1) / 2;
        int npts = values[0];
 
        Array<OneD, int> num(1, npts);
        Array<OneD, NekDouble> delta(6, 0.);
        for (int i = 0; i < dim; ++i)
        {
            delta[i + 0] = values[i + 1];
            delta[i + 3] = (values[dim + i + 1] - values[i + 1]) / (npts - 1);
        }
        m_pointDatMap[LibUtilities::ePtsLine] = delta;
        m_pointNumMap[LibUtilities::ePtsLine] = num;
    }
    else if (pParams.end() != (it = pParams.find("plane")))
    {
        vector<NekDouble> values;
        ASSERTL0(ParseUtils::GenerateVector(it->second, values),
                 "Failed to interpret plane string");
 
        ASSERTL0(values.size() > 9,
                 "plane string should contain 4 Dim+2 values "
                 "N1,N2,x0,y0,z0,x1,y1,z1,x2,y2,z2,x3,y3,z3");
 
        double tmp;
        ASSERTL0(std::modf(values[0], &tmp) == 0.0, "N1 is not an integer");
        ASSERTL0(std::modf(values[1], &tmp) == 0.0, "N2 is not an integer");
 
        ASSERTL0(values[0] > 1, "N1 is not a valid number");
        ASSERTL0(values[1] > 1, "N2 is not a valid number");
 
        int dim = (values.size() - 2) / 4;
 
        Array<OneD, int> npts(3);
        npts[0] = values[0];
        npts[1] = values[0] * values[1];
        npts[2] = values[1];
 
        Array<OneD, NekDouble> delta(12, 0.);
        for (int i = 0; i < dim; ++i)
        {
            delta[i + 0] = values[2 + i];
            delta[i + 3] = values[2 + 3 * dim + i];
            delta[i + 6] = (values[2 + 1 * dim + i] - values[2 + 0 * dim + i]) /
                           (values[0] - 1);
            delta[i + 9] = (values[2 + 2 * dim + i] - values[2 + 3 * dim + i]) /
                           (values[0] - 1);
        }
        m_pointDatMap[LibUtilities::ePtsPlane] = delta;
        m_pointNumMap[LibUtilities::ePtsPlane] = npts;
    }
    else if (pParams.end() != (it = pParams.find("box")))
    {
        vector<NekDouble> values;
        ASSERTL0(ParseUtils::GenerateVector(it->second, values),
                 "Failed to interpret box string");
 
        ASSERTL0(values.size() == 9, "box string should contain 9 values "
                                     "N1,N2,N3,xmin,xmax,ymin,ymax,zmin,zmax");
 
        int dim = 3;
        Array<OneD, int> npts(3);
        npts[0] = values[0];
        npts[1] = values[0] * values[1];
        npts[2] = values[0] * values[1] * values[2];
 
        Array<OneD, NekDouble> delta(6, 0.);
        for (int i = 0; i < dim; ++i)
        {
            delta[i + 0] = values[3 + 2 * i];
            delta[i + 3] =
                (values[4 + 2 * i] - values[3 + 2 * i]) / (values[i] - 1);
        }
        m_pointDatMap[LibUtilities::ePtsBox] = delta;
        m_pointNumMap[LibUtilities::ePtsBox] = npts;
    }
    else
    {
        ASSERTL0(false, "Missing parameter 'Points'.");
    }
}

References ASSERTL0, Nektar::LibUtilities::ePtsBox, Nektar::LibUtilities::ePtsFile, Nektar::LibUtilities::ePtsLine, Nektar::LibUtilities::ePtsPlane, Nektar::LibUtilities::Equation::Evaluate(), Nektar::ParseUtils::GenerateVector(), m_historyPointStream, m_isHomogeneous1D, m_outputFile, m_outputFrequency, m_outputOneFile, m_outputPlane, m_pointDatMap, m_pointNumMap, Nektar::SolverUtils::Filter::m_session, and m_waveSpace.

◆ ~FilterHistoryPoints()

Nektar::SolverUtils::FilterHistoryPoints::~FilterHistoryPoints ( )

override

Definition at line 231 of file FilterHistoryPoints.cpp.

232{

233}

Member Function Documentation

◆ create()

static FilterSharedPtr Nektar::SolverUtils::FilterHistoryPoints::create	(	const LibUtilities::SessionReaderSharedPtr &	pSession,
		const std::weak_ptr< EquationSystem > &	pEquation,
		const std::map< std::string, std::string > &	pParams
	)

inlinestatic

Creates an instance of this class.

Definition at line 51 of file FilterHistoryPoints.h.

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

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and CellMLToNektar.cellml_metadata::p.

◆ GetPoint()

bool Nektar::SolverUtils::FilterHistoryPoints::GetPoint	(	Array< OneD, NekDouble >	gloCoord,
		int	I
	)

protected

Definition at line 235 of file FilterHistoryPoints.cpp.

{
    if (m_pointNumMap.count(LibUtilities::ePtsFile))
    {
        m_historyPointStream >> gloCoord[0] >> gloCoord[1] >> gloCoord[2];
        return !m_historyPointStream.fail();
    }
    else if (m_pointNumMap.count(LibUtilities::ePtsLine))
    {
        if (I >= m_pointNumMap[LibUtilities::ePtsLine][0])
        {
            return false;
        }
        Array<OneD, NekDouble> values = m_pointDatMap[LibUtilities::ePtsLine];
        Array<OneD, NekDouble> delta =
            m_pointDatMap[LibUtilities::ePtsLine] + 3;
        for (int n = 0; n < 3; ++n)
        {
            gloCoord[n] = values[n] + I * delta[n];
        }
        return true;
    }
    else if (m_pointNumMap.count(LibUtilities::ePtsPlane))
    {
        if (I >= m_pointNumMap[LibUtilities::ePtsPlane][1])
        {
            return false;
        }
        Array<OneD, NekDouble> values = m_pointDatMap[LibUtilities::ePtsPlane];
        Array<OneD, NekDouble> delta =
            m_pointDatMap[LibUtilities::ePtsPlane] + 6;
        Array<OneD, int> i(2);
        i[1]      = I / m_pointNumMap[LibUtilities::ePtsPlane][0];
        i[0]      = I - i[1] * m_pointNumMap[LibUtilities::ePtsPlane][0];
        double n1 = -1. + (NekDouble)m_pointNumMap[LibUtilities::ePtsPlane][2];
        for (int n = 0; n < 3; ++n)
        {
            gloCoord[n] = (values[n] + i[0] * delta[n]) * (1.0 - i[1] / n1) +
                          (values[3 + n] + i[0] * delta[3 + n]) * (i[1] / n1);
        }
        return true;
    }
    else if (m_pointNumMap.count(LibUtilities::ePtsBox))
    {
        if (I >= m_pointNumMap[LibUtilities::ePtsBox][2])
        {
            return false;
        }
        Array<OneD, NekDouble> values = m_pointDatMap[LibUtilities::ePtsBox];
        Array<OneD, NekDouble> delta = m_pointDatMap[LibUtilities::ePtsBox] + 3;
        Array<OneD, int> i(3);
        i[2] = I / m_pointNumMap[LibUtilities::ePtsBox][1];
        I -= i[2] * m_pointNumMap[LibUtilities::ePtsBox][1];
        i[1] = I / m_pointNumMap[LibUtilities::ePtsBox][0];
        i[0] = I - i[1] * m_pointNumMap[LibUtilities::ePtsBox][0];
        for (int n = 0; n < 3; ++n)
        {
            gloCoord[n] = values[n] + i[n] * delta[n];
        }
        return true;
    }
    return false;
}

References Nektar::LibUtilities::ePtsBox, Nektar::LibUtilities::ePtsFile, Nektar::LibUtilities::ePtsLine, Nektar::LibUtilities::ePtsPlane, m_historyPointStream, m_pointDatMap, and m_pointNumMap.

Referenced by v_Initialise().

◆ v_Finalise()

void Nektar::SolverUtils::FilterHistoryPoints::v_Finalise	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		const NekDouble &	time
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::Filter.

Definition at line 685 of file FilterHistoryPoints.cpp.

{
    if (pFields[0]->GetComm()->GetRank() == 0 && m_outputOneFile)
    {
        m_outputStream.close();
    }
}

References m_outputOneFile, and m_outputStream.

◆ v_Initialise()

void Nektar::SolverUtils::FilterHistoryPoints::v_Initialise	(	const Array< OneD, const MultiRegions::ExpListSharedPtr > &	pFields,
		const NekDouble &	time
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::Filter.

Definition at line 302 of file FilterHistoryPoints.cpp.

{
    ASSERTL0(!m_historyPointStream.fail(), "No history points in stream.");
 
    LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
 
    // All processes read all the history points
    Array<OneD, NekDouble> gloCoord(3, 0.0);
    std::vector<Array<OneD, const NekDouble>> GloCoords;
    const int coordim = pFields[0]->GetGraph()->GetSpaceDimension();
 
    int i = -1;
    while (GetPoint(gloCoord, ++i))
    {
        // Overwrite gloCoord[2] for 3DH1D using m_outputPlane
        // if it is defined
        if (m_isHomogeneous1D && (m_outputPlane != -1 || m_waveSpace))
        {
            int nplanes    = pFields[0]->GetHomogeneousBasis()->GetZ().size();
            NekDouble lhom = pFields[0]->GetHomoLen();
            int plane;
            if (m_outputPlane != -1)
            {
                plane = m_outputPlane;
            }
            else
            {
                // Pick plane near the point
                plane = round((gloCoord[coordim] * nplanes) / lhom);
            }
            if (m_waveSpace)
            {
                m_planeIDs.push_back(plane);
            }
            NekDouble Z = (pFields[0]->GetHomogeneousBasis()->GetZ())[plane];
            Z           = (Z + 1) * lhom / 2;
            if (fabs(gloCoord[coordim] - Z) >
                    NekConstants::kVertexTheSameDouble &&
                vComm->GetRank() == 0)
            {
                cout << "Resetting History point from z = " << gloCoord[coordim]
                     << " to z = " << Z << endl;
            }
            gloCoord[coordim] = Z;
        }
 
        // Save a copy of the global coordinate
        GloCoords.push_back(
            Array<OneD, const NekDouble>(gloCoord.size(), gloCoord.data()));
    }
 
    // Check which history points this process is responsible for
    const int vRank       = vComm->GetRank();
    const int vColumnRank = vComm->GetColumnComm()->GetRank();
    const int vRowRank    = vComm->GetRowComm()->GetRank();
    m_historyPointsSize   = GloCoords.size();
    Array<OneD, int> procList(m_historyPointsSize, -1);
    Array<OneD, int> idList(m_historyPointsSize, -1);
 
    Array<OneD, NekDouble> locCoord(pFields[0]->GetShapeDimension(), 0.);
    std::vector<Array<OneD, const NekDouble>> LocCoords;
 
    // For each history point, find the element that contains it. If no element
    // is available on this process, idList[i] = -1 and procList[i] = -1
    for (i = 0; i < m_historyPointsSize; ++i)
    {
        // Create a copy of the global coordinate
        gloCoord = GloCoords[i];
 
        // Determine the expansion and local coordinates
        if (m_isHomogeneous1D)
        {
            idList[i] = pFields[0]->GetPlane(0)->GetExpIndex(
                gloCoord, locCoord, NekConstants::kGeomFactorsTol);
        }
        else
        {
            idList[i] = pFields[0]->GetExpIndex(gloCoord, locCoord,
                                                NekConstants::kGeomFactorsTol);
        }
 
        for (int j = 0; j < locCoord.size(); ++j)
        {
            locCoord[j] = std::max(locCoord[j], -1.0);
            locCoord[j] = std::min(locCoord[j], 1.0);
        }
 
        // If this process owns an element that contains the history point, we
        // flag this by updating the corresponding element in procList
        if (idList[i] != -1 && vColumnRank == 0)
        {
            procList[i] = vRank;
        }
 
        // Save a copy of the local coordinates
        LocCoords.push_back(
            Array<OneD, const NekDouble>(locCoord.size(), locCoord.data()));
    }
 
    // Clear list before populating it
    m_historyList.clear();
 
    // Share process ranks among all processes. If multiple processes found
    // that they own a history point, procList will only contain the rank of
    // the process with the largest rank
    vComm->AllReduce(procList, LibUtilities::ReduceMax);
    for (i = 0; i < m_historyPointsSize; ++i)
    {
        // If process owns a history point, it saves a reference to its global
        // and local coordinates. Note that, if multiple processes have found
        // that they own a history point, which can happen if the point is
        // located on a partition boundary, only the process with the highest
        // rank will store the point.
        if (procList[i] == vRowRank)
        {
            // Note that the reference counting in the Array class prevents the
            // data stored in GloCoords and LocCoords from being deallocated
            // when they go out of scope, since m_historyList will increase the
            // reference counter by 1 when it gets a reference to the data.
            m_historyList.push_back(
                std::make_tuple(GloCoords[i], LocCoords[i], i, idList[i]));
        }
    }
 
    // Root process must have access to all coordinates for writing the file
    if (vComm->TreatAsRankZero())
    {
        m_historyPoints =
            Array<OneD, Array<OneD, const NekDouble>>(m_historyPointsSize);
 
        for (i = 0; i < m_historyPointsSize; ++i)
        {
            // Create a reference to the underlying array. Note that the data
            // stored in GloCoords won't be deallocated because the Array class
            // has a reference counter.
            m_historyPoints[i] = GloCoords[i];
        }
    }
 
    // Check that each history point is allocated to a process
    if (vComm->TreatAsRankZero())
    {
        for (i = 0; i < m_historyPointsSize; ++i)
        {
            gloCoord = m_historyPoints[i];
 
            // Write an error if no process owns history point
            ASSERTL0(procList[i] != -1,
                     "History point " +
                         boost::lexical_cast<std::string>(gloCoord[0]) + ", " +
                         boost::lexical_cast<std::string>(gloCoord[1]) + ", " +
                         boost::lexical_cast<std::string>(gloCoord[2]) +
                         " cannot be found in the mesh.");
        }
 
        m_session->MatchSolverInfo("Driver", "Adaptive", m_adaptive, false);
    }
    v_Update(pFields, time);
}

References ASSERTL0, GetPoint(), Nektar::NekConstants::kGeomFactorsTol, Nektar::NekConstants::kVertexTheSameDouble, m_adaptive, m_historyList, m_historyPoints, m_historyPointsSize, m_historyPointStream, m_isHomogeneous1D, m_outputPlane, m_planeIDs, Nektar::SolverUtils::Filter::m_session, m_waveSpace, Nektar::LibUtilities::ReduceMax, and v_Update().

◆ v_IsTimeDependent()

bool Nektar::SolverUtils::FilterHistoryPoints::v_IsTimeDependent ( )

overrideprotectedvirtual

Implements Nektar::SolverUtils::Filter.

Definition at line 698 of file FilterHistoryPoints.cpp.

{
    return true;
}

◆ v_Update()

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

overrideprotectedvirtual

Implements Nektar::SolverUtils::Filter.

Reimplemented in Nektar::FilterCellHistoryPoints.

Definition at line 467 of file FilterHistoryPoints.cpp.

{
    // Only output every m_outputFrequency.
    if ((m_index++) % m_outputFrequency)
    {
        return;
    }
 
    const int numFields = pFields.size();
    const int coordim   = pFields[0]->GetGraph()->GetSpaceDimension();
    LibUtilities::CommSharedPtr vComm = pFields[0]->GetComm();
    Array<OneD, NekDouble> data(m_historyPointsSize * numFields, 0.0);
    Array<OneD, NekDouble> physvals;
    Array<OneD, NekDouble> locCoord;
 
    // Pull out data values field by field
    Array<OneD, NekDouble> gloCoord(3, 0.0);
    for (int j = 0; j < numFields; ++j)
    {
        if (m_isHomogeneous1D)
        {
            Array<OneD, const unsigned int> planes = pFields[j]->GetZIDs();
            int nPlanes    = pFields[j]->GetHomogeneousBasis()->GetZ().size();
            NekDouble lHom = pFields[j]->GetHomoLen();
            for (auto &x : m_historyList)
            {
                ASSERTL0(pFields[j]->GetWaveSpace(),
                         "HistoryPoints in Homogeneous1D require that solution "
                         "is in wavespace");
                gloCoord        = std::get<0>(x);
                locCoord        = std::get<1>(x);
                const int indx  = std::get<2>(x);
                const int expId = std::get<3>(x);
 
                NekDouble value = 0.0;
                NekDouble BetaT =
                    2. * M_PI * fmod(gloCoord[coordim], lHom) / lHom;
                for (size_t n = 0; n < planes.size(); ++n)
                {
                    if (m_waveSpace && planes[n] != m_planeIDs[indx])
                    {
                        continue;
                    }
                    physvals = pFields[j]->GetPlane(n)->UpdatePhys() +
                               pFields[j]->GetPhys_Offset(expId);
 
                    // transform elemental data if required.
                    if (pFields[j]->GetPhysState() == false)
                    {
                        pFields[j]->GetPlane(n)->GetExp(expId)->BwdTrans(
                            pFields[j]->GetPlane(n)->GetCoeffs() +
                                pFields[j]->GetCoeff_Offset(expId),
                            physvals);
                    }
                    // Interpolate data
                    NekDouble coeff =
                        pFields[j]->GetPlane(n)->GetExp(expId)->StdPhysEvaluate(
                            locCoord, physvals);
 
                    if (m_waveSpace)
                    {
                        value = coeff;
                    }
                    else
                    {
                        if (planes[n] == 0)
                        {
                            value += coeff;
                        }
                        else if (planes[n] == 1)
                        {
                            value += cos(0.5 * nPlanes * BetaT) * coeff;
                        }
                        else if (planes[n] % 2 == 0)
                        {
                            NekDouble phase = (planes[n] >> 1) * BetaT;
                            value += cos(phase) * coeff;
                        }
                        else
                        {
                            NekDouble phase = (planes[n] >> 1) * BetaT;
                            value += -sin(phase) * coeff;
                        }
                    }
                }
                // store data
                data[indx * numFields + j] = value;
            }
        }
        else
        {
            for (auto &x : m_historyList)
            {
                locCoord        = std::get<1>(x);
                const int indx  = std::get<2>(x);
                const int expId = std::get<3>(x);
 
                physvals = pFields[j]->UpdatePhys() +
                           pFields[j]->GetPhys_Offset(expId);
 
                // transform elemental data if required.
                if (pFields[j]->GetPhysState() == false)
                {
                    pFields[j]->GetExp(expId)->BwdTrans(
                        pFields[j]->GetCoeffs() +
                            pFields[j]->GetCoeff_Offset(expId),
                        physvals);
                }
 
                // interpolate point
                data[indx * numFields + j] =
                    pFields[j]->GetExp(expId)->StdPhysEvaluate(locCoord,
                                                               physvals);
            }
        }
    }
 
    // Exchange history data
    // This could be improved to reduce communication but works for now
    vComm->AllReduce(data, LibUtilities::ReduceSum);
 
    // TODO: Why not only call this routine if we are rank 0?
    v_WriteData(vComm->GetRank(), data, numFields, time);
}

References ASSERTL0, m_historyList, m_historyPointsSize, m_index, m_isHomogeneous1D, m_outputFrequency, m_planeIDs, m_waveSpace, Nektar::LibUtilities::ReduceSum, and v_WriteData().

Referenced by v_Initialise().

◆ v_WriteData()

void Nektar::SolverUtils::FilterHistoryPoints::v_WriteData	(	const int &	rank,
		const Array< OneD, NekDouble > &	data,
		const int &	numFields,
		const NekDouble &	time
	)

protected

Definition at line 594 of file FilterHistoryPoints.cpp.

{
    // Only the root process writes out history data
    if (rank == 0)
    {
        Array<OneD, NekDouble> gloCoord(3, 0.0);
        if (!m_outputOneFile || m_index == 1)
        {
            std::stringstream vOutputFilename;
            if (m_outputOneFile)
            {
                vOutputFilename << m_outputFile << ".his";
            }
            else
            {
                vOutputFilename << m_outputFile << "_" << m_outputIndex
                                << ".his";
            }
            ++m_outputIndex;
            if (m_adaptive)
            {
                m_outputStream.open(vOutputFilename.str().c_str(),
                                    ofstream::app);
            }
            else
            {
                m_outputStream.open(vOutputFilename.str().c_str());
            }
            m_outputStream << "# History data for variables (:";
 
            for (int i = 0; i < numFields; ++i)
            {
                m_outputStream << m_session->GetVariable(i) << ",";
            }
 
            if (m_isHomogeneous1D)
            {
                m_outputStream << ") at points:" << endl;
            }
            else
            {
                m_outputStream << ") at points:" << endl;
            }
 
            for (int i = 0; i < m_historyPoints.size(); ++i)
            {
                gloCoord = m_historyPoints[i];
 
                m_outputStream << "# " << boost::format("%6.0f") % i;
                m_outputStream << " " << boost::format("%25.19e") % gloCoord[0];
                m_outputStream << " " << boost::format("%25.19e") % gloCoord[1];
                m_outputStream << " " << boost::format("%25.19e") % gloCoord[2];
                m_outputStream << endl;
            }
 
            if (m_isHomogeneous1D)
            {
                if (m_waveSpace)
                {
                    m_outputStream << "# (in Wavespace)" << endl;
                }
            }
        }
 
        // Write data values point by point
        for (int k = 0; k < m_historyPoints.size(); ++k)
        {
            m_outputStream << boost::format("%25.19e") % time;
            for (int j = 0; j < numFields; ++j)
            {
                m_outputStream
                    << " "
                    << boost::format("%25.19e") % data[k * numFields + j];
            }
            m_outputStream << endl;
        }
 
        if (!m_outputOneFile)
        {
            m_outputStream.close();
            ;
        }
    }
}

References CellMLToNektar.pycml::format, m_adaptive, m_historyPoints, m_index, m_isHomogeneous1D, m_outputFile, m_outputIndex, m_outputOneFile, m_outputStream, Nektar::SolverUtils::Filter::m_session, and m_waveSpace.

Referenced by v_Update().

Friends And Related Function Documentation

◆ MemoryManager< FilterHistoryPoints >

friend class MemoryManager< FilterHistoryPoints >

friend

Definition at line 1 of file FilterHistoryPoints.h.

Member Data Documentation

◆ className

std::string Nektar::SolverUtils::FilterHistoryPoints::className

static

Initial value:

=
    GetFilterFactory().RegisterCreatorFunction("HistoryPoints",
                                               FilterHistoryPoints::create)

Name of the class.

Definition at line 63 of file FilterHistoryPoints.h.

◆ m_adaptive

bool Nektar::SolverUtils::FilterHistoryPoints::m_adaptive

protected

Definition at line 114 of file FilterHistoryPoints.h.

Referenced by v_Initialise(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_historyList

std::list<std::tuple<Array<OneD, const NekDouble>, Array<OneD, const NekDouble>, int, int> > Nektar::SolverUtils::FilterHistoryPoints::m_historyList

protected

Definition at line 109 of file FilterHistoryPoints.h.

Referenced by v_Initialise(), v_Update(), and Nektar::FilterCellHistoryPoints::v_Update().

◆ m_historyPoints

Array<OneD, Array<OneD, const NekDouble> > Nektar::SolverUtils::FilterHistoryPoints::m_historyPoints

protected

Initial value:

=

Array<OneD, Array<OneD, const NekDouble>>(0)

Definition at line 88 of file FilterHistoryPoints.h.

Referenced by v_Initialise(), Nektar::FilterCellHistoryPoints::v_Update(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_historyPointsSize

size_t Nektar::SolverUtils::FilterHistoryPoints::m_historyPointsSize = 0

protected

Definition at line 90 of file FilterHistoryPoints.h.

Referenced by v_Initialise(), and v_Update().

◆ m_historyPointStream

std::stringstream Nektar::SolverUtils::FilterHistoryPoints::m_historyPointStream

protected

Definition at line 100 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), GetPoint(), and v_Initialise().

◆ m_index

unsigned int Nektar::SolverUtils::FilterHistoryPoints::m_index = 0

protected

Definition at line 91 of file FilterHistoryPoints.h.

Referenced by v_Update(), Nektar::FilterCellHistoryPoints::v_Update(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_isHomogeneous1D

bool Nektar::SolverUtils::FilterHistoryPoints::m_isHomogeneous1D

protected

Definition at line 96 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_Initialise(), v_Update(), Nektar::FilterCellHistoryPoints::v_Update(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_outputFile

std::string Nektar::SolverUtils::FilterHistoryPoints::m_outputFile

protected

Definition at line 98 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_outputFrequency

unsigned int Nektar::SolverUtils::FilterHistoryPoints::m_outputFrequency

protected

Definition at line 92 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_Update(), and Nektar::FilterCellHistoryPoints::v_Update().

◆ m_outputIndex

unsigned int Nektar::SolverUtils::FilterHistoryPoints::m_outputIndex = 0

protected

Definition at line 112 of file FilterHistoryPoints.h.

Referenced by v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_outputOneFile

bool Nektar::SolverUtils::FilterHistoryPoints::m_outputOneFile

protected

Definition at line 113 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_Finalise(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_outputPlane

int Nektar::SolverUtils::FilterHistoryPoints::m_outputPlane

protected

plane to take history point from if using a homogeneous1D expansion

Definition at line 94 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_Initialise(), and Nektar::FilterCellHistoryPoints::v_Update().

◆ m_outputStream

std::ofstream Nektar::SolverUtils::FilterHistoryPoints::m_outputStream

protected

Definition at line 99 of file FilterHistoryPoints.h.

Referenced by v_Finalise(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

◆ m_planeIDs

std::vector<int> Nektar::SolverUtils::FilterHistoryPoints::m_planeIDs

protected

Definition at line 95 of file FilterHistoryPoints.h.

Referenced by v_Initialise(), and v_Update().

◆ m_pointDatMap

std::map<LibUtilities::PtsType, Array<OneD, NekDouble> > Nektar::SolverUtils::FilterHistoryPoints::m_pointDatMap

protected

Definition at line 110 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), and GetPoint().

◆ m_pointNumMap

std::map<LibUtilities::PtsType, Array<OneD, int> > Nektar::SolverUtils::FilterHistoryPoints::m_pointNumMap

protected

Definition at line 111 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), and GetPoint().

◆ m_waveSpace

bool Nektar::SolverUtils::FilterHistoryPoints::m_waveSpace

protected

Definition at line 97 of file FilterHistoryPoints.h.

Referenced by FilterHistoryPoints(), v_Initialise(), v_Update(), v_WriteData(), and Nektar::FilterCellHistoryPoints::v_WriteData().

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Friends

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ FilterHistoryPoints()

◆ ~FilterHistoryPoints()

Member Function Documentation

◆ create()

◆ GetPoint()

◆ v_Finalise()

◆ v_Initialise()

◆ v_IsTimeDependent()

◆ v_Update()

◆ v_WriteData()

Friends And Related Function Documentation

◆ MemoryManager< FilterHistoryPoints >

Member Data Documentation

◆ className

◆ m_adaptive

◆ m_historyList

◆ m_historyPoints

◆ m_historyPointsSize

◆ m_historyPointStream

◆ m_index

◆ m_isHomogeneous1D

◆ m_outputFile

◆ m_outputFrequency

◆ m_outputIndex

◆ m_outputOneFile

◆ m_outputPlane

◆ m_outputStream

◆ m_planeIDs

◆ m_pointDatMap

◆ m_pointNumMap

◆ m_waveSpace