crp/OpenMC
1
0
Fork 0
mirror of https://github.com/openmc-dev/openmc.git synced 2026-07-28 22:26:08 -04:00
Code Issues Projects Releases Packages Wiki Activity Actions 5

for multigroup fission matrix does not have same nonzero pattern as destruction matrix, fixed PETSC flag

Browse source
This commit is contained in:
Bryan Herman 2012-04-05 06:28:18 -07:00
parent 7e8a023381
commit bf9ae7045a
1 changed files with 1 additions and 1 deletions
Download patch file Download diff file Expand all files Collapse all files

2
src/cmfd_jacobian_operator.F90
View file

@ -223,7 +223,7 @@ contains
call MPI_BCAST(lambda,1,MPI_REAL8,n_procs_cmfd-1,PETSC_COMM_WORLD,ierr)
! compute math (M-lambda*F) M is overwritten here
call MatAXPY(ctx%loss%M,-1.0_8*lambda,ctx%prod%F,SUBSET_NONZERO_PATTERN,ierr)
call MatAXPY(ctx%loss%M,-1.0_8*lambda,ctx%prod%F,DIFFERENT_NONZERO_PATTERN,ierr)
! create tmp petsc vector for source
call VecCreateMPI(PETSC_COMM_WORLD,ctx%loss%localn,PETSC_DECIDE,source,ierr)