From 36e316bb2ca6fdccad882b590d91eecaee3f3084 Mon Sep 17 00:00:00 2001
From: Paul Romano <paul.k.romano@gmail.com>
Date: Fri, 5 Jun 2020 21:18:46 -0500
Subject: [PATCH] Make sure IPFCramSolver.__call__ doesn't modify n0

---
 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 6b6a909ec1..f5f80dfdc4 100644
--- a/openmc/deplete/cram.py
+++ b/openmc/deplete/cram.py
@@ -76,7 +76,7 @@ class IPFCramSolver(DepSystemSolver):
 
         """
         A = sp.csr_matrix(A * dt, dtype=np.float64)
-        y = np.asarray(n0, dtype=np.float64)
+        y = n0.copy()
         ident = sp.eye(A.shape[0])
         for alpha, theta in zip(self.alpha, self.theta):
             y += 2*np.real(alpha*sla.spsolve(A - theta*ident, y))