gkowal/amun-code
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

INTERPOLATIONS, REFINEMENT: Convert real precision correctly.

Browse Source 
Signed-off-by: Grzegorz Kowal <grzegorz@amuncode.org>
This commit is contained in:
Grzegorz Kowal 2020-08-06 21:39:10 -03:00
parent 036a887f39
commit 9dbc4e3430
3 changed files with 6 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

2
CMakeLists.txt
View File

@ -19,7 +19,7 @@ if(NOT CMAKE_BUILD_TYPE AND NOT CMAKE_CONFIGURATION_TYPES)
endif()
if(CMAKE_Fortran_COMPILER_ID MATCHES "GNU")
add_compile_options("$<$<CONFIG:RELEASE>:-march=native;-pipe;-pedantic;-ftree-vectorize;-fno-unsafe-math-optimizations;-frounding-math;-fsignaling-nans;-finline-limit=10000;-fdiagnostics-color=always;-Wunused-variable>")
add_compile_options("$<$<CONFIG:RELEASE>:-march=native;-pipe;-pedantic;-ftree-vectorize;-fno-unsafe-math-optimizations;-frounding-math;-fsignaling-nans;-finline-limit=10000;-fdiagnostics-color=always>")
add_compile_options("$<$<CONFIG:DEBUG>:-Og;-pedantic;-W;-Wall>")
endif()

4
sources/interpolations.F90
View File

@ -775,7 +775,7 @@ module interpolations
!
do j = 1, NDIMS
do i = 1, 2
ugp(1:dgp,i,j) = matmul(xgp(1:dgp,i,j), inv(1:dgp,1:dgp))
ugp(1:dgp,i,j) = real(matmul(xgp(1:dgp,i,j), inv(1:dgp,1:dgp)), kind=8)
end do
end do
@ -4791,7 +4791,7 @@ module interpolations
! prepare the interpolation coefficients vector
!
cgp(1:ngp) = matmul(xgp(1:ngp), agp(1:ngp,1:ngp))
cgp(1:ngp) = real(matmul(xgp(1:ngp), agp(1:ngp,1:ngp)), kind=8)
end if

6
sources/refinement.F90
View File

@ -398,7 +398,7 @@ module refinement
! local variables
!
integer :: p
real(kind=4) :: cref
real(kind=8) :: cref
!
!-------------------------------------------------------------------------------
!
@ -496,14 +496,14 @@ module refinement
! return variable
!
real(kind=4) :: error
real(kind=8) :: error
! local variables
!
integer :: i, im1, ip1
integer :: j, jm1, jp1
integer :: k = 1, km1, kp1
real(kind=4) :: fl, fr, fc, fx, fy, fz
real(kind=8) :: fl, fr, fc, fx, fy, fz
! local parameters
!