#2096: SphericalHarmonicRecon and SphericalHarmonicReconGen tests fail on Jenkins build machine ------------------------------------+--------------------------------------- Reporter: hinder | Owner: rhaas Type: defect | Status: assigned Priority: major | Milestone: ET_2018_02 Component: EinsteinToolkit thorn | Version: development version Resolution: | Keywords: ------------------------------------+---------------------------------------
Comment (by Yosef):
This is a gfortran compiler bug triggered by the compiler options "-O2 -ffast-math -fno-finite-math-only". Remove any one option, and the bug goes away. Here's a test code
program test implicit none integer, parameter :: nx = 5 integer, parameter :: ny = 5 integer, parameter :: wp = kind(1.0d0)
complex(kind=wp), dimension(nx,ny):: st_z
integer:: i,j
do j=1, ny do i=1, nx st_z(i,j) = dcmplx(-i,-j) end do end do
write(*,*) "Test output should not be zero ", minval(abs(st_z - dcmplx(1,1))) end program test