21 template <
bool strainlike = false>
22 struct VectorizeWithVoigt
24 template <
typename Derived>
25 static auto transform(
const Eigen::DenseBase<Derived>& mat) {
26 return toVoigt(mat.derived(), strainlike);
30 struct VectorizeGeneric
32 template <
typename Derived>
33 static auto transform(
const Eigen::DenseBase<Derived>& mat) {
34 return mat.derived().reshaped().eval();
38 template <
bool strainlike = false>
39 struct MatricizeWithVoigt
41 template <
typename Derived,
int RowsAtCompileTime,
int ColsAtCompileTime>
42 static auto transform(
const Eigen::DenseBase<Derived>& vec,
int rows = RowsAtCompileTime,
43 int cols = ColsAtCompileTime) {
44 assert(rows == RowsAtCompileTime && cols == ColsAtCompileTime &&
45 "Only the fixed size values work for voigt matrices and vectors");
46 static_assert(RowsAtCompileTime != Eigen::Dynamic and ColsAtCompileTime != Eigen::Dynamic,
47 "Voigt notation only available for fixed size vectors and matrices");
48 return fromVoigt(vec.derived(), strainlike);
52 struct MatricizeGeneric
54 template <
typename Derived,
int RowsAtCompileTime,
int ColsAtCompileTime>
55 static auto transform(
const Eigen::DenseBase<Derived>& vec,
int rows = RowsAtCompileTime,
56 int cols = ColsAtCompileTime) {
57 return vec.derived().reshaped(Eigen::fix<RowsAtCompileTime>(rows), Eigen::fix<ColsAtCompileTime>(cols)).eval();
63namespace ResultTypes {
64#define REGISTER_RESULTTYPE_IMPL(resultTypeName, rowsExpr, colsExpr, MaxRowsExpr, MaxColsExpr, VectorizeStruct, \
66 template <typename ScalarType, int gridDim, int worldDim> \
67 struct resultTypeName \
69 friend std::string toString(resultTypeName) { return #resultTypeName; } \
71 using type = Eigen::Matrix<ScalarType, rowsExpr, colsExpr, 0, MaxRowsExpr, MaxColsExpr>; \
72 using Vectorizer = VectorizeStruct; \
73 using Matricizer = MatricizeStruct; \
75 template <typename ScalarType_, int gridDim_, int worldDim_> \
76 using Rebind = resultTypeName<ScalarType_, gridDim_, worldDim_>; \
86#define REGISTER_SIMPLE_SYMMETRIC_RESULTTYPE(resultTypeName, rowsExpr, colsExpr, strainlike) \
87 REGISTER_RESULTTYPE_IMPL(resultTypeName, rowsExpr, colsExpr, rowsExpr, colsExpr, \
88 Ikarus::Impl::VectorizeWithVoigt<strainlike>, Ikarus::Impl::MatricizeWithVoigt<strainlike>)
97#define REGISTER_RESULTTYPE(resultTypeName, rowsExpr, colsExpr) \
98 REGISTER_RESULTTYPE_IMPL(resultTypeName, rowsExpr, colsExpr, Ikarus::Impl::VectorizeGeneric, \
99 Ikarus::Impl::MatricizeGeneric)
109#define REGISTER_RESERVED_RESULTTYPE(resultTypeName, rowsExpr, colsExpr, MaxRowsExpr, MaxColsExpr) \
110 REGISTER_RESULTTYPE_IMPL(resultTypeName, rowsExpr, colsExpr, MaxRowsExpr, MaxColsExpr, \
111 Ikarus::Impl::VectorizeGeneric, Ikarus::Impl::MatricizeGeneric)
119#define REGISTER_SIMPLE_RESULTTYPE(resultTypeName, rowsExpr, colsExpr) \
120 REGISTER_RESERVED_RESULTTYPE(resultTypeName, rowsExpr, colsExpr, rowsExpr, colsExpr)
153template <
typename RT, ResultShape storedResultShape = ResultShape::Vector>
157 using ResultTypeValueType =
typename RT::type;
158 static constexpr Eigen::Index rowsAtCompileTime = ResultTypeValueType::RowsAtCompileTime;
159 static constexpr Eigen::Index colsAtCompileTime = ResultTypeValueType::ColsAtCompileTime;
163 using VecType = std::invoke_result_t<
decltype(&RT::Vectorizer::template transform<ResultTypeValueType>),
164 const ResultTypeValueType&>;
166 std::invoke_result_t<
decltype(&RT::Matricizer::template transform<VecType, rowsAtCompileTime, colsAtCompileTime>),
168 using StoredType = std::conditional_t<storedValueIsVector, VecType, MatType>;
176 if constexpr (storedValueIsVector)
179 return RT::Vectorizer::transform(value_);
188 auto asMat(Eigen::Index rows = rowsAtCompileTime, Eigen::Index cols = colsAtCompileTime)
const {
189 if constexpr (storedValueIsVector) {
190 if constexpr (rowsAtCompileTime == Eigen::Dynamic)
191 assert(rows != rowsAtCompileTime &&
192 "For dynamic size result types you have to pass rows by hand, since it is not clear how the result "
193 "should be reshaped");
194 if constexpr (colsAtCompileTime == Eigen::Dynamic)
195 assert(cols != colsAtCompileTime &&
196 "For dynamic size result types you have to pass cols by hand, since it is not clear how the result "
197 "should be reshaped");
198 return RT::Matricizer::template transform<VecType, rowsAtCompileTime, colsAtCompileTime>(value_, rows, cols);
206 this->value_ = value;
210 this->value_ = std::move(value);
219 template <
template <
typename,
int,
int>
class RT>
220 using DummyRT = RT<double, 1, 1>;
228template <
template <
typename,
int,
int>
class RT>
230 return Impl::DummyRT<RT>{};
238template <
template <
typename,
int,
int>
class RT>
240 return toString(Impl::DummyRT<RT>{});
248template <
template <
typename,
int,
int>
class RT1,
template <
typename,
int,
int>
class RT2>
Helper for the Eigen::Tensor types.
auto fromVoigt(const Eigen::Vector< ST, size > &EVoigt, bool isStrain=true)
Converts a vector given in Voigt notation to a matrix.
Definition: tensorutils.hh:256
constexpr Eigen::Index toVoigt(Eigen::Index i, Eigen::Index j) noexcept
Converts 2D indices to Voigt notation index.
Definition: tensorutils.hh:166
Definition: simpleassemblers.hh:22
constexpr std::string toString(ScalarAffordances _e)
Definition: ferequirements.hh:37
auto makeRT()
Creates a dummy resultType which can be stored in a variable.
Definition: feresulttypes.hh:229
constexpr bool isSameResultType
Meta variable to test whether two ResultType templates are the same.
Definition: feresulttypes.hh:249
ResultShape
Definition: feresulttypes.hh:142
REGISTER_SIMPLE_SYMMETRIC_RESULTTYPE(linearStress, worldDim, worldDim, false)
REGISTER_SIMPLE_RESULTTYPE(director, worldDim, 1)
Container that is used for FE Results. It gives access to the stored value, but can also be used to a...
Definition: feresulttypes.hh:155
std::conditional_t< storedValueIsVector, VecType, MatType > StoredType
Definition: feresulttypes.hh:168
std::invoke_result_t< decltype(&RT::Matricizer::template transform< VecType, rowsAtCompileTime, colsAtCompileTime >), const VecType &, int, int > MatType
Definition: feresulttypes.hh:167
ResultWrapper(StoredType &&value)
Definition: feresulttypes.hh:203
ResultWrapper & operator=(StoredType &&value)
Definition: feresulttypes.hh:209
std::invoke_result_t< decltype(&RT::Vectorizer::template transform< ResultTypeValueType >), const ResultTypeValueType & > VecType
Definition: feresulttypes.hh:164
ResultWrapper & operator=(const StoredType &value)
Definition: feresulttypes.hh:205
auto asVec() const
Returns the stored value as Vector.
Definition: feresulttypes.hh:175
ResultWrapper(const StoredType &value)
Definition: feresulttypes.hh:204
auto asMat(Eigen::Index rows=rowsAtCompileTime, Eigen::Index cols=colsAtCompileTime) const
Returns the stored value as Matrix (if possible)
Definition: feresulttypes.hh:188
RT ResultType
Definition: feresulttypes.hh:169