From 0cb3a8838bbb1a9986866ebff44591a7ead621f1 Mon Sep 17 00:00:00 2001 From: schnetter Date: Mon, 12 Dec 2005 23:57:28 +0000 Subject: Clean up floating point constants. Make all floating point constants double precision. git-svn-id: http://svn.cactuscode.org/arrangements/CactusWave/IDScalarWaveElliptic/trunk@71 41e88fdd-2190-4c69-9c84-4659c8cf322e --- src/SourceData.F77 | 10 +++++----- 1 file changed, 5 insertions(+), 5 deletions(-) diff --git a/src/SourceData.F77 b/src/SourceData.F77 index d5fc8da..17b3efb 100644 --- a/src/SourceData.F77 +++ b/src/SourceData.F77 @@ -60,7 +60,7 @@ c Get variable indices for all grid functions c Set up all coefficients and initial guess for solution - pi = 4.0*atan(1.0) + pi = 4.0d0*atan(1.0d0) charge_factor = 4.0d0*pi*charge*3.0d0/(4.0d0*pi*radius**3) do k=1, cctk_lsh(3) @@ -74,7 +74,7 @@ c Set up all coefficients and initial guess for solution if (r(i,j,k) <= radius) then Ncoeff(i,j,k) = charge_factor else - Ncoeff(i,j,k) = 0d0 + Ncoeff(i,j,k) = 0.0d0 end if end do @@ -95,9 +95,9 @@ c Set tolerance for stopping elliptic solve c Set any options needed for the elliptic solve - call Ell_SetStrKey(ierr, "yes","EllLinFlat::Bnd::Robin") - call Ell_SetRealKey(ierr, 0d0, "EllLinFlat::Bnd::Robin::inf") - call Ell_SetIntKey (ierr, 1, "EllLinFlat::Bnd::Robin::falloff") + call Ell_SetStrKey (ierr, "yes", "EllLinFlat::Bnd::Robin") + call Ell_SetRealKey(ierr, 0.0d0, "EllLinFlat::Bnd::Robin::inf") + call Ell_SetIntKey (ierr, 1, "EllLinFlat::Bnd::Robin::falloff") c Set parameters for specific solvers if (CCTK_EQUALS(solver,"sor")) then -- cgit v1.2.3