Normalize determinants of random scaling matrices

polatory · Aug 13, 2024 · 27e01a9 · 27e01a9
1 parent dff1568
commit 27e01a9
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/test/utility.hpp b/test/utility.hpp
@@ -30,6 +30,8 @@ polatory::geometry::matrixNd<Dim> random_scaling() {
 
   Matrix scale = Matrix::Identity();
   scale.diagonal().array() *= pow(10.0, 0.5 * Vector::Random().array());
+  // Normalize the determinant.
+  scale.diagonal() *= std::pow(scale.diagonal().prod(), -1.0 / Dim);
 
   return scale;
 }