diff options
| author | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-01-12 19:14:00 +0000 |
|---|---|---|
| committer | Erik Schnetter <schnetter@cct.lsu.edu> | 2007-01-12 19:14:00 +0000 |
| commit | 74a992da87c5a312669bfffd9e6dff41dbe321ee (patch) | |
| tree | b418472c29be92c287e81000bba651307b9fc8aa /CarpetExtra | |
| parent | 2bc19d01691ddc1cd94725d330c6b6abbe24c130 (diff) | |
IDScalarWaveFO: Use double precison when calling random_number
Create random numbers explicitly in a double precision variable, since
gfortran does not support yet higher precisions for the intrinsic
random_number.
darcs-hash:20070112191430-dae7b-bf8b0a500cbf4770b588f4c4a2efa0de972934d6.gz
Diffstat (limited to 'CarpetExtra')
| -rw-r--r-- | CarpetExtra/IDScalarWaveFO/src/initialdata.F90 | 13 |
1 files changed, 9 insertions, 4 deletions
diff --git a/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 b/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 index 3530fdfbc..bd8edf30b 100644 --- a/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 +++ b/CarpetExtra/IDScalarWaveFO/src/initialdata.F90 @@ -13,6 +13,7 @@ CCTK_REAL pi parameter (pi = 3.141592653589793238462643383279502884197169399375105820974944592307816406286208998628034825342117068d0) CCTK_REAL omega + double precision rphi, rpsix, rpsiy, rpsiz integer i, j, k if (CCTK_EQUALS(initial_data, "plane")) then @@ -49,10 +50,14 @@ do k=1,cctk_lsh(3) do j=1,cctk_lsh(2) do i=1,cctk_lsh(1) - call random_number (phi(i,j,k)) - call random_number (psix(i,j,k)) - call random_number (psiy(i,j,k)) - call random_number (psiz(i,j,k)) + call random_number (rphi) + call random_number (rpsix) + call random_number (rpsiy) + call random_number (rpsiz) + phi(i,j,k) = rphi + psix(i,j,k) = rpsix + psiy(i,j,k) = rpsiy + psiz(i,j,k) = rpsiz end do end do end do |