#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) override
	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) override
	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) override

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) override
	Advects Volume Flux. More...

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) override
	Advects Trace Flux. More...

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

Protected Member Functions inherited from Nektar::SolverUtils::Advection
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 51 of file AdvectionWeakDG.cpp.

52 {

53 }

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
	)

overrideprotectedvirtual

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 233 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`
	)

overrideprotectedvirtual

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 80 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};
     }
     LibUtilities::Timer timer;
     timer.Start();
     AdvectionWeakDG::v_AdvectCoeffs(nConvectiveFields, fields, advVel, inarray,
                                     tmp, time, pFwd, pBwd);
     timer.Stop();
     timer.AccumulateRegion("AdvWeakDG:v_AdvectCoeffs", 2);
  
     // why was this broken in many loops over convective fields?
     // this is terrible for locality
  
     timer.Start();
     for (int i = 0; i < nConvectiveFields; ++i)
     {
         fields[i]->BwdTrans(tmp[i], outarray[i]);
     }
     timer.Stop();
     timer.AccumulateRegion("AdvWeakDG:_BwdTrans", 10);
     timer1.Stop();
     timer1.AccumulateRegion("AdvWeakDG:All", 10);
 }

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`
	)

overrideprotectedvirtual

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 121 of file AdvectionWeakDG.cpp.

 {
     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", 3);
  
     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", 10);
  
     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", 10);
  
         timer.Start();
         fields[i]->MultiplyByElmtInvMass(outarray[i], outarray[i]);
         timer.Stop();
         timer.AccumulateRegion("AdvWeakDG:_MultiplyByElmtInvMass", 10);
     }
 }

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	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`
	)

overrideprotectedvirtual

Advects Trace Flux.

calculate the advection flux in the cell the trace integration

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 190 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", 10);
 }

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 > &	pVolumeFlux,
		const NekDouble &	time
	)

inlineoverrideprotectedvirtual

Advects Volume Flux.

Reimplemented from Nektar::SolverUtils::Advection.

Definition at line 86 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
	)

overrideprotectedvirtual

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 62 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, "Weak DG")

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