split-up rotate tensor function

ParameterLib/CoordinateSystem.cpp

@@ -293,8 +293,7 @@ Eigen::Matrix<double, Dimension, Dimension> CoordinateSystem::rotateTensor(
     auto const tensor =
         Eigen::Map<Eigen::Matrix<double, Dimension, Dimension> const>(
             values.data(), Dimension, Dimension);
-    auto const R = transformation<Dimension>(pos);
-    return R * tensor * R.transpose();
+    return rotateTensor<Dimension>(tensor, pos);
 }
 
 template <int Dimension>

ParameterLib/CoordinateSystem.h

@@ -72,6 +72,11 @@ struct CoordinateSystem final
     Eigen::Matrix<double, Dimension, Dimension> rotateTensor(
         std::vector<double> const& values, SpatialPosition const& pos) const;
 
+    template <int Dimension, typename Derived>
+    Eigen::Matrix<double, Dimension, Dimension> rotateTensor(
+        Eigen::MatrixBase<Derived> const& tensor,
+        SpatialPosition const& pos) const;
+
     template <int Dimension>
     Eigen::Matrix<double, Dimension, Dimension> rotateDiagonalTensor(
         std::vector<double> const& values, SpatialPosition const& pos) const;
@@ -84,6 +89,14 @@ private:
         SpatialPosition const& pos) const;
 };
 
+template <int Dimension, typename Derived>
+Eigen::Matrix<double, Dimension, Dimension> CoordinateSystem::rotateTensor(
+    Eigen::MatrixBase<Derived> const& tensor, SpatialPosition const& pos) const
+{
+    auto const R = transformation<Dimension>(pos);
+    return R * tensor * R.transpose();
+}
+
 extern template Eigen::Matrix<double, 2, 2> CoordinateSystem::rotateTensor<2>(
     std::vector<double> const& values, SpatialPosition const& pos) const;
 extern template Eigen::Matrix<double, 3, 3> CoordinateSystem::rotateTensor<3>(