aboutsummaryrefslogtreecommitdiff
diff options
context:
space:
mode:
authorschnetter <schnetter@0a4070d5-58f5-498f-b6c0-2693e757fa0f>2010-03-16 22:18:33 +0000
committerschnetter <schnetter@0a4070d5-58f5-498f-b6c0-2693e757fa0f>2010-03-16 22:18:33 +0000
commit41510cd8872b2d1b88cd859c66bac55ecf19685b (patch)
treeb3fd43298f897fd909cbed7bf0f1704bfa0de6e9
parent1a551590f45aa3ee8ba33717bcc031cdb573a736 (diff)
Pass double precision instead of single precision value
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinInitialData/IDAxiBrillBH/trunk@108 0a4070d5-58f5-498f-b6c0-2693e757fa0f
-rw-r--r--src/shmgp.F7712
1 files changed, 7 insertions, 5 deletions
diff --git a/src/shmgp.F77 b/src/shmgp.F77
index f23e6d7..a269862 100644
--- a/src/shmgp.F77
+++ b/src/shmgp.F77
@@ -256,6 +256,8 @@ c .....................................................................
CCTK_REAL ru(idi),rd(idi),rc(idi)
CCTK_REAL resid(0:40),confac(0:40)
common /io/ linp,lout
+ CCTK_REAL one
+ parameter (one=1)
c
c-time tsu0=second()
if(iskip.eq.0) then
@@ -313,10 +315,10 @@ c definition below. The call passes three arguments too many.
c Maybe the three arguments pu(nic),pd(nic),pc(nic) should be
c omitted?
call CCTK_WARN (CCTK_WARN_ABORT, "Trying to find out whether this line is reached")
- 6 call ursrhs(f(n5),f(n5c),pu(nic),pd(nic),pc(nic),ru(nic),
- + rd(nic),rc(nic),im,jm,jmc,m,k+1,jr,irc)
-c 6 call ursrhs(f(n5),f(n5c),ru(nic),rd(nic),rc(nic),
-c + im,jm,jmc,m,k+1,jr,irc)
+c 6 call ursrhs(f(n5),f(n5c),pu(nic),pd(nic),pc(nic),ru(nic),
+c + rd(nic),rc(nic),im,jm,jmc,m,k+1,jr,irc)
+ 6 call ursrhs(f(n5),f(n5c),ru(nic),rd(nic),rc(nic),
+ + im,jm,jmc,m,k+1,jr,irc)
endif
c-time tsu1=second()
c-time write(lout,*) ' time for setup =',tsu1-tsu0
@@ -413,7 +415,7 @@ cdir$ noinline
call urscal(
+ ac(n5),aw(n5),as(n5),ae(n5),an(n5),asw(n9),ase(n9),
+ ane(n9),anw(n9),q(n5),f(n5),f(n5c),q(n5c),rc(nic),
- + im,jm,jmc,ifd,i9,j9,k,m,jr,1.,rmax,ipc,irc)
+ + im,jm,jmc,ifd,i9,j9,k,m,jr,one,rmax,ipc,irc)
resid(mcyc)=rmax
confac(mcyc)=rmax/rmaxo
nb=0