gcc/testsuite/gfortran.dg/graphite/pr42180.f90

   1 module mcc_m
   2   integer, parameter, private :: longreal = selected_real_kind(15,90)
   3 contains
   4   subroutine mutual_ind_cir_cir_coils (m, l12)
   5     real (kind = longreal), intent(out) :: l12
   6     real (kind = longreal), dimension(1:9), save :: zw
   7     gauss:do i = 1, 9
   8        theta_l12 = 0.0_longreal
   9        theta1:   do n1 = 1, 2*m
  10           theta_1 = pi*real(n1,longreal)/real(m,longreal)
  11           theta2:       do n2 = 1, 2*m
  12              numerator = -sin(theta_1)*tvx + cos(theta_1)*tvy
  13              theta_l12 = theta_l12 + numerator/denominator
  14           end do theta2
  15        end do theta1
  16        l12 = l12 + zw(i)*theta_l12
  17     end do gauss
  18     l12 = coefficient * l12
  19   end subroutine mutual_ind_cir_cir_coils
  20 end module mcc_m