diff options
author | miguel <miguel@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2001-04-03 15:38:53 +0000 |
---|---|---|
committer | miguel <miguel@f75ba9e5-694f-0410-ac2c-87ea7ce7132b> | 2001-04-03 15:38:53 +0000 |
commit | cc13456d06823a3d4095217f9f011cadb80d2e4e (patch) | |
tree | 1d010cdce735ed4cd6662d3b8b79e82ca502e0f3 | |
parent | 78ba61b65fb29a4b0e0217fc51fbdc59c6fb9487 (diff) |
Oooops! bug. Fixed it (divisions by zero).
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@18 f75ba9e5-694f-0410-ac2c-87ea7ce7132b
-rw-r--r-- | src/findnormals.F90 | 6 |
1 files changed, 3 insertions, 3 deletions
diff --git a/src/findnormals.F90 b/src/findnormals.F90 index 85f56c1..e9aebfe 100644 --- a/src/findnormals.F90 +++ b/src/findnormals.F90 @@ -94,9 +94,9 @@ subroutine excision_findnormals (ierr, mask, dirx, diry, dirz, ni, nj, nk) ! To do this, I add 1/2 to all components (with the ! correct sign), ang get the integer part. - sx = dble(int(dabs(sx)+0.5D0))*sx/dabs(sx) - sy = dble(int(dabs(sy)+0.5D0))*sy/dabs(sy) - sz = dble(int(dabs(sz)+0.5D0))*sz/dabs(sz) + if (sx /= 0.0D0) sx = dble(int(dabs(sx)+0.5D0))*sx/dabs(sx) + if (sy /= 0.0D0) sy = dble(int(dabs(sy)+0.5D0))*sy/dabs(sy) + if (sz /= 0.0D0) sz = dble(int(dabs(sz)+0.5D0))*sz/dabs(sz) ! Copy the vector (sx,sy,sz) into the normal arrays. |