From d4b703eb908519ab639f572cd28ecbf2fd2d15cf Mon Sep 17 00:00:00 2001
From: Dmitri Naumov <dmitri.naumov@ufz.de>
Date: Wed, 12 Dec 2018 12:04:26 +0100
Subject: [PATCH] [MaL] Specialize for dynamic Kelvin to tensor conv
---
MathLib/KelvinVector.cpp | 26 ++++++++++++++++++++++++++
1 file changed, 26 insertions(+)
diff --git a/MathLib/KelvinVector.cpp b/MathLib/KelvinVector.cpp
index 0c6718ddb8f..1359f452a8d 100644
--- a/MathLib/KelvinVector.cpp
+++ b/MathLib/KelvinVector.cpp
@@ -80,6 +80,32 @@ Eigen::Matrix<double, 3, 3> kelvinVectorToTensor(
return m;
}
+template <>
+Eigen::Matrix<double, 3, 3> kelvinVectorToTensor(Eigen::Matrix<double,
+ Eigen::Dynamic,
+ 1,
+ Eigen::ColMajor,
+ Eigen::Dynamic,
+ 1> const& v)
+{
+ if (v.size() == 4)
+ {
+ Eigen::Matrix<double, 4, 1, Eigen::ColMajor, 4, 1> v4;
+ v4 << v[0], v[1], v[2], v[3];
+ return kelvinVectorToTensor(v4);
+ }
+ if (v.size() == 6)
+ {
+ Eigen::Matrix<double, 6, 1, Eigen::ColMajor, 6, 1> v6;
+ v6 << v[0], v[1], v[2], v[3], v[4], v[5];
+ return kelvinVectorToTensor(v6);
+ }
+ OGS_FATAL(
+ "Conversion of dynamic Kelvin vector of size %d to a tensor is not "
+ "possible. Kelvin vector must be of size 4 or 6.",
+ v.size());
+}
+
template <>
KelvinVectorType<2> tensorToKelvin<2>(Eigen::Matrix<double, 3, 3> const& m)
{
--
GitLab