diff options
| author | schnetter <schnetter@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 2005-12-12 23:57:28 +0000 |
|---|---|---|
| committer | schnetter <schnetter@41e88fdd-2190-4c69-9c84-4659c8cf322e> | 2005-12-12 23:57:28 +0000 |
| commit | 0cb3a8838bbb1a9986866ebff44591a7ead621f1 (patch) | |
| tree | bc122f0927bbbc1689077cfeaeaf6a7227c0d302 | |
| parent | f1a1704e5ae19bcac56bebd8d1772acc6b3d244a (diff) | |
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
| -rw-r--r-- | src/SourceData.F77 | 10 |
1 files 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 |