summaryrefslogtreecommitdiff
path: root/gcc/testsuite/gfortran.dg/graphite/pr42180.f90
blob: bb5bc0c58db0f283afa8b6938243238f98f1bd3a (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
! { dg-options "-ffast-math -O2 -fgraphite-identity" }

module mcc_m
  integer, parameter, private :: longreal = selected_real_kind(15,90)
contains
  subroutine mutual_ind_cir_cir_coils (m, l12)
    real (kind = longreal), intent(out) :: l12
    real (kind = longreal), dimension(1:9), save :: 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
! { dg-final { cleanup-modules "mcc_m" } }