From 9830efaf2aadf68ffde9a054388bacef20ab5645 Mon Sep 17 00:00:00 2001 From: Ethan Peterson Date: Tue, 21 Nov 2023 17:22:42 -0500 Subject: [PATCH] change identity matrix format to csr in cram (#2771) --- openmc/deplete/cram.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/openmc/deplete/cram.py b/openmc/deplete/cram.py index f5f80dfdc4..8e7bb819dc 100644 --- a/openmc/deplete/cram.py +++ b/openmc/deplete/cram.py @@ -77,7 +77,7 @@ class IPFCramSolver(DepSystemSolver): """ A = sp.csr_matrix(A * dt, dtype=np.float64) y = n0.copy() - ident = sp.eye(A.shape[0]) + ident = sp.eye(A.shape[0], format='csr') for alpha, theta in zip(self.alpha, self.theta): y += 2*np.real(alpha*sla.spsolve(A - theta*ident, y)) return y * self.alpha0