crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-27 21:55:41 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

Make sure IPFCramSolver.__call__ doesn't modify n0

Browse source
This commit is contained in:
Paul Romano 2020-06-05 21:18:46 -05:00
parent 506d3e610d
commit 36e316bb2c
1 changed files with 1 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
openmc/deplete/cram.py
View file

@ -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))