diff options
Diffstat (limited to 'Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90')
| -rw-r--r-- | Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90 | 14 |
1 files changed, 7 insertions, 7 deletions
diff --git a/Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90 b/Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90 index 94a9e60f1..bd3f9b396 100644 --- a/Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90 +++ b/Carpet/CarpetLib/src/prolongate_3d_real8_eno.F90 @@ -194,17 +194,17 @@ subroutine prolongate_3d_real8_eno ( & !!$ Loop over fine region + !$omp parallel do collapse(3) private(i,j,k, i0,fi,j0,fj,k0,fk, tmp1,tmp2, ii,jj,kk) do k = 0, regkext-1 - k0 = (srckoff + k) / dstkfac - fk = mod(srckoff + k, dstkfac) - do j = 0, regjext-1 - j0 = (srcjoff + j) / dstjfac - fj = mod(srcjoff + j, dstjfac) - do i = 0, regiext-1 + i0 = (srcioff + i) / dstifac fi = mod(srcioff + i, dstifac) + j0 = (srcjoff + j) / dstjfac + fj = mod(srcjoff + j, dstjfac) + k0 = (srckoff + k) / dstkfac + fk = mod(srckoff + k, dstkfac) !!$ Where is the fine grid point w.r.t the coarse grid? @@ -279,7 +279,7 @@ subroutine prolongate_3d_real8_eno ( & eno1d(tmp2(0:3)) case default - call CCTK_WARN(0, "Internal error in ENO prolongation. Should only be used with refinement factor 2!") + call CCTK_ERROR("Internal error in ENO prolongation. Should only be used with refinement factor 2!") end select end do |