Templated bilinear form class, cf. bilinearform.?pp. More...

#include <tbilinearform.hpp>

Inheritance diagram for mfem::TBilinearForm< meshType, solFESpace, IR, IntegratorType, solVecLayout_t, complex_t, real_t, impl_traits_t >:

[legend]

Collaboration diagram for mfem::TBilinearForm< meshType, solFESpace, IR, IntegratorType, solVecLayout_t, complex_t, real_t, impl_traits_t >:

[legend]

Classes
struct	S_spec
	Contains matrix sizes, type of kernel (ElementMatrix is templated on a kernel, e.g. ElementMatrix::Compute may be AssembleGradGrad()). More...

struct	T_result

Public Types
typedef impl_traits_t	impl_traits_type

Public Types inherited from mfem::Operator
enum	DiagonalPolicy { DIAG_ZERO, DIAG_ONE, DIAG_KEEP }
	Defines operator diagonal policy upon elimination of rows and/or columns. More...

enum	Type { ANY_TYPE, MFEM_SPARSEMAT, Hypre_ParCSR, PETSC_MATAIJ, PETSC_MATIS, PETSC_MATSHELL, PETSC_MATNEST, PETSC_MATHYPRE, PETSC_MATGENERIC, Complex_Operator, MFEM_ComplexSparseMat, Complex_Hypre_ParCSR }
	Enumeration defining IDs for some classes derived from Operator. More...

Public Member Functions
	TBilinearForm (const IntegratorType &integ, const FiniteElementSpace &sol_fes)

virtual	~TBilinearForm ()

virtual const Operator *	GetProlongation () const
	Get the input finite element space prolongation matrix. More...

virtual const Operator *	GetRestriction () const
	Get the input finite element space restriction matrix. More...

virtual void	Mult (const Vector &x, Vector &y) const
	Operator application: `y=A(x)`. More...

void	MultUnassembled (const Vector &x, Vector &y) const

void	Assemble ()
	Partial assembly of quadrature point data. More...

MFEM_ALWAYS_INLINE void	ElementAddMultAssembled (int el, solFieldEval &solFEval) const

void	MultAssembled (const Vector &x, Vector &y) const

void	TestElementwiseExtractAssemble (const Vector &x, Vector &y) const

void	SerializeNodes (Vector &sNodes) const

void	AssembleFromSerializedNodes (const Vector &sNodes)
	Partial assembly from "serialized" nodes. More...

void	Serialize (const Vector &x, Vector &sx) const

void	MultAssembledSerialized (const Vector &sx, Vector &sy) const
	serialized vector sx –> serialized vector 'sy' More...

void	AssembleMatrix (SparseMatrix &M) const
	Assemble the operator in a SparseMatrix. More...

void	AssembleMatrix (DenseTensor &M) const
	Assemble element matrices and store them as a DenseTensor object. More...

void	AssembleBilinearForm (BilinearForm &a) const
	Assemble element matrices and add them to the bilinear form. More...

void	AddMult (DenseTensor &M, const Vector &x, Vector &y) const
	Multiplication using assembled element matrices stored as a DenseTensor. More...

Public Member Functions inherited from mfem::Operator
void	InitTVectors (const Operator Po, const Operator Ri, const Operator *Pi, Vector &x, Vector &b, Vector &X, Vector &B) const
	Initializes memory for true vectors of linear system. More...

	Operator (int s=0)
	Construct a square Operator with given size s (default 0). More...

	Operator (int h, int w)
	Construct an Operator with the given height (output size) and width (input size). More...

int	Height () const
	Get the height (size of output) of the Operator. Synonym with NumRows(). More...

int	NumRows () const
	Get the number of rows (size of output) of the Operator. Synonym with Height(). More...

int	Width () const
	Get the width (size of input) of the Operator. Synonym with NumCols(). More...

int	NumCols () const
	Get the number of columns (size of input) of the Operator. Synonym with Width(). More...

virtual MemoryClass	GetMemoryClass () const
	Return the MemoryClass preferred by the Operator. More...

virtual void	MultTranspose (const Vector &x, Vector &y) const
	Action of the transpose operator: `y=A^t(x)`. The default behavior in class Operator is to generate an error. More...

virtual Operator &	GetGradient (const Vector &x) const
	Evaluate the gradient operator at the point x. The default behavior in class Operator is to generate an error. More...

virtual void	AssembleDiagonal (Vector &diag) const
	Computes the diagonal entries into diag. Typically, this operation only makes sense for linear Operators. In some cases, only an approximation of the diagonal is computed. More...

virtual const Operator *	GetOutputProlongation () const
	Prolongation operator from linear algebra (linear system) vectors, to output vectors for the operator. `NULL` means identity. More...

virtual const Operator *	GetOutputRestrictionTranspose () const
	Transpose of GetOutputRestriction, directly available in this form to facilitate matrix-free RAP-type operators. More...

virtual const Operator *	GetOutputRestriction () const
	Restriction operator from output vectors for the operator to linear algebra (linear system) vectors. `NULL` means identity. More...

void	FormLinearSystem (const Array< int > &ess_tdof_list, Vector &x, Vector &b, Operator *&A, Vector &X, Vector &B, int copy_interior=0)
	Form a constrained linear system using a matrix-free approach. More...

void	FormRectangularLinearSystem (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Vector &x, Vector &b, Operator *&A, Vector &X, Vector &B)
	Form a column-constrained linear system using a matrix-free approach. More...

virtual void	RecoverFEMSolution (const Vector &X, const Vector &b, Vector &x)
	Reconstruct a solution vector x (e.g. a GridFunction) from the solution X of a constrained linear system obtained from Operator::FormLinearSystem() or Operator::FormRectangularLinearSystem(). More...

void	FormSystemOperator (const Array< int > &ess_tdof_list, Operator *&A)
	Return in A a parallel (on truedofs) version of this square operator. More...

void	FormRectangularSystemOperator (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, Operator *&A)
	Return in A a parallel (on truedofs) version of this rectangular operator (including constraints). More...

void	FormDiscreteOperator (Operator *&A)
	Return in A a parallel (on truedofs) version of this rectangular operator. More...

void	PrintMatlab (std::ostream &out, int n, int m=0) const
	Prints operator with input size n and output size m in Matlab format. More...

virtual void	PrintMatlab (std::ostream &out) const
	Prints operator in Matlab format. More...

virtual	~Operator ()
	Virtual destructor. More...

Type	GetType () const
	Return the type ID of the Operator class. More...

Protected Types
typedef complex_t	complex_type

typedef real_t	real_type

typedef meshType::FE_type	meshFE_type

typedef ShapeEvaluator < meshFE_type, IR, real_t >	meshShapeEval

typedef solFESpace::FE_type	solFE_type

typedef ShapeEvaluator < solFE_type, IR, real_t >	solShapeEval

typedef solVecLayout_t	solVecLayout_type

typedef impl_traits_t::vcomplex_t	vcomplex_t

typedef impl_traits_t::vreal_t	vreal_t

typedef kernel_t::template CoefficientEval< IR, coeff_t, impl_traits_t >::Type	coeff_eval_t

typedef TElementTransformation < meshType, IR, real_t >	Trans_t

typedef FieldEvaluator < solFESpace, solVecLayout_t, IR, complex_t, real_t >	solFieldEval

IntegratorType defines several internal types
typedef IntegratorType	integ_t

typedef integ_t::coefficient_type	coeff_t
	coeff_t might be TConstantCoefficient or TFunctionCoefficient, for example More...

typedef integ_t::template kernel< sdim, dim, vcomplex_t > ::type	kernel_t
	kernel_t may be TDiffusionKernel or TMassKernel More...

typedef kernel_t::template p_asm_data< qpts >::type	p_assembled_t
	p_assembled_t is something like a TTensor or TMatrix for partial assembly More...

typedef kernel_t::template f_asm_data< qpts >::type	f_assembled_t
	f_assembled_t is something like a TTensor or TMatrix for full assembly More...

Protected Attributes
meshType	mesh

meshShapeEval	meshEval

solFE_type	sol_fe

solShapeEval	solEval

solFESpace	solFES

solVecLayout_t	solVecLayout

IR	int_rule

coeff_t	coeff

Memory< p_assembled_t >	assembled_data

const FiniteElementSpace &	in_fes

Protected Attributes inherited from mfem::Operator
int	height
	Dimension of the output / number of rows in the matrix. More...

int	width
	Dimension of the input / number of columns in the matrix. More...

Static Protected Attributes
static const int	dim = meshType::dim

static const int	sdim = meshType::space_dim

static const int	dofs = solFE_type::dofs

static const int	vdim = solVecLayout_t::vec_dim

static const int	qpts = IR::qpts

static const int	AB = impl_traits_t::align_bytes

static const int	SS = impl_traits_t::simd_size

static const int	BE = impl_traits_t::batch_size

static const int	TE = SS*BE

Additional Inherited Members
Protected Member Functions inherited from mfem::Operator
void	FormConstrainedSystemOperator (const Array< int > &ess_tdof_list, ConstrainedOperator *&Aout)
	see FormSystemOperator() More...

void	FormRectangularConstrainedSystemOperator (const Array< int > &trial_tdof_list, const Array< int > &test_tdof_list, RectangularConstrainedOperator *&Aout)
	see FormRectangularSystemOperator() More...

Operator *	SetupRAP (const Operator Pi, const Operator Po)
	Returns RAP Operator of this, using input/output Prolongation matrices Pi corresponds to "P", Po corresponds to "Rt". More...

Detailed Description

template<typename meshType, typename solFESpace, typename IR, typename IntegratorType, typename solVecLayout_t = ScalarLayout, typename complex_t = double, typename real_t = double, typename impl_traits_t = AutoSIMDTraits<complex_t,real_t>>
class mfem::TBilinearForm< meshType, solFESpace, IR, IntegratorType, solVecLayout_t, complex_t, real_t, impl_traits_t >

Templated bilinear form class, cf. bilinearform.?pp.

Template Parameters

meshType	typically TMesh, which is templated on FE type
solFESpace	eg. H1_FiniteElementSpace
IR	integration rule, typically TIntegrationRule, which is further templated on element geometry
IntegratorType	typically a TIntegrator, which is templated on a kernel, eg. TDiffusionKernel or TMassKernel. This describes what actual problem you solve.
solVecLayout_t	describes how degrees of freedom are laid out, scalar or vector, column/row major, etc.
complex_t	data type for solution dofs
real_t	data type for mesh nodes, solution basis, and mesh basis

Definition at line 46 of file tbilinearform.hpp.