diff options
| author | schnetter <> | 2003-06-24 12:00:00 +0000 |
|---|---|---|
| committer | schnetter <> | 2003-06-24 12:00:00 +0000 |
| commit | bbfc3b0b532c073ac8aeb498db3e0fbd7ca81b41 (patch) | |
| tree | 5bc3eb38f882f8455ac79646597b33d0e85b657e /Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 | |
| parent | f4483891ba4b6dfc56ea0751011278fad85081ae (diff) | |
Add 5th order spatial interpolation.
darcs-hash:20030624120058-07bb3-01d8634b9ffbc710b6102af4d3128f32b73f21d2.gz
Diffstat (limited to 'Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77')
| -rw-r--r-- | Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 | 9 |
1 files changed, 2 insertions, 7 deletions
diff --git a/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 b/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 index 28f8c155f..e664ac47e 100644 --- a/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 +++ b/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77 @@ -1,5 +1,5 @@ c -*-Fortran-*- -c $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77,v 1.3 2004/03/11 12:03:09 schnetter Exp $ +c $Header: /home/eschnett/C/carpet/Carpet/Carpet/CarpetLib/src/prolongate_3d_real8_3tl_o5.F77,v 1.1 2003/06/24 14:00:58 schnetter Exp $ #include "cctk.h" @@ -90,11 +90,6 @@ c bbox(:,3) is stride c This could be handled, but is likely to point to an error elsewhere call CCTK_WARN (0, "Internal error: region extent is empty") end if - if (mod(srcbbox(d,2) - srcbbox(d,1), srcbbox(d,3)).ne.0 - $ .or. mod(dstbbox(d,2) - dstbbox(d,1), dstbbox(d,3)).ne.0 - $ .or. mod(regbbox(d,2) - regbbox(d,1), regbbox(d,3)).ne.0) then - call CCTK_WARN (0, "Internal error: array extents are not integer multiples of the strides") - end if regkext = (regbbox(d,2) - regbbox(d,1)) / regbbox(d,3) + 1 dstkfac = srcbbox(d,3) / dstbbox(d,3) srckoff = (regbbox(d,1) - srcbbox(d,1)) / dstbbox(d,3) @@ -154,7 +149,7 @@ c Quadratic (second order) interpolation call CCTK_WARN (0, "Internal error: arrays have same time") end if if (t.lt.min(t1,t2,t3)-eps .or. t.gt.max(t1,t2,t3)+eps) then - call CCTK_WARN (0, "Internal error: extrapolation in time") + call CCTK_WARN (0, "Internal error: extrapolation") end if s1fac = (t - t2) * (t - t3) / ((t1 - t2) * (t1 - t3)) |