diff options
| author | eschnett <eschnett@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-11 22:46:36 +0000 |
|---|---|---|
| committer | eschnett <eschnett@c83d129a-5a75-4d5a-9c4d-ed3a5855bf45> | 2012-05-11 22:46:36 +0000 |
| commit | 358d5165056ced200988620b63d158a14c1f0d3f (patch) | |
| tree | 084e704a6708b5b548f6c0562bf257adb146ce4a /src/GRHydro_Eigenproblem_Marquina.F90 | |
| parent | af6dfaa5b36603dc81d04e3c5a8e054f64c37293 (diff) | |
Support real*16 (and real*4) in GRHydro
Support using precisions other than real*8 in GRHydro by removing all
explicit references to double precision functions and constants, and
using type-generic functions and constants instead.
In particular, use "one" and "half" as constants in some places, and
use "abs", "max" etc. instead of "dabs", "dmax" etc.
git-svn-id: http://svn.einsteintoolkit.org/cactus/EinsteinEvolve/GRHydro/trunk@319 c83d129a-5a75-4d5a-9c4d-ed3a5855bf45
Diffstat (limited to 'src/GRHydro_Eigenproblem_Marquina.F90')
| -rw-r--r-- | src/GRHydro_Eigenproblem_Marquina.F90 | 18 |
1 files changed, 9 insertions, 9 deletions
diff --git a/src/GRHydro_Eigenproblem_Marquina.F90 b/src/GRHydro_Eigenproblem_Marquina.F90 index 970f7ac..9aa3765 100644 --- a/src/GRHydro_Eigenproblem_Marquina.F90 +++ b/src/GRHydro_Eigenproblem_Marquina.F90 @@ -200,11 +200,11 @@ subroutine eigenproblem_marquina(handle,rhor,velxr,velyr,& !!$ FINAL - lam1 = dmax1(dabs(lam1l),dabs(lam1r)) + lam1 = max(abs(lam1l),abs(lam1r)) lam2 = lam1 lam3 = lam1 - lamp = dmax1(dabs(lampl),dabs(lampr)) - lamm = dmax1(dabs(lamml),dabs(lammr)) + lamp = max(abs(lampl),abs(lampr)) + lamm = max(abs(lamml),abs(lammr)) !!$ lam(1) = lamm !!$ lam(2) = lam1 @@ -918,11 +918,11 @@ subroutine eigenproblem_marquina_hot(handle,rhor,velxr,velyr,& !!$ FINAL - lam1 = dmax1(dabs(lam1l),dabs(lam1r)) + lam1 = max(abs(lam1l),abs(lam1r)) lam2 = lam1 lam3 = lam1 - lamp = dmax1(dabs(lampl),dabs(lampr)) - lamm = dmax1(dabs(lamml),dabs(lammr)) + lamp = max(abs(lampl),abs(lampr)) + lamm = max(abs(lamml),abs(lammr)) !!$ lam(1) = lamm !!$ lam(2) = lam1 @@ -1610,11 +1610,11 @@ subroutine eigenproblem_marquina_general(& !!$ FINAL - lam1 = dmax1(dabs(lam1l),dabs(lam1r)) + lam1 = max(abs(lam1l),abs(lam1r)) lam2 = lam1 lam3 = lam1 - lamp = dmax1(dabs(lampl),dabs(lampr)) - lamm = dmax1(dabs(lamml),dabs(lammr)) + lamp = max(abs(lampl),abs(lampr)) + lamm = max(abs(lamml),abs(lammr)) !!$ lam(1) = lamm !!$ lam(2) = lam1 |