From dc0bae856e8bb5d0fb68fc0aafc31f0595a99f3e Mon Sep 17 00:00:00 2001
From: Antoine Cyril David Hoffmann <ahoffman@spcpc606.epfl.ch>
Date: Mon, 15 Mar 2021 16:07:45 +0100
Subject: [PATCH] Imaginary unit was not written in double precision
---
src/prec_const_mod.F90 | 38 +++++++++++++++++++-------------------
1 file changed, 19 insertions(+), 19 deletions(-)
diff --git a/src/prec_const_mod.F90 b/src/prec_const_mod.F90
index 7282015b..27da659d 100644
--- a/src/prec_const_mod.F90
+++ b/src/prec_const_mod.F90
@@ -23,31 +23,31 @@ MODULE prec_const
! Some useful constants, to avoid recomputing then too often
- REAL(dp), PARAMETER :: PI=3.141592653589793238462643383279502884197_dp
- REAL(dp), PARAMETER :: PIO2=1.57079632679489661923132169163975144209858_dp
- REAL(dp), PARAMETER :: TWOPI=6.283185307179586476925286766559005768394_dp
- REAL(dp), PARAMETER :: SQRT2=1.41421356237309504880168872420969807856967_dp
- REAL(dp), PARAMETER :: INVSQRT2=0.7071067811865475244008443621048490392848359377_dp
- REAL(dp), PARAMETER :: SQRT3=1.73205080756887729352744634150587236694281_dp
- REAL(dp), PARAMETER :: onetwelfth=0.08333333333333333333333333333333333333333333333_dp
- REAL(dp), PARAMETER :: onetwentyfourth=0.04166666666666666666666666666666666666666666667_dp
- REAL(dp), PARAMETER :: onethird=0.33333333333333333333333333333333333333333333333_dp
- REAL(dp), PARAMETER :: onesixth=0.1666666666666666666666666666666666666666666667_dp
- REAL(dp), PARAMETER :: fivesixths=0.8333333333333333333333333333333333333333333333_dp
- REAL(dp), PARAMETER :: sevensixths=1.1666666666666666666666666666666666666666666667_dp
- REAL(dp), PARAMETER :: elevensixths=1.833333333333333333333333333333333333333333333_dp
- REAL(dp), PARAMETER :: nineeighths=1.125_dp
- REAL(dp), PARAMETER :: onesixteenth=0.0625_dp
- REAL(dp), PARAMETER :: ninesixteenths=0.5625_dp
- REAL(dp), PARAMETER :: thirteentwelfths = 1.083333333333333333333333333333333333333333333_dp
- COMPLEX(dp), PARAMETER :: imagu = (0,1)
+ REAL(dp), PARAMETER :: PI=3.141592653589793238462643383279502884197_dp
+ REAL(dp), PARAMETER :: PIO2=1.57079632679489661923132169163975144209858_dp
+ REAL(dp), PARAMETER :: TWOPI=6.283185307179586476925286766559005768394_dp
+ REAL(dp), PARAMETER :: SQRT2=1.41421356237309504880168872420969807856967_dp
+ REAL(dp), PARAMETER :: INVSQRT2=0.7071067811865475244008443621048490392848359377_dp
+ REAL(dp), PARAMETER :: SQRT3=1.73205080756887729352744634150587236694281_dp
+ REAL(dp), PARAMETER :: onetwelfth=0.08333333333333333333333333333333333333333333333_dp
+ REAL(dp), PARAMETER :: onetwentyfourth=0.04166666666666666666666666666666666666666666667_dp
+ REAL(dp), PARAMETER :: onethird=0.33333333333333333333333333333333333333333333333_dp
+ REAL(dp), PARAMETER :: onesixth=0.1666666666666666666666666666666666666666666667_dp
+ REAL(dp), PARAMETER :: fivesixths=0.8333333333333333333333333333333333333333333333_dp
+ REAL(dp), PARAMETER :: sevensixths=1.1666666666666666666666666666666666666666666667_dp
+ REAL(dp), PARAMETER :: elevensixths=1.833333333333333333333333333333333333333333333_dp
+ REAL(dp), PARAMETER :: nineeighths=1.125_dp
+ REAL(dp), PARAMETER :: onesixteenth=0.0625_dp
+ REAL(dp), PARAMETER :: ninesixteenths=0.5625_dp
+ REAL(dp), PARAMETER :: thirteentwelfths = 1.083333333333333333333333333333333333333333333_dp
+ COMPLEX(dp), PARAMETER :: imagu = (0._dp,1._dp)
CONTAINS
SUBROUTINE INIT_PREC_CONST
IMPLICIT NONE
integer :: ierr,me
-
+
REAL(sp) :: a = 1_sp
REAL(dp) :: b = 1_dp
--
GitLab