A derived class of linearAlgebra::LinearSolverFunction to encapsulate the Poisson partial differential equation (PDE) discretized in a finite element (FE) basis. The Possion PDE is given as: $\nabla^2 v(\textbf{r}) = -4 \pi \rho(\textbf{r})$\f with the boundary condition on $@_fakenlv(\textbf{r})|_{\partial \Omega}=g(\textbf{r})$\f ( $\\partial Omega$\f denoting the boundary of a domain $\Omega$\f). Here $v$\f has the physical notion of a potential (e.g., Hartree potential, nuclear potential, etc.) arising due to a charge distributin $\rho$\f. More...

#include <PoissonSolverDealiiMatrixFreeFE.h>

Collaboration diagram for dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >:

[legend]

Public Types
template<typename T >
using	distributedCPUVec = basis::FEEvaluationWrapperBase::distributedCPUVec< T >
	define ValueType as the superior (bigger set) of the ValueTypeOperator and ValueTypeOperand (e.g., between double and complex<double>, complex<double> is the bigger set) More...

using	ValueType = linearAlgebra::blasLapack::scalar_type< ValueTypeOperator, ValueTypeOperand >

Public Member Functions
	PoissonSolverDealiiMatrixFreeFE (std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > > feBasisManagerField, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > feBasisDataStorageStiffnessMatrix, const std::map< std::string, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > > &feBasisDataStorageRhs, const std::map< std::string, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > & > &inpRhs, const linearAlgebra::PreconditionerType pcType, std::shared_ptr< linearAlgebra::LinAlgOpContext< memorySpace > > linAlgOpContext)
	This constructor creates an instance of a base LinearSolverFunction called PoissonSolverDealiiMatrixFreeFE. More...

	PoissonSolverDealiiMatrixFreeFE (std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > > feBasisManagerField, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > feBasisDataStorageStiffnessMatrix, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > feBasisDataStorageRhs, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > &inpRhs, const linearAlgebra::PreconditionerType pcType, std::shared_ptr< linearAlgebra::LinAlgOpContext< memorySpace > > linAlgOpContext)
	This constructor creates an instance of a base LinearSolverFunction called PoissonSolverDealiiMatrixFreeFE. More...

void	reinit (std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > > feBasisManagerField, const std::map< std::string, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > & > &inpRhs)

void	reinit (std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > > feBasisManagerField, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > &inpRhs)

	~PoissonSolverDealiiMatrixFreeFE ()=default

void	solve (const double absTolerance, const size_type maxNumberIterations)

void	getSolution (linearAlgebra::MultiVector< ValueType, memorySpaceHost > &solution)

const utils::mpi::MPIComm &	getMPIComm () const

Private Member Functions
const distributedCPUVec< ValueTypeOperand > &	getRhs () const

const distributedCPUVec< ValueType > &	getInitialGuess () const

void	setSolution (const distributedCPUVec< ValueType > &x)

void	computeRhs (distributedCPUVec< double > &rhs, const std::map< std::string, const quadrature::QuadratureValuesContainer< linearAlgebra::blasLapack::scalar_type< ValueTypeOperator, ValueTypeOperand >, memorySpaceHost > & > &inpRhs)

void	vmult (distributedCPUVec< double > &Ax, distributedCPUVec< double > &x)

void	precondition_Jacobi (distributedCPUVec< double > &dst, const distributedCPUVec< double > &src) const

void	computeDiagonalA ()

void	AX (const dealii::MatrixFree< dim, double > &matrixFreeData, distributedCPUVec< double > &dst, const distributedCPUVec< double > &src, const std::pair< size_type, size_type > &cell_range) const

void	CGsolve (const double absTolerance, const size_type maxNumberIterations, bool distributeFlag)

const linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &	getRhsDevice () const

const linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &	getInitialGuessDevice () const

void	CGsolveDevice (const double absTolerance, const size_type maxNumberIterations, bool distributeFlag)

void	computeAXDevice (linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &Ax, linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &x)

Private Attributes
size_type	d_numComponents

std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	d_feBasisManagerField

linearAlgebra::PreconditionerType	d_pcType

utils::Profiler< memorySpace >	d_p

std::shared_ptr< basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	d_feBasisManagerHomo

distributedCPUVec< ValueTypeOperator >	d_x

distributedCPUVec< ValueTypeOperator >	d_rhs

distributedCPUVec< ValueTypeOperator >	d_initial

distributedCPUVec< ValueTypeOperator >	d_diagonalA

std::shared_ptr< dealii::MatrixFree< dim, ValueTypeOperator > >	d_dealiiMatrixFree

std::shared_ptr< const dealii::DoFHandler< dim > >	d_dealiiDofHandler

const dealii::AffineConstraints< ValueTypeOperand > *	d_dealiiAffineConstraintMatrix

const dealii::AffineConstraints< ValueTypeOperand > *	d_constraintsInfo

size_type	d_num1DQuadPointsStiffnessMatrix

std::map< std::string, size_type >	d_num1DQuadPointsRhs

size_type	d_feOrder

size_type	d_dofHandlerIndex

std::map< std::string, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > >	d_feBasisDataStorageRhs

std::vector< distributedCPUVec< ValueType > >	d_nonTensorSructuredQuadeRhs

size_type	d_matrixFreeQuadCompStiffnessMatrix

std::map< dealii::CellId, size_type >	d_cellIdToCellIndexMap

std::vector< dealii::Quadrature< dim > >	d_dealiiQuadratureRuleVec

dealii::MappingQ1< dim, dim >	d_mappingDealii

utils::ConditionalOStream	pcout

linearAlgebra::MultiVector< linearAlgebra::blasLapack::scalar_type< ValueTypeOperator, ValueTypeOperand >, memorySpaceHost >	d_scratchMultiVecHost

std::shared_ptr< linearAlgebra::LinAlgOpContext< memorySpace > >	d_linAlgOpContext

std::unique_ptr< linearAlgebra::Vector< ValueTypeOperator, memorySpace > >	d_diagonalADevice

std::unique_ptr< linearAlgebra::Vector< ValueTypeOperator, memorySpace > >	d_rhsDevice

std::unique_ptr< linearAlgebra::Vector< ValueTypeOperator, memorySpace > >	d_initialDevice

linearAlgebra::Vector< ValueTypeOperator, memorySpace >	d_qvec
	define some temporary vectors for cgsolver device More...

linearAlgebra::Vector< ValueTypeOperator, memorySpace >	d_rvec

linearAlgebra::Vector< ValueTypeOperator, memorySpace >	d_dvec

Detailed Description

template<typename ValueTypeOperator, typename ValueTypeOperand, utils::MemorySpace memorySpace, size_type dim>
class dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >

A derived class of linearAlgebra::LinearSolverFunction to encapsulate the Poisson partial differential equation (PDE) discretized in a finite element (FE) basis. The Possion PDE is given as: $\nabla^2 v(\textbf{r}) = -4 \pi \rho(\textbf{r})$\f with the boundary condition on $@_fakenlv(\textbf{r})|_{\partial \Omega}=g(\textbf{r})$\f ( $\\partial Omega$\f denoting the boundary of a domain $\Omega$\f). Here $v$\f has the physical notion of a potential (e.g., Hartree potential, nuclear potential, etc.) arising due to a charge distributin $\rho$\f.

Template Parameters

ValueTypeOperator	The datatype (float, double, complex<double>, etc.) for the underlying operator
ValueTypeOperand	The datatype (float, double, complex<double>, etc.) of the vector, matrices, etc. on which the operator will act
memorySpace	The meory sapce (HOST, DEVICE, HOST_PINNES, etc.) in which the data of the operator and its operands reside
dim	Dimension of the Poisson problem

Member Typedef Documentation

◆ distributedCPUVec

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

template<typename T >

using dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::distributedCPUVec = basis::FEEvaluationWrapperBase::distributedCPUVec<T>

define ValueType as the superior (bigger set) of the ValueTypeOperator and ValueTypeOperand (e.g., between double and complex<double>, complex<double> is the bigger set)

◆ ValueType

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

using dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::ValueType = linearAlgebra::blasLapack::scalar_type<ValueTypeOperator, ValueTypeOperand>

Constructor & Destructor Documentation

◆ PoissonSolverDealiiMatrixFreeFE() [1/2]

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

dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::PoissonSolverDealiiMatrixFreeFE	(	std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	feBasisManagerField,
		std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > >	feBasisDataStorageStiffnessMatrix,
		const std::map< std::string, std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > > > &	feBasisDataStorageRhs,
		const std::map< std::string, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > & > &	inpRhs,
		const linearAlgebra::PreconditionerType	pcType,
		std::shared_ptr< linearAlgebra::LinAlgOpContext< memorySpace > >	linAlgOpContext
	)

This constructor creates an instance of a base LinearSolverFunction called PoissonSolverDealiiMatrixFreeFE.

Here is the call graph for this function:

◆ PoissonSolverDealiiMatrixFreeFE() [2/2]

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

dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::PoissonSolverDealiiMatrixFreeFE	(	std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	feBasisManagerField,
		std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > >	feBasisDataStorageStiffnessMatrix,
		std::shared_ptr< const basis::FEBasisDataStorage< ValueTypeOperator, memorySpaceHost > >	feBasisDataStorageRhs,
		const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > &	inpRhs,
		const linearAlgebra::PreconditionerType	pcType,
		std::shared_ptr< linearAlgebra::LinAlgOpContext< memorySpace > >	linAlgOpContext
	)

This constructor creates an instance of a base LinearSolverFunction called PoissonSolverDealiiMatrixFreeFE.

◆ ~PoissonSolverDealiiMatrixFreeFE()

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

dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::~PoissonSolverDealiiMatrixFreeFE ( )

default

Member Function Documentation

◆ AX()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::AX	(	const dealii::MatrixFree< dim, double > &	matrixFreeData,
		distributedCPUVec< double > &	dst,
		const distributedCPUVec< double > &	src,
		const std::pair< size_type, size_type > &	cell_range
	)		const

private

Here is the call graph for this function:

◆ CGsolve()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::CGsolve	(	const double	absTolerance,
		const size_type	maxNumberIterations,
		bool	distributeFlag
	)

private

Here is the call graph for this function:

◆ CGsolveDevice()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::CGsolveDevice	(	const double	absTolerance,
		const size_type	maxNumberIterations,
		bool	distributeFlag
	)

private

Here is the call graph for this function:

◆ computeAXDevice()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::computeAXDevice	(	linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &	Ax,
		linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > &	x
	)

private

Here is the call graph for this function:

◆ computeDiagonalA()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::computeDiagonalA

private

Here is the call graph for this function:

Here is the caller graph for this function:

◆ computeRhs()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::computeRhs	(	distributedCPUVec< double > &	rhs,
		const std::map< std::string, const quadrature::QuadratureValuesContainer< linearAlgebra::blasLapack::scalar_type< ValueTypeOperator, ValueTypeOperand >, memorySpaceHost > & > &	inpRhs
	)

private

Here is the call graph for this function:

◆ getInitialGuess()

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

const basis::FEEvaluationWrapperBase::distributedCPUVec< linearAlgebra::blasLapack::scalar_type< ValueTypeOperator, ValueTypeOperand > > & dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getInitialGuess

private

◆ getInitialGuessDevice()

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

const linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > & dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getInitialGuessDevice

private

◆ getMPIComm()

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

const utils::mpi::MPIComm & dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getMPIComm

Here is the caller graph for this function:

◆ getRhs()

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

const basis::FEEvaluationWrapperBase::distributedCPUVec< ValueTypeOperand > & dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getRhs

private

◆ getRhsDevice()

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

const linearAlgebra::Vector< ValueTypeOperator, utils::MemorySpace::DEVICE > & dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getRhsDevice

private

◆ getSolution()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::getSolution ( linearAlgebra::MultiVector< ValueType, memorySpaceHost > & solution )

Here is the call graph for this function:

◆ precondition_Jacobi()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::precondition_Jacobi	(	distributedCPUVec< double > &	dst,
		const distributedCPUVec< double > &	src
	)		const

private

◆ reinit() [1/2]

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::reinit	(	std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	feBasisManagerField,
		const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > &	inpRhs
	)

◆ reinit() [2/2]

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::reinit	(	std::shared_ptr< const basis::FEBasisManager< ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim > >	feBasisManagerField,
		const std::map< std::string, const quadrature::QuadratureValuesContainer< ValueType, memorySpaceHost > & > &	inpRhs
	)

Here is the call graph for this function:

Here is the caller graph for this function:

◆ setSolution()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::setSolution ( const distributedCPUVec< ValueType > & x )

private

◆ solve()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::solve	(	const double	absTolerance,
		const size_type	maxNumberIterations
	)

◆ vmult()

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

void dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::vmult	(	distributedCPUVec< double > &	Ax,
		distributedCPUVec< double > &	x
	)

private

Member Data Documentation

◆ d_cellIdToCellIndexMap

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

std::map<dealii::CellId, size_type> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_cellIdToCellIndexMap

private

◆ d_constraintsInfo

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

const dealii::AffineConstraints<ValueTypeOperand>* dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_constraintsInfo

private

◆ d_dealiiAffineConstraintMatrix

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

const dealii::AffineConstraints<ValueTypeOperand>* dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dealiiAffineConstraintMatrix

private

◆ d_dealiiDofHandler

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

std::shared_ptr<const dealii::DoFHandler<dim> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dealiiDofHandler

private

◆ d_dealiiMatrixFree

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

std::shared_ptr<dealii::MatrixFree<dim, ValueTypeOperator> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dealiiMatrixFree

private

◆ d_dealiiQuadratureRuleVec

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

std::vector<dealii::Quadrature<dim> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dealiiQuadratureRuleVec

private

◆ d_diagonalA

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

distributedCPUVec<ValueTypeOperator> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_diagonalA

private

◆ d_diagonalADevice

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

std::unique_ptr<linearAlgebra::Vector<ValueTypeOperator, memorySpace> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_diagonalADevice

private

◆ d_dofHandlerIndex

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

size_type dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dofHandlerIndex

private

◆ d_dvec

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

linearAlgebra::Vector<ValueTypeOperator, memorySpace> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_dvec

private

◆ d_feBasisDataStorageRhs

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

std::map< std::string, std::shared_ptr< const basis::FEBasisDataStorage<ValueTypeOperator, memorySpaceHost> > > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_feBasisDataStorageRhs

private

◆ d_feBasisManagerField

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

std::shared_ptr<const basis::FEBasisManager<ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_feBasisManagerField

private

◆ d_feBasisManagerHomo

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

std::shared_ptr<basis::FEBasisManager<ValueTypeOperand, ValueTypeOperator, memorySpaceHost, dim> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_feBasisManagerHomo

private

◆ d_feOrder

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

size_type dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_feOrder

private

◆ d_initial

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

distributedCPUVec<ValueTypeOperator> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_initial

private

◆ d_initialDevice

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

std::unique_ptr<linearAlgebra::Vector<ValueTypeOperator, memorySpace> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_initialDevice

private

◆ d_linAlgOpContext

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

std::shared_ptr<linearAlgebra::LinAlgOpContext<memorySpace> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_linAlgOpContext

private

◆ d_mappingDealii

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

dealii::MappingQ1<dim, dim> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_mappingDealii

private

◆ d_matrixFreeQuadCompStiffnessMatrix

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

size_type dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_matrixFreeQuadCompStiffnessMatrix

private

◆ d_nonTensorSructuredQuadeRhs

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

std::vector<distributedCPUVec<ValueType> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_nonTensorSructuredQuadeRhs

private

◆ d_num1DQuadPointsRhs

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

std::map<std::string, size_type> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_num1DQuadPointsRhs

private

◆ d_num1DQuadPointsStiffnessMatrix

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

size_type dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_num1DQuadPointsStiffnessMatrix

private

◆ d_numComponents

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

size_type dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_numComponents

private

◆ d_p

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

utils::Profiler<memorySpace> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_p

private

◆ d_pcType

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

linearAlgebra::PreconditionerType dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_pcType

private

◆ d_qvec

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

linearAlgebra::Vector<ValueTypeOperator, memorySpace> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_qvec

private

define some temporary vectors for cgsolver device

◆ d_rhs

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

distributedCPUVec<ValueTypeOperator> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_rhs

private

◆ d_rhsDevice

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

std::unique_ptr<linearAlgebra::Vector<ValueTypeOperator, memorySpace> > dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_rhsDevice

private

◆ d_rvec

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

linearAlgebra::Vector<ValueTypeOperator, memorySpace> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_rvec

private

◆ d_scratchMultiVecHost

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

linearAlgebra::MultiVector< linearAlgebra::blasLapack::scalar_type<ValueTypeOperator, ValueTypeOperand>, memorySpaceHost> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_scratchMultiVecHost

private

◆ d_x

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

distributedCPUVec<ValueTypeOperator> dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::d_x

private

◆ pcout

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

utils::ConditionalOStream dftefe::electrostatics::PoissonSolverDealiiMatrixFreeFE< ValueTypeOperator, ValueTypeOperand, memorySpace, dim >::pcout

private

The documentation for this class was generated from the following files:

src/electrostatics/PoissonSolverDealiiMatrixFreeFE.h
src/electrostatics/PoissonSolverDealiiMatrixFreeFE.t.cpp

Public Types

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ distributedCPUVec

◆ ValueType

Constructor & Destructor Documentation

◆ PoissonSolverDealiiMatrixFreeFE() [1/2]

◆ PoissonSolverDealiiMatrixFreeFE() [2/2]

◆ ~PoissonSolverDealiiMatrixFreeFE()

Member Function Documentation

◆ AX()

◆ CGsolve()

◆ CGsolveDevice()

◆ computeAXDevice()

◆ computeDiagonalA()

◆ computeRhs()

◆ getInitialGuess()

◆ getInitialGuessDevice()

◆ getMPIComm()

◆ getRhs()

◆ getRhsDevice()

◆ getSolution()

◆ precondition_Jacobi()

◆ reinit() [1/2]

◆ reinit() [2/2]

◆ setSolution()

◆ solve()

◆ vmult()

Member Data Documentation

◆ d_cellIdToCellIndexMap

◆ d_constraintsInfo

◆ d_dealiiAffineConstraintMatrix

◆ d_dealiiDofHandler

◆ d_dealiiMatrixFree

◆ d_dealiiQuadratureRuleVec

◆ d_diagonalA

◆ d_diagonalADevice

◆ d_dofHandlerIndex

◆ d_dvec

◆ d_feBasisDataStorageRhs

◆ d_feBasisManagerField

◆ d_feBasisManagerHomo

◆ d_feOrder

◆ d_initial

◆ d_initialDevice

◆ d_linAlgOpContext

◆ d_mappingDealii

◆ d_matrixFreeQuadCompStiffnessMatrix

◆ d_nonTensorSructuredQuadeRhs

◆ d_num1DQuadPointsRhs

◆ d_num1DQuadPointsStiffnessMatrix

◆ d_numComponents

◆ d_p

◆ d_pcType

◆ d_qvec

◆ d_rhs

◆ d_rhsDevice

◆ d_rvec

◆ d_scratchMultiVecHost

◆ d_x

◆ pcout