change solve to splu then solve

openmc-dev · Nov 28, 2023 · 7dbd5e5 · 7dbd5e5
1 parent 8b2698f
commit 7dbd5e5
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/openmc/deplete/cram.py b/openmc/deplete/cram.py
@@ -79,7 +79,7 @@ def __call__(self, A, n0, dt):
         y = n0.copy()
         ident = sp.eye(A.shape[0], format='csc')
         for alpha, theta in zip(self.alpha, self.theta):
-            y += 2*np.real(alpha*sla.spsolve(A - theta*ident, y))
+            y += 2*np.real(alpha*sla.splu(A - theta*ident).solve(y))
         return y * self.alpha0