COMP: Avoid type conversion errors on Mac with new ITK

dzenanz · Jan 17, 2024 · e18b1d2 · e18b1d2
1 parent 62550fc
commit e18b1d2
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/Tensor/TensorFunctions.h b/Tensor/TensorFunctions.h
@@ -438,7 +438,7 @@ GetMetricTensorCost(TVectorType dpath, TTensorType dtv, unsigned int matrixpower
   vec(0, 0) = dpath[0];
   vec(1, 0) = dpath[1];
   vec(2, 0) = dpath[2];
-  MatrixType inv = vnl_matrix_inverse<double>(DT);
+  MatrixType inv = vnl_matrix_inverse<double>(DT).inverse();
   for (unsigned int lo = 1; lo < matrixpower; lo++)
   {
     inv = inv * inv;
@@ -837,7 +837,7 @@ GetMetricTensorCost(itk::Vector<float, 3> dpath, TTensorType dtv)
   vec(0, 0) = dpath[0];
   vec(1, 0) = dpath[1];
   vec(2, 0) = dpath[2];
-  MatrixType inv = vnl_matrix_inverse<double>(DT);
+  MatrixType inv = vnl_matrix_inverse<double>(DT).inverse();
 
   MatrixType sol = vec.transpose() * inv * vec;
   float      cost = sol(0, 0) / etot;