Don't allocate memory for eigendecomposition under COMPUTATION_ACTION

beagle-dev · Feb 13, 2024 · 33c0534 · 33c0534
1 parent 1a5d3ba
commit 33c0534
Showing 1 changed file with 3 additions and 1 deletion.
diff --git a/libhmsbeagle/CPU/BeagleCPUImpl.hpp b/libhmsbeagle/CPU/BeagleCPUImpl.hpp
@@ -335,7 +335,9 @@ int BeagleCPUImpl<BEAGLE_CPU_GENERIC>::createInstance(int tipCount,
     else
         kFlags |= BEAGLE_FLAG_THREADING_NONE;
 
-    if (kFlags & BEAGLE_FLAG_EIGEN_COMPLEX)
+    if (preferenceFlags & BEAGLE_FLAG_COMPUTATION_ACTION)
+	gEigenDecomposition = nullptr;
+    else if (kFlags & BEAGLE_FLAG_EIGEN_COMPLEX)
         gEigenDecomposition = new EigenDecompositionSquare<BEAGLE_CPU_EIGEN_GENERIC>(kEigenDecompCount,
                 kStateCount,kCategoryCount,kFlags);
     else