#include <AdvectionWeakDG.h>

Inheritance diagram for Nektar::SolverUtils::AdvectionWeakDG:

Static Public Member Functions
static AdvectionSharedPtr	create (std::string advType)

Static Public Attributes
static std::string	type

Protected Member Functions
	AdvectionWeakDG ()

virtual void	v_InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
	Initialise AdvectionWeakDG objects and store them before starting the time-stepping. More...

virtual void	v_Advect (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Compute the advection term at each time-step using the Discontinuous Galerkin approach (DG). More...

virtual void	v_AdvectCoeffs (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)

virtual void	v_AdvectVolumeFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &advVel, const Array< OneD, Array< OneD, NekDouble >> &inarray, TensorOfArray3D< NekDouble > &VolumeFlux, const NekDouble &time)

virtual void	v_AdvectTraceFlux (const int nConvective, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble >> &advVel, const Array< OneD, Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &TraceFlux, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)

virtual void	v_AddVolumJacToMat (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int &nConvectiveFields, const TensorOfArray5D< NekDouble > &ElmtJacArray, Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &gmtxarray)

Protected Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT void	v_AdvectVolumeFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, TensorOfArray3D< NekDouble > &pVolumeFlux, const NekDouble &time)
	Advects Volume Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_AdvectTraceFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble >> &pTraceFlux, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Advects Trace Flux. More...

virtual SOLVER_UTILS_EXPORT void	v_SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
	Overrides the base flow used during linearised advection. More...

Additional Inherited Members
Public Member Functions inherited from Nektar::SolverUtils::Advection
virtual SOLVER_UTILS_EXPORT	~Advection ()

SOLVER_UTILS_EXPORT void	InitObject (LibUtilities::SessionReaderSharedPtr pSession, Array< OneD, MultiRegions::ExpListSharedPtr > pFields)
	Interface function to initialise the advection object. More...

SOLVER_UTILS_EXPORT void	Advect (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Interface function to advect the vector field. More...

SOLVER_UTILS_EXPORT void	AdvectVolumeFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, TensorOfArray3D< NekDouble > &pVolumeFlux, const NekDouble &pTime)
	Interface function to advect the Volume field. More...

SOLVER_UTILS_EXPORT void	AdvectTraceFlux (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, Array< OneD, NekDouble >> &pAdvVel, const Array< OneD, Array< OneD, NekDouble >> &pInarray, Array< OneD, Array< OneD, NekDouble >> &pTraceFlux, const NekDouble &pTime, const Array< OneD, Array< OneD, NekDouble >> &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble >> &pBwd=NullNekDoubleArrayOfArray)
	Interface function to advect the Trace field. More...

SOLVER_UTILS_EXPORT void	AdvectCoeffs (const int nConvectiveFields, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields, const Array< OneD, Array< OneD, NekDouble > > &advVel, const Array< OneD, Array< OneD, NekDouble > > &inarray, Array< OneD, Array< OneD, NekDouble > > &outarray, const NekDouble &time, const Array< OneD, Array< OneD, NekDouble > > &pFwd=NullNekDoubleArrayOfArray, const Array< OneD, Array< OneD, NekDouble > > &pBwd=NullNekDoubleArrayOfArray)
	Similar with Advection::Advect(): calculate the advection flux The difference is in the outarray: it is the coefficients of basis for AdvectCoeffs() it is the physics on quadrature points for Advect() More...

template<typename FuncPointerT , typename ObjectPointerT >
void	SetFluxVector (FuncPointerT func, ObjectPointerT obj)
	Set the flux vector callback function. More...

void	SetRiemannSolver (RiemannSolverSharedPtr riemann)
	Set a Riemann solver object for this advection object. More...

void	SetFluxVector (AdvectionFluxVecCB fluxVector)
	Set the flux vector callback function. More...

void	SetBaseFlow (const Array< OneD, Array< OneD, NekDouble > > &inarray, const Array< OneD, MultiRegions::ExpListSharedPtr > &fields)
	Set the base flow used for linearised advection objects. More...

template<typename DataType , typename TypeNekBlkMatSharedPtr >
SOLVER_UTILS_EXPORT void	AddTraceJacToMat (const int nConvectiveFields, const int nSpaceDim, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacCons, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr > > &gmtxarray, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacGrad, const Array< OneD, Array< OneD, DataType > > &TracePntJacGradSign)

template<typename DataType , typename TypeNekBlkMatSharedPtr >
void	CalcJacobTraceInteg (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int m, const int n, const Array< OneD, const TypeNekBlkMatSharedPtr > &PntJac, const Array< OneD, const Array< OneD, DataType > > &PntJacSign, Array< OneD, DNekMatSharedPtr > &TraceJacFwd, Array< OneD, DNekMatSharedPtr > &TraceJacBwd)

SOLVER_UTILS_EXPORT void	AddVolumJacToMat (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int &nConvectiveFields, const TensorOfArray5D< NekDouble > &ElmtJacArray, Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &gmtxarray)

template<typename DataType , typename TypeNekBlkMatSharedPtr >
void	AddTraceJacToMat (const int nConvectiveFields, const int nSpaceDim, const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacCons, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr > > &gmtxarray, const Array< OneD, TypeNekBlkMatSharedPtr > &TracePntJacGrad, const Array< OneD, Array< OneD, DataType > > &TracePntJacGradSign)

Protected Attributes inherited from Nektar::SolverUtils::Advection
AdvectionFluxVecCB	m_fluxVector
	Callback function to the flux vector (set when advection is in conservative form). More...

RiemannSolverSharedPtr	m_riemann
	Riemann solver for DG-type schemes. More...

int	m_spaceDim
	Storage for space dimension. Used for homogeneous extension. More...

Detailed Description

Definition at line 46 of file AdvectionWeakDG.h.

Constructor & Destructor Documentation

◆ AdvectionWeakDG()

Nektar::SolverUtils::AdvectionWeakDG::AdvectionWeakDG ( )

protected

Definition at line 50 of file AdvectionWeakDG.cpp.

51 {

52 }

Referenced by create().

Member Function Documentation

◆ create()

static AdvectionSharedPtr Nektar::SolverUtils::AdvectionWeakDG::create ( std::string advType )

inlinestatic

Definition at line 49 of file AdvectionWeakDG.h.

             {
                 boost::ignore_unused(advType);
                 return AdvectionSharedPtr(new AdvectionWeakDG());
             }

References AdvectionWeakDG().

◆ v_AddVolumJacToMat()

void Nektar::SolverUtils::AdvectionWeakDG::v_AddVolumJacToMat	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const int &	nConvectiveFields,
		const TensorOfArray5D< NekDouble > &	ElmtJacArray,
		Array< OneD, Array< OneD, SNekBlkMatSharedPtr >> &	gmtxarray
	)

protectedvirtual

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 244 of file AdvectionWeakDG.cpp.

         {
             MultiRegions::ExpListSharedPtr explist = pFields[0];
                 std::shared_ptr<LocalRegions::ExpansionVector> pexp = explist->GetExp();
             int nTotElmt            = (*pexp).size();
             int nElmtPnt,nElmtCoef;
  
             SNekMatSharedPtr        tmpGmtx;
             DNekMatSharedPtr        ElmtMat;
  
             Array<OneD,NekSingle> GMat_data;
             Array<OneD,NekDouble> Elmt_data;
             Array<OneD,NekSingle> Elmt_dataSingle;
  
             Array<OneD, DNekMatSharedPtr>  mtxPerVar(nTotElmt);
             Array<OneD, int > elmtpnts(nTotElmt);
             Array<OneD, int > elmtcoef(nTotElmt);
             for(int nelmt = 0; nelmt < nTotElmt; nelmt++)
             {
                 nElmtCoef           = (*pexp)[nelmt]->GetNcoeffs();
                 nElmtPnt            = (*pexp)[nelmt]->GetTotPoints();
                 elmtpnts[nelmt]     =   nElmtPnt;
                 elmtcoef[nelmt]     =   nElmtCoef;
                 mtxPerVar[nelmt]    =MemoryManager<DNekMat>
                                     ::AllocateSharedPtr(nElmtCoef, nElmtPnt);
             }
  
             Array<OneD, DNekMatSharedPtr>  mtxPerVarCoeff(nTotElmt);
             for(int nelmt = 0; nelmt < nTotElmt; nelmt++)
             {
                 nElmtCoef   =   elmtcoef[nelmt];
                 mtxPerVarCoeff[nelmt]    =MemoryManager<DNekMat>
                                     ::AllocateSharedPtr(nElmtCoef, nElmtCoef);
             }
  
             for(int m = 0; m < nConvectiveFields; m++)
             {
                 for(int n = 0; n < nConvectiveFields; n++)
                 {
                     for(int nelmt = 0; nelmt < nTotElmt; nelmt++)
                     {
                         (*mtxPerVarCoeff[nelmt])   =    0.0;
                         (*mtxPerVar[nelmt])        =    0.0;
                     }
  
                     explist->GetMatIpwrtDeriveBase(ElmtJacArray[m][n], 
                         mtxPerVar);
                     // Memory can be saved by reusing mtxPerVar
                     explist->AddRightIPTBaseMatrix(mtxPerVar,mtxPerVarCoeff);
  
                     for(int nelmt = 0; nelmt < nTotElmt; nelmt++)
                     {
                         nElmtCoef = elmtcoef[nelmt];
                         nElmtPnt = elmtpnts[nelmt];
                         int ntotDofs = nElmtCoef*nElmtCoef;
  
                         if(Elmt_dataSingle.size()<ntotDofs)
                         {
                             Elmt_dataSingle = Array<OneD, NekSingle> (ntotDofs);
                         }
  
                         tmpGmtx = gmtxarray[m][n]->GetBlock(nelmt,nelmt);
                         ElmtMat = mtxPerVarCoeff[nelmt];
  
                         GMat_data = tmpGmtx->GetPtr();
                         Elmt_data = ElmtMat->GetPtr();
  
                         for(int i=0;i<ntotDofs;i++)
                         {
                             Elmt_dataSingle[i]  =   NekSingle( Elmt_data[i] );
                         }
  
                         Vmath::Vadd(ntotDofs,GMat_data,1,Elmt_dataSingle,1,GMat_data,1);
                     }
                 }
             }
         }

References Nektar::MemoryManager< DataType >::AllocateSharedPtr(), and Vmath::Vadd().

◆ v_Advect()

void Nektar::SolverUtils::AdvectionWeakDG::v_Advect	(	const int	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	advVel,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const NekDouble &	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

protectedvirtual

Compute the advection term at each time-step using the Discontinuous Galerkin approach (DG).

Parameters

nConvectiveFields	Number of fields.
fields	Pointer to fields.
advVel	Advection velocities.
inarray	Solution at the previous time-step.
outarray	Advection term to be passed at the time integration class.

Implements Nektar::SolverUtils::Advection.

Definition at line 79 of file AdvectionWeakDG.cpp.

         {
             LibUtilities::Timer timer1;
             timer1.Start();
             
             boost::ignore_unused(advVel, time);
             size_t nCoeffs         = fields[0]->GetNcoeffs();
  
             Array<OneD, Array<OneD, NekDouble> > tmp{size_t(nConvectiveFields)};
             for (int i = 0; i < nConvectiveFields; ++i)
             {
                 tmp[i] = Array<OneD, NekDouble> {nCoeffs, 0.0};
             }
  
             AdvectionWeakDG::v_AdvectCoeffs(
                 nConvectiveFields, fields, advVel, inarray, tmp, time,
                 pFwd, pBwd);
  
             // why was this broken in many loops over convective fields?
             // this is terrible for locality
             LibUtilities::Timer timer;
             timer.Start();
             for (int i = 0; i < nConvectiveFields; ++i)
             {
                 fields[i]->BwdTrans(tmp[i], outarray[i]);
             }
             timer.Stop();
             timer.AccumulateRegion("AdvWeakDG:_BwdTrans",1);
             timer1.Stop();
             timer1.AccumulateRegion("AdvWeakDG:All");
         }

References Nektar::LibUtilities::Timer::AccumulateRegion(), Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), and v_AdvectCoeffs().

◆ v_AdvectCoeffs()

void Nektar::SolverUtils::AdvectionWeakDG::v_AdvectCoeffs	(	const int	nConvective,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble > > &	advVel,
		const Array< OneD, Array< OneD, NekDouble > > &	inarray,
		Array< OneD, Array< OneD, NekDouble > > &	outarray,
		const NekDouble &	time,
		const Array< OneD, Array< OneD, NekDouble > > &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble > > &	pBwd = `NullNekDoubleArrayOfArray`
	)

protectedvirtual

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 119 of file AdvectionWeakDG.cpp.

         {
             LibUtilities::Timer timer1;
             timer1.Start();
             size_t nPointsTot      = fields[0]->GetTotPoints();
             size_t nCoeffs         = fields[0]->GetNcoeffs();
             size_t nTracePointsTot = fields[0]->GetTrace()->GetTotPoints();
  
             Array<OneD, Array<OneD, Array<OneD, NekDouble> > >
                 fluxvector(nConvectiveFields);
             // Allocate storage for flux vector F(u).
             for (int i = 0; i < nConvectiveFields; ++i)
             {
                 fluxvector[i] =
                     Array<OneD, Array<OneD, NekDouble> > {size_t(m_spaceDim)};
                 for (int j = 0; j < m_spaceDim; ++j)
                 {
                     fluxvector[i][j] = Array<OneD, NekDouble>(nPointsTot,0.0);
                 }
             }
  
             LibUtilities::Timer timer;
             timer.Start();
             v_AdvectVolumeFlux(nConvectiveFields, fields, advVel, inarray,
                                 fluxvector, time);
             timer.Stop();
             timer.AccumulateRegion("AdvWeakDG:_fluxVector",1);
  
  
             timer.Start();
             // Get the advection part (without numerical flux)
             for (int i = 0; i < nConvectiveFields; ++i)
             {
                 Vmath::Fill(outarray[i].size(), 0.0, outarray[i], 1);
                 fields[i]->IProductWRTDerivBase(fluxvector[i], outarray[i]);
             }
             timer.Stop();
             timer.AccumulateRegion("AdvWeakDG:_IProductWRTDerivBase",1);
  
             Array<OneD, Array<OneD, NekDouble> >
                 numflux{size_t(nConvectiveFields)};
  
             for (int i = 0; i < nConvectiveFields; ++i)
             {
                 numflux[i] = Array<OneD, NekDouble> {nTracePointsTot, 0.0};
             }
  
             v_AdvectTraceFlux(nConvectiveFields, fields, advVel, inarray,
                 numflux, time, pFwd, pBwd);
  
             // Evaulate <\phi, \hat{F}\cdot n> - OutField[i]
             for(int i = 0; i < nConvectiveFields; ++i)
             {
                 Vmath::Neg                      (nCoeffs, outarray[i], 1);
  
                 timer.Start();
                 fields[i]->AddTraceIntegral     (numflux[i], outarray[i]);
                 timer.Stop();
                 timer.AccumulateRegion("AdvWeakDG:_AddTraceIntegral",1);
  
                 timer.Start();
                 fields[i]->MultiplyByElmtInvMass(outarray[i], outarray[i]);
                 timer.Stop();
                 timer.AccumulateRegion("AdvWeakDG:_MultiplyByElmtInvMass",1);
             }
  
             timer1.Stop();
             timer1.AccumulateRegion("AdvWeakDG: Coeff All");
         }

References Nektar::LibUtilities::Timer::AccumulateRegion(), Vmath::Fill(), Nektar::SolverUtils::Advection::m_spaceDim, Vmath::Neg(), Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), v_AdvectTraceFlux(), and v_AdvectVolumeFlux().

Referenced by v_Advect().

◆ v_AdvectTraceFlux()

void Nektar::SolverUtils::AdvectionWeakDG::v_AdvectTraceFlux	(	const int	nConvective,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	advVel,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		Array< OneD, Array< OneD, NekDouble >> &	TraceFlux,
		const NekDouble &	time,
		const Array< OneD, Array< OneD, NekDouble >> &	pFwd = `NullNekDoubleArrayOfArray`,
		const Array< OneD, Array< OneD, NekDouble >> &	pBwd = `NullNekDoubleArrayOfArray`
	)

protectedvirtual

Definition at line 197 of file AdvectionWeakDG.cpp.

         {
             boost::ignore_unused(advVel, time);
             int nTracePointsTot = fields[0]->GetTrace()->GetTotPoints();
  
             ASSERTL1(m_riemann,
                 "Riemann solver must be provided for AdvectionWeakDG.");
  
             // Store forwards/backwards space along trace space
             Array<OneD, Array<OneD, NekDouble>> Fwd(nConvectiveFields);
             Array<OneD, Array<OneD, NekDouble>> Bwd(nConvectiveFields);
  
             if (pFwd == NullNekDoubleArrayOfArray ||
                 pBwd == NullNekDoubleArrayOfArray)
             {
                 for (int i = 0; i < nConvectiveFields; ++i)
                 {
                     Fwd[i] = Array<OneD, NekDouble>(nTracePointsTot, 0.0);
                     Bwd[i] = Array<OneD, NekDouble>(nTracePointsTot, 0.0);
                     fields[i]->GetFwdBwdTracePhys(inarray[i], Fwd[i], Bwd[i]);
                 }
             }
             else
             {
                 for (int i = 0; i < nConvectiveFields; ++i)
                 {
                     Fwd[i] = pFwd[i];
                     Bwd[i] = pBwd[i];
                 }
             }
  
             LibUtilities::Timer timer;
             timer.Start();
             m_riemann->Solve(m_spaceDim, Fwd, Bwd, TraceFlux);
             timer.Stop();
             timer.AccumulateRegion("AdvWeakDG:_Riemann",1);
  
         }

References Nektar::LibUtilities::Timer::AccumulateRegion(), ASSERTL1, Nektar::SolverUtils::Advection::m_riemann, Nektar::SolverUtils::Advection::m_spaceDim, Nektar::NullNekDoubleArrayOfArray, Nektar::LibUtilities::Timer::Start(), and Nektar::LibUtilities::Timer::Stop().

Referenced by v_AdvectCoeffs().

◆ v_AdvectVolumeFlux()

virtual void Nektar::SolverUtils::AdvectionWeakDG::v_AdvectVolumeFlux	(	const int	nConvectiveFields,
		const Array< OneD, MultiRegions::ExpListSharedPtr > &	fields,
		const Array< OneD, Array< OneD, NekDouble >> &	advVel,
		const Array< OneD, Array< OneD, NekDouble >> &	inarray,
		TensorOfArray3D< NekDouble > &	VolumeFlux,
		const NekDouble &	time
	)

inlineprotectedvirtual

Definition at line 88 of file AdvectionWeakDG.h.

             {
                 boost::ignore_unused(nConvectiveFields, fields, advVel, inarray,
                                         VolumeFlux, time);
                 m_fluxVector(inarray, VolumeFlux);
             }

References Nektar::SolverUtils::Advection::m_fluxVector.

Referenced by v_AdvectCoeffs().

◆ v_InitObject()

void Nektar::SolverUtils::AdvectionWeakDG::v_InitObject	(	LibUtilities::SessionReaderSharedPtr	pSession,
		Array< OneD, MultiRegions::ExpListSharedPtr >	pFields
	)

protectedvirtual

Initialise AdvectionWeakDG objects and store them before starting the time-stepping.

Parameters

pSession	Pointer to session reader.
pFields	Pointer to fields.

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 61 of file AdvectionWeakDG.cpp.

         {
             Advection::v_InitObject(pSession, pFields);
         }

References Nektar::SolverUtils::Advection::v_InitObject().

Member Data Documentation

◆ type

std::string Nektar::SolverUtils::AdvectionWeakDG::type

static

Initial value:

= GetAdvectionFactory().

RegisterCreatorFunction("WeakDG", AdvectionWeakDG::create)

Nektar::SolverUtils::AdvectionWeakDG::create

static AdvectionSharedPtr create(std::string advType)

Definition: AdvectionWeakDG.h:49

Nektar::SolverUtils::GetAdvectionFactory

AdvectionFactory & GetAdvectionFactory()

Gets the factory for initialising advection objects.

Definition: Advection.cpp:47

Definition at line 55 of file AdvectionWeakDG.h.

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ AdvectionWeakDG()

Member Function Documentation

◆ create()

◆ v_AddVolumJacToMat()

◆ v_Advect()

◆ v_AdvectCoeffs()

◆ v_AdvectTraceFlux()

◆ v_AdvectVolumeFlux()

◆ v_InitObject()

Member Data Documentation

◆ type