From 9442aa0fe03afd77dd26a74a7c2aff12421133e0 Mon Sep 17 00:00:00 2001 From: Bryan Herman Date: Wed, 4 Apr 2012 07:39:52 -0700 Subject: [PATCH] fixed multigroup index mapping --- src/cmfd_loss_operator.F90 | 6 +++--- src/cmfd_prod_operator.F90 | 6 +++--- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/src/cmfd_loss_operator.F90 b/src/cmfd_loss_operator.F90 index 28633e7c6f..0388906a5c 100644 --- a/src/cmfd_loss_operator.F90 +++ b/src/cmfd_loss_operator.F90 @@ -460,9 +460,9 @@ contains ! get indices from indexmap g = mod(irow,ng) + 1 - i = cmfd % indexmap(irow+1,1) - j = cmfd % indexmap(irow+1,2) - k = cmfd % indexmap(irow+1,3) + i = cmfd % indexmap(irow/ng+1,1) + j = cmfd % indexmap(irow/ng+1,2) + k = cmfd % indexmap(irow/ng+1,3) else diff --git a/src/cmfd_prod_operator.F90 b/src/cmfd_prod_operator.F90 index 31450370e8..2b43c2ea83 100644 --- a/src/cmfd_prod_operator.F90 +++ b/src/cmfd_prod_operator.F90 @@ -281,9 +281,9 @@ contains ! get indices from indexmap g = mod(irow,ng) + 1 - i = cmfd % indexmap(irow+1,1) - j = cmfd % indexmap(irow+1,2) - k = cmfd % indexmap(irow+1,3) + i = cmfd % indexmap(irow/ng+1,1) + j = cmfd % indexmap(irow/ng+1,2) + k = cmfd % indexmap(irow/ng+1,3) else