Tensor fix

Another Kelvin issue

MaterialLib/MPL/Properties/EmbeddedFracturePermeability.cpp

@@ -58,7 +58,8 @@ PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::value(
         {
             return _n;
         }
-        auto const sigma = formEigenTensor<3>(std::get<SymmetricTensor>(
+        auto const sigma = MathLib::KelvinVector::kelvinVectorToTensor(
+        std::get<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
             variable_array[static_cast<int>(Variable::total_stress)]));
         Eigen::SelfAdjointEigenSolver<Eigen::Matrix<double, 3, 3>> e_s(sigma);
         return (Eigen::Matrix<double, 3, 1>)e_s.eigenvectors().col(2);
@@ -97,7 +98,8 @@ PropertyDataType EmbeddedFracturePermeability<DisplacementDim>::dValue(
         {
             return _n;
         }
-        auto const sigma = formEigenTensor<3>(std::get<SymmetricTensor>(
+        auto const sigma = MathLib::KelvinVector::kelvinVectorToTensor(
+        std::get<MathLib::KelvinVector::KelvinVectorType<DisplacementDim>>(
             variable_array[static_cast<int>(Variable::total_stress)]));
         Eigen::SelfAdjointEigenSolver<Eigen::Matrix<double, 3, 3>> e_s(sigma);
         return (Eigen::Matrix<double, 3, 1>)e_s.eigenvectors().col(2);