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

change identity matrix format to csr in cram (#2771)

Browse source
This commit is contained in:
Ethan Peterson 2023-11-21 17:22:42 -05:00 • committed by GitHub
parent 24e1c95161
commit 9830efaf2a
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
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

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