Detailed Description

template<class T>
class gismo::gsGenericAssembler< T >

Assembles mass and stiffness matrices and right-hand sides on a given domain.

Inheritance diagram for gsGenericAssembler< T >:

Collaboration diagram for gsGenericAssembler< T >:

Public Member Functions
const std::vector< gsMatrix< T > > &	allFixedDofs () const
	Returns all the Dirichlet values (if applicable)

virtual void	assemble ()
	Main assemble routine, to be implemented in derived classes.

virtual void	assemble (const gsMultiPatch< T > &curSolution)
	Main non-linear assemble routine with input from current solution.

const gsSparseMatrix< T > &	assembleDG (const boundaryInterface &iFace, bool refresh=true)
	Assemble dG interface terms.

const gsSparseMatrix< T > &	assembleMass (const index_t patchIndex=-1, bool refresh=true)
	Mass assembly routine.

const gsSparseMatrix< T > &	assembleMass2 ()
	Mass assembly routine.

const gsMatrix< T > &	assembleMoments (const gsFunction< T > &func, index_t patchIndex=-1, bool refresh=true)
	Moments assembly routine.

const gsSparseMatrix< T > &	assembleNeumann (const boundary_condition< T > &bc, bool refresh=true)
	Assemble Neumann boundary terms.

const gsSparseMatrix< T > &	assembleNitsche (const boundary_condition< T > &bc, bool refresh=true)
	Assemble Nitsche terms for weakly imposing Dirichlet conditions.

const gsSparseMatrix< T > &	assembleStiffness (const index_t patchIndex=-1, const bool refresh=true)
	Stiffness assembly routine.

bool	check ()
	checks for consistency and legal values of the stored members.

virtual gsAssembler *	clone () const
	Clone this Assembler, making a deep copy.

void	computeDirichletDofs (short_t unk=0)
	Triggers computation of the Dirichlet dofs.

virtual void	constructSolution (const gsMatrix< T > &solVector, gsMultiPatch< T > &result, const gsVector< index_t > &unknowns) const
	Construct solution from computed solution vector for a set of unknowns. The result is a vectorfield, where each component is given the corresponding entry of.

virtual void	constructSolution (const gsMatrix< T > &solVector, gsMultiPatch< T > &result, short_t unk=0) const
	Construct solution from computed solution vector for a single unknows.

virtual gsAssembler *	create () const
	Create an empty Assembler of the derived type and return a pointer to it. Call the initialize functions to set the members.

void	finalize ()
	finishes the assembling of the system matrix, i.e. calls its .makeCompressed() method.

const gsMatrix< T > &	fixedDofs (short_t unk=0) const
	Returns the Dirichlet values for a unknown (if applicable)

gsSparseMatrix< T >::fullView	fullMatrix ()
	Returns an expression of the "full" assembled sparse matrix.

const gsSparseMatrix< T >::constFullView	fullMatrix () const
	Returns an expression of the "full" assembled sparse matrix.

	gsGenericAssembler (const gsMultiPatch< T > &patches, const gsMultiBasis< T > &bases, const gsOptionList &opt=defaultOptions(), const gsBoundaryConditions< T > *bc=NULL)
	Constructor.

void	homogenizeFixedDofs (short_t unk=0)
	Sets any Dirichlet values to homogeneous (if applicable)

void	initialize (const gsPde< T > &pde, const gsBasisRefs< T > &basis, const gsOptionList &opt=defaultOptions())
	Intitialize function for, sets data fields using the pde, a vector of bases and assembler options.

void	initialize (const gsPde< T > &pde, const gsMultiBasis< T > &bases, const gsOptionList &opt=defaultOptions())
	Intitialize function for, sets data fields using the pde, a multi-basis and assembler options.

void	initialize (const gsPde< T > &pde, const gsStdVectorRef< gsMultiBasis< T > > &bases, const gsOptionList &opt=defaultOptions())
	Intitialize function for, sets data fields using the pde, a vector of multi-basis and assembler options.

void	initialize (typename gsPde< T >::Ptr pde, const gsStdVectorRef< gsMultiBasis< T > > &bases, const gsOptionList &opt=defaultOptions())
	Intitialize function for, sets data fields using the pde, a vector of multi-basis and assembler options.

const gsSparseMatrix< T > &	matrix () const
	Returns the left-hand global matrix.

gsMultiBasis< T > &	multiBasis (index_t k=0)
	Return the multi-basis. Note: if the basis is altered, then refresh() should be called.

const gsMultiBasis< T > &	multiBasis (index_t k=0) const
	Return the multi-basis.

index_t	numColNz () const
	Provides an estimation of the number of non-zero matrix entries per column. This value can be used for sparse matrix memory allocation.

index_t	numDofs () const
	Returns the number of (free) degrees of freedom.

size_t	numMultiBasis () const
	Returns the number of multi-bases.

const gsMultiPatch< T > &	patches () const
	Return the multipatch.

const gsPde< T > &	pde () const
	Return the Pde.

void	penalizeDirichletDofs (short_t unk=0)

T	penalty (index_t k) const
	Penalty constant for patch k, used for Nitsche and / Discontinuous Galerkin methods.

template<class ElementVisitor >
void	push ()
	Iterates over all elements of the domain and applies the ElementVisitor.

template<class BElementVisitor >
void	push (const bcContainer &BCs)
	Iterates over all elements of the boundaries BCs and applies the BElementVisitor.

template<class BElementVisitor >
void	push (const BElementVisitor &visitor, const boundary_condition< T > &BC)
	Applies the BElementVisitor to the boundary condition BC.

template<class ElementVisitor >
void	push (const ElementVisitor &visitor)
	Iterates over all elements of the domain and applies the ElementVisitor.

template<class InterfaceVisitor >
void	pushInterface ()
	Iterates over all elements of interfaces and applies the InterfaceVisitor.

void	refresh ()
	Refreshes the sparse system (deletes assembled matrices and vectors)

void	refresh (gsDofMapper mapper)
	Refreshes the sparse system based on given dof mapper.

const gsMatrix< T > &	rhs () const
	Returns the left-hand side vector(s) ( multiple right hand sides possible )

void	setFixedDofs (const gsMatrix< T > &coefMatrix, short_t unk=0, size_t patch=0)
	the user can manually set the dirichlet Dofs for a given patch and unknown, based on the Basis coefficients

void	setFixedDofVector (gsMatrix< T > vals, short_t unk=0)
	the user can manually set the dirichlet Dofs for a given patch and unknown.

void	setSparseSystem (gsSparseSystem< T > &sys)
	Swaps the actual sparse system with the given one.

const gsSparseSystem< T > &	system () const
	Returns the left-hand global matrix.

virtual void	updateSolution (const gsMatrix< T > &solVector, gsMultiPatch< T > &result, T theta=(T)(1)) const
	Update solution by adding the computed solution vector to the current solution specified by.

Static Public Member Functions
static gsOptionList	defaultOptions ()
	Returns the list of default options for assembly.

Protected Member Functions
template<class ElementVisitor >
void	apply (ElementVisitor &visitor, size_t patchIndex=0, boxSide side=boundary::none)
	Generic assembly routine for volume or boundary integrals.

template<class InterfaceVisitor >
void	apply (InterfaceVisitor &visitor, const boundaryInterface &bi)
	Generic assembly routine for patch-interface integrals.

void	computeDirichletDofsIntpl (const gsDofMapper &mapper, const gsMultiBasis< T > &mbasis, const short_t unk_=0)
	calculates the values of the eliminated dofs based on Interpolation.

void	computeDirichletDofsL2Proj (const gsDofMapper &mapper, const gsMultiBasis< T > &mbasis, const short_t unk_=0)
	calculates the values of the eliminated dofs based on L2 Projection.

void	scalarProblemGalerkinRefresh ()
	A prototype of the refresh function for a "standard" scalar problem. Creats one mapper based on the set options and initializes the sparse system (without allocating memory.

Private Attributes
std::vector< gsMultiBasis< T > >	m_bases

std::vector< gsMatrix< T > >	m_ddof

gsOptionList	m_options
	Options.

memory::shared_ptr< gsPde< T > >	m_pde_ptr

gsSparseSystem< T >	m_system
	Global sparse linear system.

Constructor & Destructor Documentation

◆ gsGenericAssembler()

template<class T >

gsGenericAssembler	(	const gsMultiPatch< T > &	patches,
		const gsMultiBasis< T > &	bases,
		const gsOptionList &	opt = `defaultOptions()`,
		const gsBoundaryConditions< T > *	bc = `NULL`
	)

inline

Constructor.

Parameters

patches	The geometry as gsMultiPatch
bases	The bases as gsMultiBasis
opt	The assembler options
bc	The boundary conditons (used for setting up gsDofMapper)

Member Function Documentation

◆ apply()

template<class T >

template<class ElementVisitor >

void apply	(	ElementVisitor &	visitor,
		size_t	patchIndex = `0`,
		boxSide	side = `boundary::none`
	)

protectedinherited

Generic assembly routine for volume or boundary integrals.

Parameters

[in]	visitor	The visitor for the boundary or volume integral
[in]	patchIndex	The considered patch
[in]	side	The considered boundary side, only necessary for boundary integrals.

◆ assembleDG()

template<class T >

const gsSparseMatrix< T > & assembleDG	(	const boundaryInterface &	iFace,
		bool	refresh = `true`
	)

Assemble dG interface terms.

Parameters

iFace	The interface for which assembling is done
refresh	Calles member refresh See gsVisiorDg for possible options

◆ assembleMass()

template<class T >

const gsSparseMatrix< T > & assembleMass	(	const index_t	patchIndex = `-1`,
		bool	refresh = `true`
	)

Mass assembly routine.

Parameters

patchIndex	If non-negative, only assemble matrix for specified patch
refresh	Calles member refresh

◆ assembleMass2()

template<class T >

const gsSparseMatrix< T > & assembleMass2 ( )

inline

Mass assembly routine.

This routine uses the expression assembler and always overwrites already assembled matrices.

◆ assembleMoments()

template<class T >

const gsMatrix< T > & assembleMoments	(	const gsFunction< T > &	func,
		index_t	patchIndex = `-1`,
		bool	refresh = `true`
	)

Moments assembly routine.

Parameters

func	Right-hand-side (source) function
patchIndex	If non-negative, only assemble matrix for specified patch
refresh	Calles member refresh

◆ assembleNeumann()

template<class T >

const gsSparseMatrix< T > & assembleNeumann	(	const boundary_condition< T > &	bc,
		bool	refresh = `true`
	)

Assemble Neumann boundary terms.

Parameters

bc	The boundary condition for which assembling is done
refresh	Calles member refresh

◆ assembleNitsche()

template<class T >

const gsSparseMatrix< T > & assembleNitsche	(	const boundary_condition< T > &	bc,
		bool	refresh = `true`
	)

Assemble Nitsche terms for weakly imposing Dirichlet conditions.

Parameters

bc	The boundary condition for which assembling is done
refresh	Calles member refresh

See gsVisiorNitsche for possible options

◆ assembleStiffness()

template<class T >

const gsSparseMatrix< T > & assembleStiffness	(	const index_t	patchIndex = `-1`,
		const bool	refresh = `true`
	)

Stiffness assembly routine.

Parameters

patchIndex	If non-negative, only assemble matrix for specified patch
refresh	Calles member refresh

◆ computeDirichletDofs()

template<class T >

void computeDirichletDofs ( short_t unk = 0 )

inherited

Triggers computation of the Dirichlet dofs.

Parameters

[in] unk the considered unknown

◆ computeDirichletDofsIntpl()

template<class T >

void computeDirichletDofsIntpl	(	const gsDofMapper &	mapper,
		const gsMultiBasis< T > &	mbasis,
		const short_t	unk_ = `0`
	)

protectedinherited

calculates the values of the eliminated dofs based on Interpolation.

Parameters

[in]	mapper	the dofMapper for the considered unknown
[in]	mbasis	the multipabasis for the considered unknown
[in]	unk_	the considered unknown

◆ computeDirichletDofsL2Proj()

template<class T >

void computeDirichletDofsL2Proj	(	const gsDofMapper &	mapper,
		const gsMultiBasis< T > &	mbasis,
		const short_t	unk_ = `0`
	)

protectedinherited

calculates the values of the eliminated dofs based on L2 Projection.

Parameters

[in]	mapper	the dofMapper for the considered unknown
[in]	mbasis	the multipabasis for the considered unknown
[in]	unk_	the considered unknown

◆ constructSolution() [1/2]

template<class T >

void constructSolution	(	const gsMatrix< T > &	solVector,
		gsMultiPatch< T > &	result,
		const gsVector< index_t > &	unknowns
	)		const

virtualinherited

Construct solution from computed solution vector for a set of unknowns. The result is a vectorfield, where each component is given the corresponding entry of.

unknowns. This method assumes that the specified unknowns have the: same basis.

Parameters

[in]	solVector	the solution vector obtained from the linear system
[out]	result	the solution seen as vectorfield in form of a gsMultiBasis
[in]	unknowns	the considered vector of unknowns

Reimplemented in gsBaseAssembler< T >, gsElasticityAssembler< T >, gsElTimeIntegrator< T >, gsNsAssembler< T >, gsNsTimeIntegrator< T >, and gsBiharmonicAssembler< T, bhVisitor >.

◆ constructSolution() [2/2]

template<class T >

void constructSolution	(	const gsMatrix< T > &	solVector,
		gsMultiPatch< T > &	result,
		short_t	unk = `0`
	)		const

virtualinherited

Construct solution from computed solution vector for a single unknows.

Parameters

[in]	solVector	the solution vector obtained from the linear system
[out]	result	the solution in form of a gsMultiBasis
[in]	unk	the considered unknown

Reimplemented in gsBaseAssembler< T >, gsElasticityAssembler< T >, gsElTimeIntegrator< T >, gsNsAssembler< T >, gsNsTimeIntegrator< T >, and gsBiharmonicAssembler< T, bhVisitor >.

◆ fixedDofs()

template<class T >

const gsMatrix< T > & fixedDofs ( short_t unk = 0 ) const

inlineinherited

Returns the Dirichlet values for a unknown (if applicable)

Parameters

[in] unk the considered unknown

◆ fullMatrix() [1/2]

template<class T >

gsSparseMatrix< T >::fullView fullMatrix ( )

inline

Returns an expression of the "full" assembled sparse matrix.

Note that matrix() might return a lower diagonal matrix, if we exploit possible symmetry during assembly (check: m_matrix.symmetry() == true )

◆ fullMatrix() [2/2]

template<class T >

const gsSparseMatrix< T >::constFullView fullMatrix ( ) const

inline

Returns an expression of the "full" assembled sparse matrix.

Note that matrix() might return a lower diagonal matrix, if we exploit possible symmetry during assembly (check: m_matrix.symmetry() == true )

◆ homogenizeFixedDofs()

template<class T >

void homogenizeFixedDofs ( short_t unk = 0 )

inlineinherited

Sets any Dirichlet values to homogeneous (if applicable)

Parameters

[in] unk the considered unknown

◆ penalizeDirichletDofs()

template<class T >

void penalizeDirichletDofs ( short_t unk = 0 )

inherited

Enforce Dirichlet boundary conditions by diagonal penalization

Parameters

[in] unk the considered unknown

◆ setFixedDofs()

template<class T >

void setFixedDofs	(	const gsMatrix< T > &	coefMatrix,
		short_t	unk = `0`,
		size_t	patch = `0`
	)

inherited

the user can manually set the dirichlet Dofs for a given patch and unknown, based on the Basis coefficients

Parameters

[in]	coefMatrix	the coefficients of the function
[in]	unk	the consideren unknown
[in]	patch	the patch index

◆ setFixedDofVector()

template<class T >

void setFixedDofVector	(	gsMatrix< T >	vals,
		short_t	unk = `0`
	)

inherited

the user can manually set the dirichlet Dofs for a given patch and unknown.

Parameters

[in]	vals	the values of the eliminated dofs.
[in]	unk	the considered unknown

◆ updateSolution()

template<class T >

void updateSolution	(	const gsMatrix< T > &	solVector,
		gsMultiPatch< T > &	result,
		T	theta = `(T)(1)`
	)		const

virtualinherited

Update solution by adding the computed solution vector to the current solution specified by.

result. This method assumes that all: unknowns have the same basis.

Parameters

[in]	solVector	the solution vector obtained from the linear system
[out]	result	the solution in form of a gsMultiBasis,

solVector is added to the: coefficients of result.

Parameters

[in] theta damping factor for update, theta = 1 corresponds to a full step.

Member Data Documentation

◆ m_bases

template<class T >

std::vector< gsMultiBasis<T> > m_bases

private

The discretization bases corresponding to the patches and to the number of solution fields that are to be computed. One entry of the vector corresponds to one or more unknown (in case of a system of PDEs)

◆ m_ddof

template<class T >

std::vector<gsMatrix<T> > m_ddof

private

Fixed DoF values (if applicable, for instance eliminated Dirichlet DoFs). One entry of the vector corresponds to exactly one unknown, i.e. the size of m_ddof must fit m_system.colBlocks().

◆ m_pde_ptr

template<class T >

memory::shared_ptr<gsPde<T> > m_pde_ptr

private

The PDE: contains multi-patch domain, boundary conditions and coeffcient functions

Detailed Description

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Constructor & Destructor Documentation

◆ gsGenericAssembler()

Member Function Documentation

◆ apply()

◆ assembleDG()

◆ assembleMass()

◆ assembleMass2()

◆ assembleMoments()

◆ assembleNeumann()

◆ assembleNitsche()

◆ assembleStiffness()

◆ computeDirichletDofs()

◆ computeDirichletDofsIntpl()

◆ computeDirichletDofsL2Proj()

◆ constructSolution() [1/2]

◆ constructSolution() [2/2]

◆ fixedDofs()

◆ fullMatrix() [1/2]

◆ fullMatrix() [2/2]

◆ homogenizeFixedDofs()

◆ penalizeDirichletDofs()

◆ setFixedDofs()

◆ setFixedDofVector()

◆ updateSolution()

Member Data Documentation

◆ m_bases

◆ m_ddof

◆ m_pde_ptr