------- Comment #2 from spop at gcc dot gnu dot org 2009-12-18 07:43 -------
Reduced testcase:
module mcc_m
integer, parameter, private :: longreal = selected_real_kind(15,90)
contains
subroutine mutual_ind_cir_cir_coils (r1, r2, x12, y12, z12, l1, l2,
turns1, turns2, &
mu, rotate, m, l12)
real (kind = longreal), intent(out) :: l12
real (kind = longreal), dimension(1:9), save :: z1g, z2g, zw
gauss:do i = 1, 9
theta_l12 = 0.0_longreal
theta1: do n1 = 1, 2*m
theta_1 = pi*real(n1,longreal)/real(m,longreal)
theta2: do n2 = 1, 2*m
numerator = -sin(theta_1)*tvx + cos(theta_1)*tvy
theta_l12 = theta_l12 + numerator/denominator
end do theta2
end do theta1
l12 = l12 + zw(i)*theta_l12
end do gauss
l12 = coefficient * l12
end subroutine mutual_ind_cir_cir_coils
end module mcc_m
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=42180
