Element.h
No OneTemporary
Actions

Subscribers

None

File Metadata

	/*

	(c - GPLv3) T.W.J. de Geus (Tom) \| tom@geus.me \| www.geus.me \| github.com/tdegeus/GooseFEM

	*/

	#ifndef GOOSEFEM_ELEMENT_H
	#define GOOSEFEM_ELEMENT_H

	#include "config.h"

	namespace GooseFEM {
	namespace Element {

	// Convert nodal vector [nnode, ndim] to nodal vector stored per element [nelem, nne, ndim]
	inline xt::xtensor<double, 3> asElementVector(
	const xt::xtensor<size_t, 2>& conn, const xt::xtensor<double, 2>& nodevec);

	// Assemble nodal vector stored per element [nelem, nne, ndim] to nodal vector [nnode, ndim]
	inline xt::xtensor<double, 2> assembleNodeVector(
	const xt::xtensor<size_t, 2>& conn, const xt::xtensor<double, 3>& elemvec);

	// Check that DOFs leave no holes
	template <class E>
	inline bool isSequential(const E& dofs);

	// Check structure of the matrices stored per element [nelem, nnendim, nnendim]
	bool isDiagonal(const xt::xtensor<double, 3>& elemmat);

	} // namespace Element
	} // namespace GooseFEM

	#include "Element.hpp"

	#endif