gsExprAssembler< T > Class Template Reference

Detailed Description

template<class T>
class gismo::gsExprAssembler< T >

Assembler class for generating matrices and right-hand sides based on isogeometric expressions

Collaboration diagram for gsExprAssembler< T >:

Public Types
typedef gsExprHelper< T >::element	element
	Current element.
typedef gsExprHelper< T >::geometryMap	geometryMap
	Geometry map type.
typedef expr::gsFeSolution< T >	solution
	Solution type.
typedef gsExprHelper< T >::space	space
	Space type.
typedef gsExprHelper< T >::variable	variable
	Variable type.

Public Member Functions
template<class... expr>
void	assemble (const expr &... args)
	Adds the expressions args to the system matrix/rhs The arguments are considered as integrals over the whole domain.
template<class... expr>
void	assembleBdr (const bcRefList &BCs, expr &... args)
	Adds the expressions args to the system matrix/rhs The arguments are considered as integrals over the boundary parts in BCs.
template<class expr >
void	assembleJacobian (const expr residual, solution &u)
	Assembles the Jacobian matrix of the expression args with.
void	clearMatrix (const bool &save_sparsety_pattern=true)
	Re-Init Matrix (set zero by default)
template<class... expr>
void	computePattern (const expr &... args)
	Initializes the pattern of the sparse matrix.
template<class... expr>
void	computePatternBdr (const bcRefList &BCs, const expr &... args)
	Initializes the pattern of the sparse matrix at boundary integrals.
template<class... expr>
void	computePatternIfc (const ifContainer &iFaces, expr... args)
	Initializes the pattern of the sparse matrix at boundary integrals.
const FiberMatrix &	fiberMatrix () const
	Returns the internally stored sparse fiber matrix.
variable	getCoeff (const gsFunctionSet< T > &func)
expr::gsComposition< T >	getCoeff (const gsFunctionSet< T > &func, geometryMap &G)
geometryMap	getMap (const gsFunctionSet< T > &g)
	Registers g as an isogeometric geometry map and return a handle to it.
solution	getSolution (const expr::gsFeSpace< T > &s, gsMatrix< T > &cf) const
	Registers a representation of a solution variable from space s, based on the vector cf.
space	getSpace (const gsFunctionSet< T > &mp, index_t dim=1, index_t id=0)
space	getTestSpace (const gsFunctionSet< T > &mp, index_t dim=1, index_t id=0)
space	getTestSpace (space u, const gsFunctionSet< T > &mp, index_t dim=-1)
	Registers mp as an isogeometric test space corresponding to trial space u and return a handle to it.
	gsExprAssembler (index_t _rBlocks=1, index_t _cBlocks=1)
void	initMatrix ()
	Initializes the sparse matrix only.
void	initSystem (const index_t numRhs=1)
	Initializes the sparse system (sparse matrix and rhs)
void	initVector (const index_t numRhs=1)
	Initializes the right-hand side vector only.
const gsMultiBasis< T > &	integrationElements () const
	Returns the domain of integration.
const gsSparseMatrix< T > &	makeMatrix () const
const gsSparseMatrix< T > &	matrix () const
	Returns the left-hand global matrix.
void	matrix_into (gsSparseMatrix< T > &out)
	Writes the resulting matrix in out. The internal matrix is moved.
matBlockView	matrixBlockView ()
matConstBlockView	matrixBlockView () const
index_t	numBlocks () const
index_t	numDofs () const
	Returns the number of degrees of freedom (after initialization)
index_t	numTestDofs () const
	Returns the number of test functions (after initialization)
gsOptionList &	options ()
	Returns a reference to the options structure.
const gsMatrix< T > &	rhs () const
	Returns the right-hand side vector(s)
void	rhs_into (gsMatrix< T > &out)
	Writes the resulting vector in out. The internal data is moved.
void	setGeometryMap (const gsMultiPatch< T > &gMap)
	Set the geometrymap ( used for interface assembly)
void	setIntegrationElements (const gsMultiBasis< T > &mesh)
	Sets the domain of integration.
void	setOptions (gsOptionList opt)
	Set the assembler options.
space	testSpace (const index_t id)
space	testSpace (space u) const
	Return the test space of a pre-existing trial space u.
space	trialSpace (const index_t id) const
space	trialSpace (space &v) const
	Return the trial space of a pre-existing test space v.

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

Private Member Functions
void	resetDimensions ()
	Reset the dimensions of all involved spaces. Called internally by the init* functions.

Constructor & Destructor Documentation

gsExprAssembler()

template<class T >

gsExprAssembler ( index_t  _rBlocks = 1, index_t  _cBlocks = 1  )

inline

Constructor

Parameters

_rBlocks	Number of spaces for test functions
_cBlocks	Number of spaces for solution variables

Member Function Documentation

assemble()

template<class T >

template<class... expr>

void assemble

(

const expr &...

args

)

Adds the expressions args to the system matrix/rhs The arguments are considered as integrals over the whole domain.

See also

gsExprAssembler::setIntegrationElements

clearMatrix()

template<class T >

void clearMatrix ( const bool &  save_sparsety_pattern = true )

inline

Re-Init Matrix (set zero by default)

Parameters

save_sparsety_pattern

only modify values but keep sparsety information by multiplying matrix by zero in-place

defaultOptions()

template<class T >

gsOptionList defaultOptions ( )

static

Returns the list of default options for assembly.

dirichlet treatment? elimination ????

getCoeff() [1/2]

template<class T >

variable getCoeff ( const gsFunctionSet< T > &  func )

inline

Registers func as a variable and returns a handle to it

getCoeff() [2/2]

template<class T >

expr::gsComposition< T > getCoeff ( const gsFunctionSet< T > &  func, geometryMap &  G  )

inline

Registers func as a variable defined on G and returns a handle to it

getSolution()

template<class T >

solution getSolution ( const expr::gsFeSpace< T > &  s, gsMatrix< T > &  cf  ) const

inline

Registers a representation of a solution variable from space s, based on the vector cf.

The vector cf should have the structure of the columns of the system matrix this->matrix(). The returned handle corresponds to a function in the space s

getSpace()

template<class T >

space getSpace ( const gsFunctionSet< T > &  mp, index_t  dim = 1, index_t  id = 0  )

inline

Registers mp as an isogeometric (both trial and test) space and returns a handle to it

getTestSpace() [1/2]

template<class T >

space getTestSpace ( const gsFunctionSet< T > &  mp, index_t  dim = 1, index_t  id = 0  )

inline

Registers mp as an isogeometric test (row) space and returns a handle to it

getTestSpace() [2/2]

template<class T >

space getTestSpace ( space  u, const gsFunctionSet< T > &  mp, index_t  dim = -1  )

inline

Registers mp as an isogeometric test space corresponding to trial space u and return a handle to it.

Note

Both test and trial spaces are registered at once by gsExprAssembler::getSpace.

Use this function after calling gsExprAssembler::getSpace when a distinct test space is requred (eg. Petrov-Galerkin methods).

Note

The dimension is set to the same as u, unless the caller sets as a third argument a new value.

makeMatrix()

template<class T >

const gsSparseMatrix< T > & makeMatrix ( ) const

inline

When calling assemble, the matrix is not filled (but only stored internally). Call this function to fill the sparsematrix with all the assemblies so far

matrixBlockView() [1/2]

template<class T >

matBlockView matrixBlockView ( )

inline

Returns a block view of the system matrix, each block corresponding to a different space, or to different groups of dofs, in case of calar problems

matrixBlockView() [2/2]

template<class T >

matConstBlockView matrixBlockView ( ) const

inline

Returns a const block view of the system matrix, each block corresponding to a different space, or to different groups of dofs, in case of calar problems

numBlocks()

template<class T >

index_t numBlocks ( ) const

inline

Returns the number of blocks in the matrix, corresponding to variables/components

setGeometryMap()

template<class T >

void setGeometryMap ( const gsMultiPatch< T > &  gMap )

inline

Set the geometrymap ( used for interface assembly)

Warning

Must be called before any computation is requested

setIntegrationElements()

template<class T >

void setIntegrationElements ( const gsMultiBasis< T > &  mesh )

inline

Sets the domain of integration.

Warning

Must be called before any computation is requested

testSpace()

template<class T >

space testSpace ( const index_t  id )

inline

Return the variable (previously created by getTrialSpace) with the given id

trialSpace()

template<class T >

space trialSpace ( const index_t  id ) const

inline

Return a variable handle (previously created by getSpace) for unknown id