From cc13456d06823a3d4095217f9f011cadb80d2e4e Mon Sep 17 00:00:00 2001 From: miguel Date: Tue, 3 Apr 2001 15:38:53 +0000 Subject: Oooops! bug. Fixed it (divisions by zero). git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/LegoExcision/trunk@18 f75ba9e5-694f-0410-ac2c-87ea7ce7132b --- src/findnormals.F90 | 6 +++--- 1 file 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. -- cgit v1.2.3