[PL] Fixed IP data getter

ProcessLib/Utils/SetOrGetIntegrationPointData.h

@@ -16,7 +16,7 @@
 
 #include "BaseLib/DynamicSpan.h"
 #include "MathLib/KelvinVector.h"
-#include "MathLib/LinAlg/Eigen/EigenMapTools.h"
+#include "TransposeInPlace.h"
 
 namespace ProcessLib
 {
@@ -65,17 +65,25 @@ std::vector<double> const& getIntegrationPointKelvinVectorData(
     return cache;
 }
 
+//! Overload without \c cache argument.
+//!
+//! \note This function returns the data in transposed storage order compared to
+//! the overloads that have a \c cache argument.
 template <int DisplacementDim, typename IntegrationPointDataVector,
           typename MemberType>
 std::vector<double> getIntegrationPointKelvinVectorData(
     IntegrationPointDataVector const& ip_data_vector, MemberType member)
 {
-    std::vector<double> ip_kelvin_vector_values;
-
-    getIntegrationPointKelvinVectorData<DisplacementDim>(
-        ip_data_vector, member, ip_kelvin_vector_values);
+    constexpr int kelvin_vector_size =
+        MathLib::KelvinVector::kelvin_vector_dimensions(DisplacementDim);
 
-    return ip_kelvin_vector_values;
+    return transposeInPlace<kelvin_vector_size>(
+        [&](std::vector<double>& values)
+        {
+            return getIntegrationPointKelvinVectorData<DisplacementDim>(
+                ip_data_vector, member, values);
+            ;
+        });
 }
 
 template <int DisplacementDim, typename IntegrationPointDataVector,