procyberian/gcc
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

re PR middle-end/33794 (Wrong code w/ -ffast-math)

Browse source 
PR middle-end/33794
        * gfortran.dg/pr33794.f90: New testcase.


Co-Authored-By: Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>

From-SVN: r129410
This commit is contained in:
Uros Bizjak 2007-10-17 19:09:58 +02:00 • committed by Uros Bizjak
parent 590e9a42f4
commit 2e23a5b3e3
2 changed files with 56 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

6
gcc/testsuite/ChangeLog
View file

@ -1,3 +1,9 @@
2007-10-17 Uros Bizjak <ubizjak@gmail.com>
Francois-Xavier Coudert <fxcoudert@gcc.gnu.org>
PR middle-end/33794
* gfortran.dg/pr33794.f90: New testcase.
2007-10-17 Uros Bizjak <ubizjak@gmail.com>
PR target/32961

50
gcc/testsuite/gfortran.dg/pr33794.f90 Normal file
View file

@ -0,0 +1,50 @@
! { dg-do run }
! { dg-options "-O2 -ffast-math -mfpmath=387" { target { { i?86-*-* x86_64-*-* } && lp64 } } }
! { dg-options "-O2 -ffast-math" }
module scc_m
implicit none
integer, parameter :: dp = selected_real_kind(15,90)
contains
subroutine self_ind_cir_coil (r, l, turns, mu, self_l)
implicit none
real (kind = dp), intent(in) :: r, l, turns, mu
real (kind = dp), intent(out) :: self_l
real (kind = dp) :: alpha, modulus, pk, ak, bk, ae, be, elliptice, elliptick
real (kind = dp) :: expected
alpha = atan(2.0_dp*r/l)
modulus = sin(alpha)
pk = 1.0_dp - modulus**2
ak = (((0.01451196212_dp*pk+0.03742563713_dp)*pk+ &
0.03590092383_dp)*pk+0.09666344259_dp)*pk+1.38629436112_dp
bk = (((0.00441787012_dp*pk+0.03328355346_dp)*pk+ &
0.06880248576_dp)*pk+0.12498593597_dp)*pk+0.5_dp
elliptick = ak - bk * log(pk)
ae = (((0.01736506451_dp*pk+0.04757383546_dp)*pk+ &
0.0626060122_dp)*pk+0.44325141463_dp)*pk+1.0_dp
be = (((0.00526449639_dp*pk+0.04069697526_dp)*pk+ &
0.09200180037_dp)*pk+0.2499836831_dp)*pk
elliptice = ae - be * log(pk)
self_l = (mu * turns**2 * l**2 * 2.0_dp * r)/3.0_dp * &
(((tan(alpha)**2-1.0_dp)*elliptice+elliptick)/sin(alpha) - &
tan(alpha)**2)
expected = 3.66008420600434162E-002_dp
if (abs(self_l - expected) / expected > 1e-3) &
call abort
end subroutine self_ind_cir_coil
end module scc_m
program test
use scc_m
implicit none
real (kind = dp) :: mu, turns, r, l, self_l
mu = 1.25663706143591729E-006_dp
turns = 166666.66666666666_dp
l = 3.00000000000000006E-003_dp
r = 2.99999999999999989E-002_dp
call self_ind_cir_coil (r, l, turns, mu, self_l)
end program test
! { dg-final { cleanup-modules "scc_m" } }