#include <PreconCfsBRJ.h>

Inheritance diagram for Nektar::PreconCfsBRJ:

Public Member Functions
	PreconCfsBRJ (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm)

	~PreconCfsBRJ ()

virtual bool	UpdatePreconMatCheck (const Array< OneD, const NekDouble > &res, const NekDouble dtLambda)

Public Member Functions inherited from Nektar::PreconCfsOp
	PreconCfsOp (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm)

virtual	~PreconCfsOp ()

void	SetOperators (const NekPreconCfsOperators &in)

Public Member Functions inherited from Nektar::PreconCfs
	PreconCfs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm)

virtual	~PreconCfs ()

void	DoPreconCfs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const bool &flag)

void	BuildPreconCfs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, const Array< OneD, NekDouble >> &intmp, const NekDouble time, const NekDouble lambda)

void	InitObject ()

Static Public Member Functions
static PreconCfsOpSharedPtr	create (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const LibUtilities::SessionReaderSharedPtr &pSession, const LibUtilities::CommSharedPtr &vComm)
	Creates an instance of this class. More...

Static Public Attributes
static std::string	className
	Name of the class. More...

Protected Member Functions
virtual void	v_InitObject ()

Protected Member Functions inherited from Nektar::PreconCfs
void	DoNullPrecon (const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const bool &flag)

Protected Attributes
int	m_PreconItsStep

int	m_BRJRelaxParam

Array< OneD, Array< OneD, SNekBlkMatSharedPtr > >	m_PreconMatVarsSingle

SNekBlkMatSharedPtr	m_PreconMatSingle

Array< OneD, SNekBlkMatSharedPtr >	m_TraceJacSingle

TensorOfArray4D< NekSingle >	m_TraceJacArraySingle

Array< OneD, SNekBlkMatSharedPtr >	m_TraceJacDerivSingle

TensorOfArray4D< NekSingle >	m_TraceJacDerivArraySingle

Array< OneD, Array< OneD, NekSingle > >	m_TraceJacDerivSignSingle

TensorOfArray5D< NekSingle >	m_TraceIPSymJacArraySingle

PrecType	m_PreconMatStorage

Protected Attributes inherited from Nektar::PreconCfsOp
NekPreconCfsOperators	m_operator

Protected Attributes inherited from Nektar::PreconCfs
LibUtilities::CommSharedPtr	m_Comm

bool	m_root

bool	m_verbose

int	m_spacedim

int	m_PreconMatFreezNumb

int	m_PreconTimesCounter

NekDouble	m_DtLambdaPreconMat = -1.0

NekDouble	m_BndEvaluateTime

bool	m_CalcPreconMatFlag = false

Private Member Functions
virtual void	v_DoPreconCfs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, NekDouble > &pInput, Array< OneD, NekDouble > &pOutput, const bool &flag)

virtual void	v_BuildPreconCfs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, const Array< OneD, NekDouble >> &intmp, const NekDouble time, const NekDouble lambda)

template<typename DataType , typename TypeNekBlkMatSharedPtr >
void	PreconBlkDiag (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const Array< OneD, NekDouble > &inarray, Array< OneD, NekDouble > &outarray, const TypeNekBlkMatSharedPtr &PreconMatVars, const DataType &tmpDataType)

template<typename DataType >
void	MinusOffDiag2Rhs (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, const int nvariables, const int nCoeffs, const Array< OneD, const Array< OneD, NekDouble >> &inarray, Array< OneD, Array< OneD, NekDouble >> &outarray, bool flagUpdateDervFlux, Array< OneD, Array< OneD, NekDouble >> &FwdFluxDeriv, Array< OneD, Array< OneD, NekDouble >> &BwdFluxDeriv, TensorOfArray3D< NekDouble > &qfield, TensorOfArray3D< NekDouble > &wspTrace, Array< OneD, Array< OneD, DataType >> &wspTraceDataType, const TensorOfArray4D< DataType > &TraceJacArray, const TensorOfArray4D< DataType > &TraceJacDerivArray, const Array< OneD, const Array< OneD, DataType >> &TraceJacDerivSign, const TensorOfArray5D< DataType > &TraceIPSymJacArray)

template<typename TypeNekBlkMatSharedPtr >
void	AllocatePreconBlkDiagCoeff (const Array< OneD, MultiRegions::ExpListSharedPtr > &pFields, Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr >> &gmtxarray, const int &nscale=1)

void	AllocateNekBlkMatDig (SNekBlkMatSharedPtr &mat, const Array< OneD, unsigned int > nrow, const Array< OneD, unsigned int > ncol)

Friends
class	MemoryManager< PreconCfsBRJ >

Detailed Description

Block Relaxed(weighted) Jacobi iterative (BRJ) Preconditioner for CFS

Solves a linear system using iterative methods.

Definition at line 46 of file PreconCfsBRJ.h.

Constructor & Destructor Documentation

◆ PreconCfsBRJ()

Nektar::PreconCfsBRJ::PreconCfsBRJ	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const LibUtilities::SessionReaderSharedPtr &	pSession,
		const LibUtilities::CommSharedPtr &	vComm
	)

Definition at line 53 of file PreconCfsBRJ.cpp.

     : PreconCfsOp(pFields, pSession, vComm)
 {
     pSession->LoadParameter("PreconItsStep", m_PreconItsStep, 7);
     pSession->LoadParameter("BRJRelaxParam", m_BRJRelaxParam, 1.0);
  
     int nvariables     = pFields.size();
     m_PreconMatStorage = eDiagonal;
  
     m_PreconMatVarsSingle = TensorOfArray2D<SNekBlkMatSharedPtr>(nvariables);
     for (int i = 0; i < nvariables; i++)
     {
         m_PreconMatVarsSingle[i] = Array<OneD, SNekBlkMatSharedPtr>(nvariables);
     }
     AllocatePreconBlkDiagCoeff(pFields, m_PreconMatVarsSingle);
  
     int nelmts = pFields[0]->GetNumElmts();
     int nelmtcoef;
     Array<OneD, unsigned int> nelmtmatdim(nelmts);
     for (int i = 0; i < nelmts; i++)
     {
         nelmtcoef      = pFields[0]->GetExp(i)->GetNcoeffs();
         nelmtmatdim[i] = nelmtcoef * nvariables;
     }
     AllocateNekBlkMatDig(m_PreconMatSingle, nelmtmatdim, nelmtmatdim);
 }

References AllocateNekBlkMatDig(), AllocatePreconBlkDiagCoeff(), Nektar::eDiagonal, m_BRJRelaxParam, m_PreconItsStep, m_PreconMatSingle, m_PreconMatStorage, and m_PreconMatVarsSingle.

◆ ~PreconCfsBRJ()

Nektar::PreconCfsBRJ::~PreconCfsBRJ ( )

inline

Definition at line 68 of file PreconCfsBRJ.h.

68 {};

Member Function Documentation

◆ AllocateNekBlkMatDig()

void Nektar::PreconCfsBRJ::AllocateNekBlkMatDig	(	SNekBlkMatSharedPtr &	mat,
		const Array< OneD, unsigned int >	nrow,
		const Array< OneD, unsigned int >	ncol
	)

inlineprivate

Definition at line 130 of file PreconCfsBRJ.h.

     {
         mat =
             MemoryManager<SNekBlkMat>::AllocateSharedPtr(nrow, ncol, eDIAGONAL);
         SNekMatSharedPtr loc_matNvar;
         for (int nelm = 0; nelm < nrow.size(); ++nelm)
         {
             int nrowsVars = nrow[nelm];
             int ncolsVars = ncol[nelm];
  
             loc_matNvar = MemoryManager<SNekMat>::AllocateSharedPtr(
                 nrowsVars, ncolsVars, 0.0);
             mat->SetBlock(nelm, nelm, loc_matNvar);
         }
     }

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

Referenced by AllocatePreconBlkDiagCoeff(), and PreconCfsBRJ().

◆ AllocatePreconBlkDiagCoeff()

template<typename TypeNekBlkMatSharedPtr >

void Nektar::PreconCfsBRJ::AllocatePreconBlkDiagCoeff	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		Array< OneD, Array< OneD, TypeNekBlkMatSharedPtr >> &	gmtxarray,
		const int &	nscale = `1`
	)

private

Definition at line 422 of file PreconCfsBRJ.cpp.

 {
  
     int nvars  = pFields.size();
     int nelmts = pFields[0]->GetNumElmts();
     int nelmtcoef;
     Array<OneD, unsigned int> nelmtmatdim(nelmts);
     for (int i = 0; i < nelmts; i++)
     {
         nelmtcoef      = pFields[0]->GetExp(i)->GetNcoeffs();
         nelmtmatdim[i] = nelmtcoef * nscale;
     }
  
     for (int i = 0; i < nvars; i++)
     {
         for (int j = 0; j < nvars; j++)
         {
             AllocateNekBlkMatDig(gmtxarray[i][j], nelmtmatdim, nelmtmatdim);
         }
     }
 }

References AllocateNekBlkMatDig().

Referenced by PreconCfsBRJ().

◆ create()

static PreconCfsOpSharedPtr Nektar::PreconCfsBRJ::create	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const LibUtilities::SessionReaderSharedPtr &	pSession,
		const LibUtilities::CommSharedPtr &	vComm
	)

inlinestatic

Creates an instance of this class.

Definition at line 52 of file PreconCfsBRJ.h.

     {
         PreconCfsOpSharedPtr p = MemoryManager<PreconCfsBRJ>::AllocateSharedPtr(
             pFields, pSession, vComm);
         return p;
     }

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

◆ MinusOffDiag2Rhs()

template<typename DataType >

void Nektar::PreconCfsBRJ::MinusOffDiag2Rhs	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const int	nvariables,
		const int	nCoeffs,
		const Array< OneD, const Array< OneD, NekDouble >> &	inarray,
		Array< OneD, Array< OneD, NekDouble >> &	outarray,
		bool	flagUpdateDervFlux,
		Array< OneD, Array< OneD, NekDouble >> &	FwdFluxDeriv,
		Array< OneD, Array< OneD, NekDouble >> &	BwdFluxDeriv,
		TensorOfArray3D< NekDouble > &	qfield,
		TensorOfArray3D< NekDouble > &	wspTrace,
		Array< OneD, Array< OneD, DataType >> &	wspTraceDataType,
		const TensorOfArray4D< DataType > &	TraceJacArray,
		const TensorOfArray4D< DataType > &	TraceJacDerivArray,
		const Array< OneD, const Array< OneD, DataType >> &	TraceJacDerivSign,
		const TensorOfArray5D< DataType > &	TraceIPSymJacArray
	)

private

Definition at line 300 of file PreconCfsBRJ.cpp.

 {
     boost::ignore_unused(flagUpdateDervFlux, qfield, TraceJacDerivArray,
                          TraceJacDerivSign, FwdFluxDeriv, BwdFluxDeriv,
                          TraceIPSymJacArray);
  
     int nTracePts = pFields[0]->GetTrace()->GetNpoints();
     int npoints   = pFields[0]->GetNpoints();
     int nDim      = m_spacedim;
  
     Array<OneD, Array<OneD, NekDouble>> outpnts(nvariables);
     for (int i = 0; i < nvariables; i++)
     {
         outpnts[i] = Array<OneD, NekDouble>(npoints, 0.0);
         pFields[i]->BwdTrans(outarray[i], outpnts[i]);
     }
  
     // Store forwards/backwards space along trace space
     Array<OneD, Array<OneD, NekDouble>> Fwd;
     Array<OneD, Array<OneD, NekDouble>> Bwd;
     Array<OneD, Array<OneD, NekDouble>> FwdFlux;
     Array<OneD, Array<OneD, NekDouble>> BwdFlux;
     TensorOfArray3D<NekDouble> numDerivBwd(nDim);
     TensorOfArray3D<NekDouble> numDerivFwd(nDim);
     int indexwspTrace = 0;
     Fwd               = wspTrace[indexwspTrace], indexwspTrace++;
     Bwd               = wspTrace[indexwspTrace], indexwspTrace++;
     FwdFlux           = wspTrace[indexwspTrace], indexwspTrace++;
     BwdFlux           = wspTrace[indexwspTrace], indexwspTrace++;
  
     LibUtilities::Timer timer;
     for (int i = 0; i < nvariables; ++i)
     {
         timer.Start();
         pFields[i]->GetFwdBwdTracePhys(outpnts[i], Fwd[i], Bwd[i]);
         timer.Stop();
         timer.AccumulateRegion("ExpList::GetFwdBwdTracePhys", 10);
     }
  
     int indexwspTraceDataType = 0;
     Array<OneD, Array<OneD, DataType>> Fwdarray(nvariables);
     for (int m = 0; m < nvariables; ++m)
     {
         Fwdarray[m] = wspTraceDataType[indexwspTraceDataType],
         indexwspTraceDataType++;
     }
     Array<OneD, DataType> Fwdreslt;
     Fwdreslt = wspTraceDataType[indexwspTraceDataType], indexwspTraceDataType++;
  
     for (int m = 0; m < nvariables; ++m)
     {
         for (int i = 0; i < nTracePts; ++i)
         {
             Fwdarray[m][i] = DataType(Fwd[m][i]);
         }
     }
     for (int m = 0; m < nvariables; ++m)
     {
         Vmath::Zero(nTracePts, Fwdreslt, 1);
         for (int n = 0; n < nvariables; ++n)
         {
             Vmath::Vvtvp(nTracePts, TraceJacArray[0][m][n], 1, Fwdarray[n], 1,
                          Fwdreslt, 1, Fwdreslt, 1);
         }
  
         for (int i = 0; i < nTracePts; ++i)
         {
             FwdFlux[m][i] = NekDouble(Fwdreslt[i]);
         }
     }
  
     for (int m = 0; m < nvariables; ++m)
     {
         for (int i = 0; i < nTracePts; ++i)
         {
             Fwdarray[m][i] = DataType(Bwd[m][i]);
         }
     }
     for (int m = 0; m < nvariables; ++m)
     {
         Vmath::Zero(nTracePts, Fwdreslt, 1);
         for (int n = 0; n < nvariables; ++n)
         {
             Vmath::Vvtvp(nTracePts, TraceJacArray[1][m][n], 1, Fwdarray[n], 1,
                          Fwdreslt, 1, Fwdreslt, 1);
         }
         for (int i = 0; i < nTracePts; ++i)
         {
             BwdFlux[m][i] = NekDouble(Fwdreslt[i]);
         }
     }
  
     for (int i = 0; i < nvariables; ++i)
     {
         Vmath::Fill(nCoeffs, 0.0, outarray[i], 1);
         timer.Start();
         pFields[i]->AddTraceIntegralToOffDiag(FwdFlux[i], BwdFlux[i],
                                               outarray[i]);
         timer.Stop();
         timer.AccumulateRegion("ExpList::AddTraceIntegralToOffDiag", 10);
     }
  
     for (int i = 0; i < nvariables; ++i)
     {
         Vmath::Svtvp(nCoeffs, -m_DtLambdaPreconMat, outarray[i], 1, inarray[i],
                      1, outarray[i], 1);
     }
 }

References Nektar::LibUtilities::Timer::AccumulateRegion(), Vmath::Fill(), Nektar::PreconCfs::m_DtLambdaPreconMat, Nektar::PreconCfs::m_spacedim, Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), Vmath::Svtvp(), Vmath::Vvtvp(), and Vmath::Zero().

Referenced by v_DoPreconCfs().

◆ PreconBlkDiag()

template<typename DataType , typename TypeNekBlkMatSharedPtr >

void Nektar::PreconCfsBRJ::PreconBlkDiag	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const Array< OneD, NekDouble > &	inarray,
		Array< OneD, NekDouble > &	outarray,
		const TypeNekBlkMatSharedPtr &	PreconMatVars,
		const DataType &	tmpDataType
	)

private

Definition at line 245 of file PreconCfsBRJ.cpp.

 {
     boost::ignore_unused(tmpDataType);
  
     unsigned int nvariables = pFields.size();
     unsigned int npoints    = pFields[0]->GetNcoeffs();
     unsigned int npointsVar = nvariables * npoints;
     Array<OneD, DataType> Sinarray(npointsVar);
     Array<OneD, DataType> Soutarray(npointsVar);
     NekVector<DataType> tmpVect(npointsVar, Sinarray, eWrapper);
     NekVector<DataType> outVect(npointsVar, Soutarray, eWrapper);
  
     std::shared_ptr<LocalRegions::ExpansionVector> expvect =
         pFields[0]->GetExp();
     int nTotElmt = (*expvect).size();
  
     for (int m = 0; m < nvariables; m++)
     {
         int nVarOffset = m * npoints;
         for (int ne = 0; ne < nTotElmt; ne++)
         {
             int nCoefOffset = pFields[0]->GetCoeff_Offset(ne);
             int nElmtCoef   = pFields[0]->GetNcoeffs(ne);
             int inOffset    = nVarOffset + nCoefOffset;
             int outOffset   = nCoefOffset * nvariables + m * nElmtCoef;
             for (int i = 0; i < nElmtCoef; i++)
             {
                 Sinarray[outOffset + i] = DataType(inarray[inOffset + i]);
             }
         }
     }
  
     outVect = (*PreconMatVars) * tmpVect;
  
     for (int m = 0; m < nvariables; m++)
     {
         int nVarOffset = m * npoints;
         for (int ne = 0; ne < nTotElmt; ne++)
         {
             int nCoefOffset = pFields[0]->GetCoeff_Offset(ne);
             int nElmtCoef   = pFields[0]->GetNcoeffs(ne);
             int inOffset    = nVarOffset + nCoefOffset;
             int outOffset   = nCoefOffset * nvariables + m * nElmtCoef;
             for (int i = 0; i < nElmtCoef; i++)
             {
                 outarray[inOffset + i] = NekDouble(Soutarray[outOffset + i]);
             }
         }
     }
 }

References Nektar::eWrapper.

Referenced by v_DoPreconCfs().

◆ UpdatePreconMatCheck()

bool Nektar::PreconCfsBRJ::UpdatePreconMatCheck	(	const Array< OneD, const NekDouble > &	res,
		const NekDouble	dtLambda
	)

virtual

Reimplemented from Nektar::PreconCfs.

Definition at line 223 of file PreconCfsBRJ.cpp.

 {
     boost::ignore_unused(res);
  
     bool flag = false;
  
     if (m_CalcPreconMatFlag || (m_DtLambdaPreconMat != dtLambda))
     {
         flag = true;
     }
  
     if (m_PreconMatFreezNumb < m_PreconTimesCounter)
     {
         flag = true;
     }
  
     m_CalcPreconMatFlag = flag;
     return flag;
 }

References Nektar::PreconCfs::m_CalcPreconMatFlag, Nektar::PreconCfs::m_DtLambdaPreconMat, Nektar::PreconCfs::m_PreconMatFreezNumb, and Nektar::PreconCfs::m_PreconTimesCounter.

◆ v_BuildPreconCfs()

void Nektar::PreconCfsBRJ::v_BuildPreconCfs	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const Array< OneD, const Array< OneD, NekDouble >> &	intmp,
		const NekDouble	time,
		const NekDouble	lambda
	)

privatevirtual

Reimplemented from Nektar::PreconCfsOp.

Definition at line 197 of file PreconCfsBRJ.cpp.

 {
     if (0 < m_PreconItsStep)
     {
         m_operator.DoCalcPreconMatBRJCoeff(
             intmp, m_PreconMatVarsSingle, m_PreconMatSingle, m_TraceJacSingle,
             m_TraceJacDerivSingle, m_TraceJacDerivSignSingle,
             m_TraceJacArraySingle, m_TraceJacDerivArraySingle,
             m_TraceIPSymJacArraySingle);
  
         if (m_verbose && m_root)
         {
             cout << "     ## CalcuPreconMat " << endl;
         }
     }
  
     m_BndEvaluateTime   = time;
     m_DtLambdaPreconMat = lambda;
  
     m_CalcPreconMatFlag  = false;
     m_PreconTimesCounter = 1;
 }

References Nektar::NekPreconCfsOperators::DoCalcPreconMatBRJCoeff(), Nektar::PreconCfs::m_BndEvaluateTime, Nektar::PreconCfs::m_CalcPreconMatFlag, Nektar::PreconCfs::m_DtLambdaPreconMat, Nektar::PreconCfsOp::m_operator, m_PreconItsStep, m_PreconMatSingle, m_PreconMatVarsSingle, Nektar::PreconCfs::m_PreconTimesCounter, Nektar::PreconCfs::m_root, m_TraceIPSymJacArraySingle, m_TraceJacArraySingle, m_TraceJacDerivArraySingle, m_TraceJacDerivSignSingle, m_TraceJacDerivSingle, m_TraceJacSingle, and Nektar::PreconCfs::m_verbose.

◆ v_DoPreconCfs()

void Nektar::PreconCfsBRJ::v_DoPreconCfs	(	const Array< OneD, MultiRegions::ExpListSharedPtr > &	pFields,
		const Array< OneD, NekDouble > &	pInput,
		Array< OneD, NekDouble > &	pOutput,
		const bool &	flag
	)

privatevirtual

Reimplemented from Nektar::PreconCfsOp.

Definition at line 88 of file PreconCfsBRJ.cpp.

 {
     boost::ignore_unused(flag);
  
     int nBRJIterTot = m_PreconItsStep;
     if (0 == nBRJIterTot)
     {
         DoNullPrecon(inarray, outarray, flag);
     }
     else
     {
         const NekDouble BRJParam   = m_BRJRelaxParam;
         const NekDouble OmBRJParam = 1.0 - BRJParam;
  
         unsigned int nvariables = pFields.size();
         unsigned int npoints    = pFields[0]->GetNcoeffs();
         unsigned int ntotpnt    = inarray.size();
  
         ASSERTL0(nvariables * npoints == ntotpnt,
                  "nvariables*npoints!=ntotpnt in PreconCoeff");
  
         Array<OneD, NekDouble> rhs(ntotpnt);
  
         Array<OneD, NekDouble> outN(ntotpnt);
         Array<OneD, NekDouble> outTmp(ntotpnt);
         Array<OneD, Array<OneD, NekDouble>> rhs2d(nvariables);
         Array<OneD, Array<OneD, NekDouble>> out_2d(nvariables);
         Array<OneD, Array<OneD, NekDouble>> outTmp_2d(nvariables);
         for (int m = 0; m < nvariables; m++)
         {
             int moffset  = m * npoints;
             rhs2d[m]     = rhs + moffset;
             out_2d[m]    = outarray + moffset;
             outTmp_2d[m] = outTmp + moffset;
             pFields[m]->MultiplyByMassMatrix(inarray + moffset, rhs2d[m]);
         }
  
         int nphysic   = pFields[0]->GetNpoints();
         int nTracePts = pFields[0]->GetTrace()->GetNpoints();
         TensorOfArray3D<NekDouble> qfield(m_spacedim);
         for (int i = 0; i < m_spacedim; i++)
         {
             qfield[i] = Array<OneD, Array<OneD, NekDouble>>(nvariables);
             for (int j = 0; j < nvariables; j++)
             {
                 qfield[i][j] = Array<OneD, NekDouble>(nphysic);
             }
         }
         int ntmpTrace = 4 + 2 * m_spacedim;
         TensorOfArray3D<NekDouble> tmpTrace(ntmpTrace);
         for (int i = 0; i < ntmpTrace; i++)
         {
             tmpTrace[i] = Array<OneD, Array<OneD, NekDouble>>(nvariables);
             for (int j = 0; j < nvariables; j++)
             {
                 tmpTrace[i][j] = Array<OneD, NekDouble>(nTracePts);
             }
         }
         Array<OneD, Array<OneD, NekDouble>> FwdFluxDeriv(nvariables);
         Array<OneD, Array<OneD, NekDouble>> BwdFluxDeriv(nvariables);
         for (int j = 0; j < nvariables; j++)
         {
             FwdFluxDeriv[j] = Array<OneD, NekDouble>(nTracePts);
             BwdFluxDeriv[j] = Array<OneD, NekDouble>(nTracePts);
         }
  
         bool flagUpdateDervFlux = false;
  
         const int nwspTraceDataType = nvariables + 1;
         NekSingle tmpSingle;
         Array<OneD, Array<OneD, NekSingle>> wspTraceDataType(nwspTraceDataType);
         for (int m = 0; m < nwspTraceDataType; m++)
         {
             wspTraceDataType[m] = Array<OneD, NekSingle>(nTracePts);
         }
  
         LibUtilities::Timer timer;
         timer.Start();
         PreconBlkDiag(pFields, rhs, outarray, m_PreconMatSingle, tmpSingle);
         timer.Stop();
         timer.AccumulateRegion("PreconCfsBRJ::PreconBlkDiag", 2);
  
         for (int nrelax = 0; nrelax < nBRJIterTot - 1; nrelax++)
         {
             Vmath::Smul(ntotpnt, OmBRJParam, outarray, 1, outN, 1);
  
             timer.Start();
             MinusOffDiag2Rhs(
                 pFields, nvariables, npoints, rhs2d, out_2d, flagUpdateDervFlux,
                 FwdFluxDeriv, BwdFluxDeriv, qfield, tmpTrace, wspTraceDataType,
                 m_TraceJacArraySingle, m_TraceJacDerivArraySingle,
                 m_TraceJacDerivSignSingle, m_TraceIPSymJacArraySingle);
             timer.Stop();
             timer.AccumulateRegion("PreconCfsBRJ::MinusOffDiag2Rhs", 2);
  
             timer.Start();
             PreconBlkDiag(pFields, outarray, outTmp, m_PreconMatSingle,
                           tmpSingle);
             timer.Stop();
             timer.AccumulateRegion("PreconCfsBRJ::PreconBlkDiag", 2);
  
             Vmath::Svtvp(ntotpnt, BRJParam, outTmp, 1, outN, 1, outarray, 1);
         }
     }
 }

References Nektar::LibUtilities::Timer::AccumulateRegion(), ASSERTL0, Nektar::PreconCfs::DoNullPrecon(), m_BRJRelaxParam, m_PreconItsStep, m_PreconMatSingle, Nektar::PreconCfs::m_spacedim, m_TraceIPSymJacArraySingle, m_TraceJacArraySingle, m_TraceJacDerivArraySingle, m_TraceJacDerivSignSingle, MinusOffDiag2Rhs(), PreconBlkDiag(), Vmath::Smul(), Nektar::LibUtilities::Timer::Start(), Nektar::LibUtilities::Timer::Stop(), and Vmath::Svtvp().

◆ v_InitObject()

void Nektar::PreconCfsBRJ::v_InitObject ( )

protectedvirtual

Reimplemented from Nektar::PreconCfsOp.

Definition at line 83 of file PreconCfsBRJ.cpp.

 {
     PreconCfsOp::v_InitObject();
 }

References Nektar::PreconCfsOp::v_InitObject().

Friends And Related Function Documentation

◆ MemoryManager< PreconCfsBRJ >

friend class MemoryManager< PreconCfsBRJ >

friend

Definition at line 1 of file PreconCfsBRJ.h.

Member Data Documentation

◆ className

std::string Nektar::PreconCfsBRJ::className

static

Initial value:

=
    GetPreconCfsOpFactory().RegisterCreatorFunction(
        "PreconCfsBRJ", PreconCfsBRJ::create,
        "Block Relaxed Jacobi Preconditioner for CFS.")

Name of the class.

Definition at line 63 of file PreconCfsBRJ.h.

◆ m_BRJRelaxParam

int Nektar::PreconCfsBRJ::m_BRJRelaxParam

protected

Definition at line 75 of file PreconCfsBRJ.h.

Referenced by PreconCfsBRJ(), and v_DoPreconCfs().

◆ m_PreconItsStep

int Nektar::PreconCfsBRJ::m_PreconItsStep

protected

Definition at line 74 of file PreconCfsBRJ.h.

Referenced by PreconCfsBRJ(), v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_PreconMatSingle

SNekBlkMatSharedPtr Nektar::PreconCfsBRJ::m_PreconMatSingle

protected

Definition at line 78 of file PreconCfsBRJ.h.

Referenced by PreconCfsBRJ(), v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_PreconMatStorage

PrecType Nektar::PreconCfsBRJ::m_PreconMatStorage

protected

Definition at line 86 of file PreconCfsBRJ.h.

Referenced by PreconCfsBRJ().

◆ m_PreconMatVarsSingle

Array<OneD, Array<OneD, SNekBlkMatSharedPtr> > Nektar::PreconCfsBRJ::m_PreconMatVarsSingle

protected

Definition at line 77 of file PreconCfsBRJ.h.

Referenced by PreconCfsBRJ(), and v_BuildPreconCfs().

◆ m_TraceIPSymJacArraySingle

TensorOfArray5D<NekSingle> Nektar::PreconCfsBRJ::m_TraceIPSymJacArraySingle

protected

Definition at line 84 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_TraceJacArraySingle

TensorOfArray4D<NekSingle> Nektar::PreconCfsBRJ::m_TraceJacArraySingle

protected

Definition at line 80 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_TraceJacDerivArraySingle

TensorOfArray4D<NekSingle> Nektar::PreconCfsBRJ::m_TraceJacDerivArraySingle

protected

Definition at line 82 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_TraceJacDerivSignSingle

Array<OneD, Array<OneD, NekSingle> > Nektar::PreconCfsBRJ::m_TraceJacDerivSignSingle

protected

Definition at line 83 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs(), and v_DoPreconCfs().

◆ m_TraceJacDerivSingle

Array<OneD, SNekBlkMatSharedPtr> Nektar::PreconCfsBRJ::m_TraceJacDerivSingle

protected

Definition at line 81 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs().

◆ m_TraceJacSingle

Array<OneD, SNekBlkMatSharedPtr> Nektar::PreconCfsBRJ::m_TraceJacSingle

protected

Definition at line 79 of file PreconCfsBRJ.h.

Referenced by v_BuildPreconCfs().

Public Member Functions

Static Public Member Functions

Static Public Attributes

Protected Member Functions

Protected Attributes

Private Member Functions

Friends

Detailed Description

Constructor & Destructor Documentation

◆ PreconCfsBRJ()

◆ ~PreconCfsBRJ()

Member Function Documentation

◆ AllocateNekBlkMatDig()

◆ AllocatePreconBlkDiagCoeff()

◆ create()

◆ MinusOffDiag2Rhs()

◆ PreconBlkDiag()

◆ UpdatePreconMatCheck()

◆ v_BuildPreconCfs()

◆ v_DoPreconCfs()

◆ v_InitObject()

Friends And Related Function Documentation

◆ MemoryManager< PreconCfsBRJ >

Member Data Documentation

◆ className

◆ m_BRJRelaxParam

◆ m_PreconItsStep

◆ m_PreconMatSingle

◆ m_PreconMatStorage

◆ m_PreconMatVarsSingle

◆ m_TraceIPSymJacArraySingle

◆ m_TraceJacArraySingle

◆ m_TraceJacDerivArraySingle

◆ m_TraceJacDerivSignSingle

◆ m_TraceJacDerivSingle

◆ m_TraceJacSingle