gcc/testsuite/gfortran.dg/pr80494.f90

   1 ! { dg-do compile }
   2 ! { dg-options "-std=gnu -O2" }
   3
   4 subroutine CalcCgr(C,rmax,ordgr_max)
   5   integer, intent(in) :: rmax,ordgr_max
   6   double complex :: Zadj(2,2), Zadj2(2,2)
   7   double complex, intent(out) :: C(0:rmax,0:rmax,0:rmax)
   8   double complex, allocatable :: Cexpgr(:,:,:,:)
   9   double complex :: Caux
  10   integer :: rmaxB,rmaxExp,r,n0,n1,n2,k,l,i,j,m,n,nn
  11
  12   rmaxB = 2*rmax
  13   rmaxExp = rmaxB
  14   allocate(Cexpgr(0:rmaxExp/2,0:rmaxExp,0:rmaxExp,0:ordgr_max))
  15
  16   rloop: do r=0,rmaxExp/2
  17     do n0=r,1,-1
  18       do nn=r-n0,0,-1
  19         do i=1,2
  20           Caux = Caux - Zadj(i,l)
  21         end do
  22         Cexpgr(n0,0,0,0) = Caux/(2*(nn+1))
  23       end do
  24     end do
  25     do n1=0,r
  26       n2 = r-n1
  27       if (r.le.rmax) then
  28         C(0,n1,n2) = Cexpgr(0,n1,n2,0)
  29       end if
  30     end do
  31   end do rloop
  32 end subroutine CalcCgr