#include <FilterLagrangianPoints.h>

Inheritance diagram for Nektar::SolverUtils::StatLagrangianPoints:

Public Member Functions
	StatLagrangianPoints (int rank, int dim, int intOrder, int idOffset, NekDouble dt, const std::vector< int > &Np, const std::vector< NekDouble > &Box, std::vector< std::string > extraVars)

	~StatLagrangianPoints () override

Public Member Functions inherited from Nektar::SolverUtils::StationaryPoints
virtual	~StationaryPoints ()=default

SOLVER_UTILS_EXPORT void	OutputData (std::string filename, bool verbose, std::map< std::string, NekDouble > &params)

void	ReSize (int Np)

void	AssignPoint (int id, int pid, const Array< OneD, NekDouble > &gcoords)

void	GetCoordsByPID (int pid, Array< OneD, NekDouble > &gcoords)

void	SetCoordsByPID (int pid, const Array< OneD, NekDouble > &gcoords)

void	GetPhysicsByPID (int pid, Array< OneD, NekDouble > &data)

void	SetPhysicsByPID (int pid, const Array< OneD, NekDouble > &data)

void	GetCoords (int pid, Array< OneD, NekDouble > &gcoords)

void	SetCoords (int pid, const Array< OneD, NekDouble > &gcoords)

void	GetPhysics (int pid, Array< OneD, NekDouble > &data)

void	SetPhysics (int pid, const Array< OneD, NekDouble > &data)

void	TimeAdvance (int order)

int	GetTotPoints ()

int	GetDim ()

int	LocalToGlobal (int id)

int	GlobalToLocal (int id)

Static Public Member Functions
static StationaryPointsSharedPtr	create (int rank, int dim, int intOrder, int idOffset, NekDouble dt, const std::vector< int > &Np, const std::vector< NekDouble > &Box, std::vector< std::string > extraVars)
	Creates an instance of this class. More...

Protected Member Functions
void	v_OutputData (std::string filename, bool verbose, std::map< std::string, NekDouble > &params) override

void	v_TimeAdvance (int order) override

void	v_GetCoords (int pid, Array< OneD, NekDouble > &gcoords) override

void	v_SetCoords (int pid, const Array< OneD, NekDouble > &gcoords) override

void	v_GetPhysics (int pid, Array< OneD, NekDouble > &data) override

void	v_SetPhysics (int pid, const Array< OneD, NekDouble > &data) override

void	v_ReSize (int Np) override

void	v_AssignPoint (int id, int pid, const Array< OneD, NekDouble > &gcoords) override

virtual void	v_OutputData (std::string filename, bool verbose, std::map< std::string, NekDouble > &params)=0

virtual void	v_TimeAdvance (int order)=0

virtual void	v_GetCoords (int pid, Array< OneD, NekDouble > &gcoords)=0

virtual void	v_SetCoords (int pid, const Array< OneD, NekDouble > &gcoords)=0

virtual void	v_GetPhysics (int pid, Array< OneD, NekDouble > &data)=0

virtual void	v_SetPhysics (int pid, const Array< OneD, NekDouble > &data)=0

virtual void	v_ReSize (int Np)=0

virtual void	v_AssignPoint (int id, int pid, const Array< OneD, NekDouble > &gcoords)

Protected Attributes
Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_coords

Array< OneD, Array< OneD, Array< OneD, NekDouble > > >	m_velocity

Array< OneD, Array< OneD, NekDouble > >	m_extraPhysics

std::vector< std::string >	m_extraPhysVars

NekDouble	m_dt

int	m_idOffset

std::vector< int >	m_N

int	m_intOrder

Protected Attributes inherited from Nektar::SolverUtils::StationaryPoints
int	m_dim

int	m_totPts

std::map< int, int >	m_localIDToGlobal

std::map< int, int >	m_globalIDToLocal

Friends
class	MemoryManager< StatLagrangianPoints >

Detailed Description

Definition at line 49 of file FilterLagrangianPoints.h.

Constructor & Destructor Documentation

◆ StatLagrangianPoints()

Nektar::SolverUtils::StatLagrangianPoints::StatLagrangianPoints	(	int	rank,
		int	dim,
		int	intOrder,
		int	idOffset,
		NekDouble	dt,
		const std::vector< int > &	Np,
		const std::vector< NekDouble > &	Box,
		std::vector< std::string >	extraVars
	)

Definition at line 366 of file FilterLagrangianPoints.cpp.

{
    m_idOffset = idOffset;
    m_dim      = dim;
    if (intOrder < 0)
    {
        m_intOrder = -intOrder;
        m_dt       = -dt;
    }
    else
    {
        m_intOrder = intOrder;
        m_dt       = dt;
    }
    if (Np.size() < m_dim || Box.size() < m_dim * 2)
    {
        throw ErrorUtil::NekError("Not enough input for StationaryPoints.");
    }
    Array<OneD, NekDouble> delta(3, 1.);
    m_N.resize(3, 1);
    m_totPts = 1;
    for (int d = 0; d < m_dim; ++d)
    {
        m_N[d] = Np[d];
        m_totPts *= Np[d];
        delta[d] = Np[d] < 2 ? 1. : (Box[d * 2 + 1] - Box[d * 2]) / (Np[d] - 1);
    }
    m_coords   = Array<OneD, Array<OneD, Array<OneD, NekDouble>>>(m_intOrder);
    m_velocity = Array<OneD, Array<OneD, Array<OneD, NekDouble>>>(m_intOrder);
    for (int i = 0; i < m_intOrder; ++i)
    {
        m_coords[i]   = Array<OneD, Array<OneD, NekDouble>>(m_dim);
        m_velocity[i] = Array<OneD, Array<OneD, NekDouble>>(m_dim);
        for (int d = 0; d < m_dim; ++d)
        {
            m_coords[i][d]   = Array<OneD, NekDouble>(m_totPts, 0.);
            m_velocity[i][d] = Array<OneD, NekDouble>(m_totPts, 0.);
        }
    }
    m_extraPhysics = Array<OneD, Array<OneD, NekDouble>>(extraVars.size());
    for (size_t i = 0; i < extraVars.size(); ++i)
    {
        m_extraPhysVars.push_back(extraVars[i]);
        m_extraPhysics[i] = Array<OneD, NekDouble>(m_totPts, 0.);
    }
    // initialise points
    int count = 0;
    for (int k = 0; k < m_N[2]; ++k)
    {
        for (int j = 0; j < m_N[1]; ++j)
        {
            for (int i = 0; i < m_N[0]; ++i)
            {
                if (m_dim > 0)
                {
                    m_coords[0][0][count] = Box[0] + delta[0] * i;
                }
                if (m_dim > 1)
                {
                    m_coords[0][1][count] = Box[2] + delta[1] * j;
                }
                if (m_dim > 2)
                {
                    m_coords[0][2][count] = Box[4] + delta[2] * k;
                }
                m_localIDToGlobal[count]              = m_idOffset + count;
                m_globalIDToLocal[m_idOffset + count] = count;
                ++count;
            }
        }
    }
}

References Nektar::UnitTests::d(), m_coords, Nektar::SolverUtils::StationaryPoints::m_dim, m_dt, m_extraPhysics, m_extraPhysVars, Nektar::SolverUtils::StationaryPoints::m_globalIDToLocal, m_idOffset, m_intOrder, Nektar::SolverUtils::StationaryPoints::m_localIDToGlobal, m_N, Nektar::SolverUtils::StationaryPoints::m_totPts, and m_velocity.

◆ ~StatLagrangianPoints()

Nektar::SolverUtils::StatLagrangianPoints::~StatLagrangianPoints ( )

inlineoverride

Definition at line 72 of file FilterLagrangianPoints.h.

73 {

74 }

Member Function Documentation

◆ create()

static StationaryPointsSharedPtr Nektar::SolverUtils::StatLagrangianPoints::create	(	int	rank,
		int	dim,
		int	intOrder,
		int	idOffset,
		NekDouble	dt,
		const std::vector< int > &	Np,
		const std::vector< NekDouble > &	Box,
		std::vector< std::string >	extraVars
	)

inlinestatic

Creates an instance of this class.

Definition at line 55 of file FilterLagrangianPoints.h.

    {
        StationaryPointsSharedPtr p =
            MemoryManager<StatLagrangianPoints>::AllocateSharedPtr(
                rank, dim, intOrder, idOffset, dt, Np, Box, extraVars);
        return p;
    }

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

◆ v_AssignPoint()

void Nektar::SolverUtils::StatLagrangianPoints::v_AssignPoint	(	int	id,
		int	pid,
		const Array< OneD, NekDouble > &	gcoords
	)

overrideprotectedvirtual

Reimplemented from Nektar::SolverUtils::StationaryPoints.

Definition at line 356 of file FilterLagrangianPoints.cpp.

{
    StationaryPoints::v_AssignPoint(id, pid, gcoords);
    for (int d = 0; d < m_dim; ++d)
    {
        m_coords[0][d][id] = gcoords[d];
    }
}

References Nektar::UnitTests::d(), m_coords, Nektar::SolverUtils::StationaryPoints::m_dim, and Nektar::SolverUtils::StationaryPoints::v_AssignPoint().

◆ v_GetCoords()

void Nektar::SolverUtils::StatLagrangianPoints::v_GetCoords	(	int	pid,
		Array< OneD, NekDouble > &	gcoords
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 51 of file FilterLagrangianPoints.cpp.

{
    for (int d = 0; d < m_dim; ++d)
    {
        gcoords[d] = m_coords[0][d][id];
    }
}

References Nektar::UnitTests::d(), m_coords, and Nektar::SolverUtils::StationaryPoints::m_dim.

◆ v_GetPhysics()

void Nektar::SolverUtils::StatLagrangianPoints::v_GetPhysics	(	int	pid,
		Array< OneD, NekDouble > &	data
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 68 of file FilterLagrangianPoints.cpp.

{
    for (int d = 0; d < m_dim; ++d)
    {
        data[d] = m_velocity[0][d][id];
    }
}

References Nektar::UnitTests::d(), Nektar::SolverUtils::StationaryPoints::m_dim, and m_velocity.

◆ v_OutputData()

void Nektar::SolverUtils::StatLagrangianPoints::v_OutputData	(	std::string	filename,
		bool	verbose,
		std::map< std::string, NekDouble > &	params
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 279 of file FilterLagrangianPoints.cpp.

{
    std::vector<std::string> COORDS = {"x", "y", "z"};
    std::vector<std::string> VELOCI = {"u", "v", "w"};
    std::vector<std::string> variables;
    std::vector<Array<OneD, NekDouble>> data;
    for (int d = 0; d < m_dim; ++d)
    {
        variables.push_back(COORDS[d]);
    }
    for (int d = 0; d < m_dim; ++d)
    {
        variables.push_back(VELOCI[d]);
    }
    if (params.size() > 0)
    {
        for (int d = 0; d < m_dim; ++d)
        {
            Array<OneD, NekDouble> x(m_totPts, 0.);
            Vmath::Sadd(m_totPts, params[COORDS[d]], m_coords[0][d], 1, x, 1);
            data.push_back(x);
        }
        for (int d = 0; d < m_dim; ++d)
        {
            Array<OneD, NekDouble> x(m_totPts, 0.);
            Vmath::Sadd(m_totPts, params[VELOCI[d]], m_velocity[0][d], 1, x, 1);
            data.push_back(x);
        }
    }
    else
    {
        for (int d = 0; d < m_dim; ++d)
        {
            data.push_back(m_coords[0][d]);
        }
        for (int d = 0; d < m_dim; ++d)
        {
            data.push_back(m_velocity[0][d]);
        }
    }
    for (size_t d = 0; d < m_extraPhysics.size(); ++d)
    {
        data.push_back(m_extraPhysics[d]);
        variables.push_back(m_extraPhysVars[d]);
    }
    OutputTec360_binary(filename, variables, m_N, data, 1);
    if (verbose)
    {
        int Np = m_coords[0][0].size();
        for (int d = 0; d < m_dim; ++d)
        {
            NekDouble value = 0.;
            for (int i = 0; i < Np; ++i)
            {
                value += m_coords[0][d][i] * m_coords[0][d][i];
            }
            value = sqrt(value / Np);
            cout << "L 2 error (variable L" << COORDS[d] << ") : " << value
                 << endl;
        }
        for (int d = 0; d < m_dim; ++d)
        {
            NekDouble value = 0.;
            for (int i = 0; i < Np; ++i)
            {
                value += m_velocity[0][d][i] * m_velocity[0][d][i];
            }
            value = sqrt(value / Np);
            cout << "L 2 error (variable L" << VELOCI[d] << ") : " << value
                 << endl;
        }
    }
    cout << "Write file " << filename << endl;
}

References Nektar::UnitTests::d(), m_coords, Nektar::SolverUtils::StationaryPoints::m_dim, m_extraPhysics, m_extraPhysVars, m_N, Nektar::SolverUtils::StationaryPoints::m_totPts, m_velocity, Nektar::SolverUtils::OutputTec360_binary(), Vmath::Sadd(), and tinysimd::sqrt().

◆ v_ReSize()

void Nektar::SolverUtils::StatLagrangianPoints::v_ReSize ( int Np )

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 135 of file FilterLagrangianPoints.cpp.

{
    if (Np > m_totPts)
    {
        for (int i = 0; i < m_coords.size(); ++i)
        {
            for (int d = 0; d < m_dim; ++d)
            {
                m_coords[i][d]   = Array<OneD, NekDouble>(Np);
                m_velocity[i][d] = Array<OneD, NekDouble>(Np);
            }
        }
    }
    m_totPts = Np;
    m_localIDToGlobal.clear();
    m_globalIDToLocal.clear();
}

References Nektar::UnitTests::d(), m_coords, Nektar::SolverUtils::StationaryPoints::m_dim, Nektar::SolverUtils::StationaryPoints::m_globalIDToLocal, Nektar::SolverUtils::StationaryPoints::m_localIDToGlobal, Nektar::SolverUtils::StationaryPoints::m_totPts, and m_velocity.

◆ v_SetCoords()

void Nektar::SolverUtils::StatLagrangianPoints::v_SetCoords	(	int	pid,
		const Array< OneD, NekDouble > &	gcoords
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 59 of file FilterLagrangianPoints.cpp.

{
    for (int d = 0; d < m_dim; ++d)
    {
        m_coords[0][d][id] = gcoords[d];
    }
}

References Nektar::UnitTests::d(), m_coords, and Nektar::SolverUtils::StationaryPoints::m_dim.

◆ v_SetPhysics()

void Nektar::SolverUtils::StatLagrangianPoints::v_SetPhysics	(	int	pid,
		const Array< OneD, NekDouble > &	data
	)

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 76 of file FilterLagrangianPoints.cpp.

{
    for (int d = 0; d < m_dim; ++d)
    {
        m_velocity[0][d][id] = data[d];
    }
    for (size_t d = 0; d < m_extraPhysVars.size(); ++d)
    {
        m_extraPhysics[d][id] = data[d + m_dim];
    }
}

References Nektar::UnitTests::d(), Nektar::SolverUtils::StationaryPoints::m_dim, m_extraPhysics, m_extraPhysVars, and m_velocity.

◆ v_TimeAdvance()

void Nektar::SolverUtils::StatLagrangianPoints::v_TimeAdvance ( int order )

overrideprotectedvirtual

Implements Nektar::SolverUtils::StationaryPoints.

Definition at line 104 of file FilterLagrangianPoints.cpp.

{
    order = min(order, m_intOrder);
    if (order == 1)
    {
        for (int d = 0; d < m_dim; ++d)
        {
            Vmath::Svtvp(m_totPts, m_dt, m_velocity[0][d], 1, m_coords[0][d], 1,
                         m_coords[0][d], 1);
        }
        RollOver(m_velocity);
    }
    else if (order == 2)
    {
        for (int d = 0; d < m_dim; ++d)
        {
            Array<OneD, NekDouble> res = m_coords[1][d];
            Vmath::Svtvm(m_totPts, 3., m_velocity[0][d], 1, m_velocity[1][d], 1,
                         res, 1);
            Vmath::Svtvp(m_totPts, 0.5 * m_dt, res, 1, m_coords[0][d], 1, res,
                         1);
        }
        RollOver(m_coords);
        RollOver(m_velocity);
    }
    else
    {
        ASSERTL0(false, "Integration order not supported.");
    }
}

References ASSERTL0, Nektar::UnitTests::d(), m_coords, Nektar::SolverUtils::StationaryPoints::m_dim, m_dt, m_intOrder, Nektar::SolverUtils::StationaryPoints::m_totPts, m_velocity, Nektar::SolverUtils::RollOver(), Vmath::Svtvm(), and Vmath::Svtvp().

Friends And Related Function Documentation

◆ MemoryManager< StatLagrangianPoints >

friend class MemoryManager< StatLagrangianPoints >

friend

Definition at line 1 of file FilterLagrangianPoints.h.

Member Data Documentation

◆ m_coords

Array<OneD, Array<OneD, Array<OneD, NekDouble> > > Nektar::SolverUtils::StatLagrangianPoints::m_coords

protected

Definition at line 87 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), v_AssignPoint(), v_GetCoords(), v_OutputData(), v_ReSize(), v_SetCoords(), and v_TimeAdvance().

◆ m_dt

NekDouble Nektar::SolverUtils::StatLagrangianPoints::m_dt

protected

Definition at line 91 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), and v_TimeAdvance().

◆ m_extraPhysics

Array<OneD, Array<OneD, NekDouble> > Nektar::SolverUtils::StatLagrangianPoints::m_extraPhysics

protected

Definition at line 89 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), v_OutputData(), and v_SetPhysics().

◆ m_extraPhysVars

std::vector<std::string> Nektar::SolverUtils::StatLagrangianPoints::m_extraPhysVars

protected

Definition at line 90 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), v_OutputData(), and v_SetPhysics().

◆ m_idOffset

int Nektar::SolverUtils::StatLagrangianPoints::m_idOffset

protected

Definition at line 92 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints().

◆ m_intOrder

int Nektar::SolverUtils::StatLagrangianPoints::m_intOrder

protected

Definition at line 94 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), and v_TimeAdvance().

◆ m_N

std::vector<int> Nektar::SolverUtils::StatLagrangianPoints::m_N

protected

Definition at line 93 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), and v_OutputData().

◆ m_velocity

Array<OneD, Array<OneD, Array<OneD, NekDouble> > > Nektar::SolverUtils::StatLagrangianPoints::m_velocity

protected

Definition at line 88 of file FilterLagrangianPoints.h.

Referenced by StatLagrangianPoints(), v_GetPhysics(), v_OutputData(), v_ReSize(), v_SetPhysics(), and v_TimeAdvance().

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Constructor & Destructor Documentation

◆ StatLagrangianPoints()

◆ ~StatLagrangianPoints()

Member Function Documentation

◆ create()

◆ v_AssignPoint()

◆ v_GetCoords()

◆ v_GetPhysics()

◆ v_OutputData()

◆ v_ReSize()

◆ v_SetCoords()

◆ v_SetPhysics()

◆ v_TimeAdvance()

Friends And Related Function Documentation

◆ MemoryManager< StatLagrangianPoints >

Member Data Documentation

◆ m_coords

◆ m_dt

◆ m_extraPhysics

◆ m_extraPhysVars

◆ m_idOffset

◆ m_intOrder

◆ m_N

◆ m_velocity