Class GooseFEM::Vector

class Vector

Class to switch between storage types.

In particular:

• “dofval”: DOF values [ndof].

• ”nodevec”: nodal vectors [nnode, ndim].

• ”elemvec”: nodal vectors stored per element [nelem, nne, ndim].

Subclassed by GooseFEM::VectorPartitioned, GooseFEM::VectorPartitionedTyings

Public Functions

template<class S, class T>
inline Vector(const S &conn, const T &dofs)

Constructor.

Parameters:

• conn – connectivity [nelem, nne].

• dofs – DOFs per node [nnode, ndim].

inline size_t nelem() const

Returns:

Number of elements.

inline size_t nne() const

Returns:

Number of nodes per element.

inline size_t nnode() const

Returns:

Number of nodes.

inline size_t ndim() const

Returns:

Number of dimensions.

inline size_t ndof() const

Returns:

Number of DOFs.

inline const array_type::tensor<size_t, 2> &conn() const

Returns:

Connectivity (nodes per element) [nelem, nne].

inline const array_type::tensor<size_t, 2> &dofs() const

Returns:

DOFs per node [nnode, ndim]

template<class T>
inline T Copy(const T &nodevec_src, const T &nodevec_dest) const

Copy “nodevec” to another “nodevec”.

Parameters:

• nodevec_src – input [nnode, ndim]

• nodevec_dest – input [nnode, ndim]

Returns:

nodevec output [nnode, ndim]

template<class T>
inline void copy(const T &nodevec_src, T &nodevec_dest) const

Copy “nodevec” to another “nodevec”.

Parameters:

• nodevec_src – input [nnode, ndim]

• nodevec_dest – output [nnode, ndim]

template<class T>
inline array_type::tensor<double, 1> AsDofs(const T &arg) const

Convert “nodevec” or “elemvec” to “dofval” (overwrite entries that occur more than once).

Parameters:

arg – nodevec [nnode, ndim] or elemvec [nelem, nne, ndim]

Returns:

dofval [ndof]

template<class T, class R>
inline void asDofs(const T &arg, R &ret) const

Convert “nodevec” or “elemvec” to “dofval” (overwrite entries that occur more than once).

Parameters:

• arg – nodevec [nnode, ndim] or elemvec [nelem, nne, ndim]

• ret – dofval (output) [ndof]

template<class T>
inline array_type::tensor<double, 2> AsNode(const T &arg) const

Convert “dofval” or “elemvec” to “nodevec” (overwrite entries that occur more than once).

Parameters:

arg – dofval [ndof] or elemvec [nelem, nne, ndim]

Returns:

nodevec output [nnode, ndim]

template<class T, class R>
inline void asNode(const T &arg, R &ret) const

Convert “dofval” or “elemvec” to “nodevec” (overwrite entries that occur more than once).

Parameters:

• arg – dofval [ndof] or elemvec [nelem, nne, ndim]

• ret – nodevec, output [nnode, ndim]

template<class T>
inline array_type::tensor<double, 3> AsElement(const T &arg) const

Convert “dofval” or “nodevec” to “elemvec” (overwrite entries that occur more than once).

Parameters:

arg – dofval [ndof] or nodevec [nnode, ndim].

Returns:

elemvec output [nelem, nne, ndim].

template<class T, class R>
inline void asElement(const T &arg, R &ret) const

Convert “dofval” or “nodevec” to “elemvec” (overwrite entries that occur more than once).

Parameters:

• arg – dofval [ndof] or nodevec [nnode, ndim].

• ret – elemvec, output [nelem, nne, ndim].

template<class T>
inline array_type::tensor<double, 1> AssembleDofs(const T &arg) const

Assemble “nodevec” or “elemvec” to “dofval” (adds entries that occur more that once).

Parameters:

arg – nodevec [nnode, ndim] or elemvec [nelem, nne, ndim]

Returns:

dofval output [ndof]

template<class T, class R>
inline void assembleDofs(const T &arg, R &ret) const

Assemble “nodevec” or “elemvec” to “dofval” (adds entries that occur more that once).

Parameters:

• arg – nodevec [nnode, ndim] or elemvec [nelem, nne, ndim]

• ret – dofval, output [ndof]

template<class T>
inline array_type::tensor<double, 2> AssembleNode(const T &arg) const

Assemble “elemvec” to “nodevec” (adds entries that occur more that once.

Parameters:

arg – elemvec [nelem, nne, ndim]

Returns:

nodevec output [nnode, ndim]

template<class T, class R>
inline void assembleNode(const T &arg, R &ret) const

Assemble “elemvec” to “nodevec” (adds entries that occur more that once.

Parameters:

• arg – elemvec [nelem, nne, ndim]

• ret – nodevec, output [nnode, ndim]

inline std::array<size_t, 1> shape_dofval() const

Shape of “dofval”.

Returns:

[ndof]

inline std::array<size_t, 2> shape_nodevec() const

Shape of “nodevec”.

Returns:

[nnode, ndim]

inline std::array<size_t, 3> shape_elemvec() const

Shape of “elemvec”.

Returns:

[nelem, nne, ndim]

inline std::array<size_t, 3> shape_elemmat() const

Shape of “elemmat”.

Returns:

[nelem, nne * ndim, nne * ndim]

inline array_type::tensor<double, 1> allocate_dofval() const

Allocated “dofval”.

Returns:

[ndof]

inline array_type::tensor<double, 1> allocate_dofval(double val) const

Allocated and initialised “dofval”.

Parameters:

val – value to which to initialise.

Returns:

[ndof]

inline array_type::tensor<double, 2> allocate_nodevec() const

Allocated “nodevec”.

Returns:

[nnode, ndim]

inline array_type::tensor<double, 2> allocate_nodevec(double val) const

Allocated and initialised “nodevec”.

Parameters:

val – value to which to initialise.

Returns:

[nnode, ndim]

inline array_type::tensor<double, 3> allocate_elemvec() const

Allocated “elemvec”.

Returns:

[nelem, nne, ndim]

inline array_type::tensor<double, 3> allocate_elemvec(double val) const

Allocated and initialised “elemvec”.

Parameters:

val – value to which to initialise.

Returns:

[nelem, nne, ndim]

inline array_type::tensor<double, 3> allocate_elemmat() const

Allocated “elemmat”.

Returns:

[nelem, nne * ndim, nne * ndim]

inline array_type::tensor<double, 3> allocate_elemmat(double val) const

Allocated and initialised “elemmat”.

Parameters:

val – value to which to initialise.

Returns:

[nelem, nne * ndim, nne * ndim]