dftefe::basis::OrthoEFEOverlapInverseOpContextGLLInternal Namespace Reference

Classes
class	OverlapMatrixInverseLinearSolverFunctionFE

Functions
template<typename ValueTypeOperator , typename ValueTypeOperand , utils::MemorySpace memorySpace, size_type dim>
void	computeBasisOverlapMatrix (const FEBasisDataStorage< ValueTypeOperator, memorySpace > &classicalBlockGLLBasisDataStorage, const FEBasisDataStorage< ValueTypeOperator, memorySpace > &enrichmentBlockEnrichmentBasisDataStorage, const FEBasisDataStorage< ValueTypeOperator, memorySpace > &enrichmentBlockClassicalBasisDataStorage, std::shared_ptr< const FEBasisDofHandler< ValueTypeOperand, memorySpace, dim > > cfeBDH, std::shared_ptr< const EFEBasisDofHandler< ValueTypeOperand, ValueTypeOperator, memorySpace, dim > > efeBDH, utils::MemoryStorage< ValueTypeOperator, memorySpace > &NiNjInAllCells, linearAlgebra::LinAlgOpContext< memorySpace > &linAlgOpContext)

Function Documentation

computeBasisOverlapMatrix()

template<typename ValueTypeOperator , typename ValueTypeOperand , utils::MemorySpace memorySpace, size_type dim>

void dftefe::basis::OrthoEFEOverlapInverseOpContextGLLInternal::computeBasisOverlapMatrix	(	const FEBasisDataStorage< ValueTypeOperator, memorySpace > &	classicalBlockGLLBasisDataStorage,
		const FEBasisDataStorage< ValueTypeOperator, memorySpace > &	enrichmentBlockEnrichmentBasisDataStorage,
		const FEBasisDataStorage< ValueTypeOperator, memorySpace > &	enrichmentBlockClassicalBasisDataStorage,
		std::shared_ptr< const FEBasisDofHandler< ValueTypeOperand, memorySpace, dim > >	cfeBDH,
		std::shared_ptr< const EFEBasisDofHandler< ValueTypeOperand, ValueTypeOperator, memorySpace, dim > >	efeBDH,
		utils::MemoryStorage< ValueTypeOperator, memorySpace > &	NiNjInAllCells,
		linearAlgebra::LinAlgOpContext< memorySpace > &	linAlgOpContext
	)

Ni_pristine*Ni_pristine at quadpoints for (unsigned int qPoint = 0; qPoint < nQuadPointInCellEnrichmentBlockEnrichment; qPoint++) { basisOverlapTmpIter += (cumulativeEnrichmentBlockEnrichmentDofQuadPoints

• nQuadPointInCellEnrichmentBlockEnrichment * iNode + qPoint) * (cumulativeEnrichmentBlockEnrichmentDofQuadPoints
• nQuadPointInCellEnrichmentBlockEnrichment * jNode + qPoint) * cellJxWValuesEnrichmentBlockEnrichment[qPoint]; }

                    ValueTypeOperator NpiNpj     = (ValueTypeOperator)0,
                                      ciNciNpj   = (ValueTypeOperator)0,
                                      NpicjNcj   = (ValueTypeOperator)0,
                                      ciNcicjNcj = (ValueTypeOperator)0;

Ni_pristine*Ni_pristine at quadpoints for (unsigned int qPoint = 0; qPoint < nQuadPointInCellEnrichmentBlockEnrichment; qPoint++) { NpiNpj += (enrichmentValuesVec.data() + (iNode - dofsPerCellCFE) * nQuadPointInCellEnrichmentBlockEnrichment + qPoint) * (enrichmentValuesVec.data() + (jNode - dofsPerCellCFE) * nQuadPointInCellEnrichmentBlockEnrichment + qPoint) * cellJxWValuesEnrichmentBlockEnrichment[qPoint]; } Ni_pristine* interpolated ci's in Ni_classicalQuadratureOfPristine at quadpoints for (unsigned int qPoint = 0; qPoint < nQuadPointInCellEnrichmentBlockEnrichment; qPoint++) { ciNciNpj += classicalComponentInQuadValuesEE [numEnrichmentIdsInCell * qPoint + (iNode - dofsPerCellCFE)] * (enrichmentValuesVec.data() + (jNode - dofsPerCellCFE) * nQuadPointInCellEnrichmentBlockEnrichment + qPoint) * cellJxWValuesEnrichmentBlockEnrichment[qPoint]; } Ni_pristine* interpolated ci's in Ni_classicalQuadratureOfPristine at quadpoints for (unsigned int qPoint = 0; qPoint < nQuadPointInCellEnrichmentBlockEnrichment; qPoint++) { NpicjNcj += (enrichmentValuesVec.data() + (iNode - dofsPerCellCFE) * nQuadPointInCellEnrichmentBlockEnrichment + qPoint) * classicalComponentInQuadValuesEE [numEnrichmentIdsInCell * qPoint + (jNode - dofsPerCellCFE)] * cellJxWValuesEnrichmentBlockEnrichment[qPoint]; } interpolated ci's in Ni_classicalQuadrature of Mc = d interpolated ci's in Ni_classicalQuadrature of Mc = d for (unsigned int qPoint = 0; qPoint < nQuadPointInCellEnrichmentBlockClassical; qPoint++) { ciNcicjNcj += classicalComponentInQuadValuesEC [numEnrichmentIdsInCell * qPoint + (iNode - dofsPerCellCFE)] * classicalComponentInQuadValuesEC [numEnrichmentIdsInCell * qPoint + (jNode - dofsPerCellCFE)] * cellJxWValuesEnrichmentBlockClassical[qPoint]; } basisOverlapTmpIter += NpiNpj - NpicjNcj - ciNciNpj + ciNcicjNcj;

Here is the call graph for this function:

Here is the caller graph for this function: