2 ! PR32489 Endless loop when compiling.
3 ! Derived from fft257.f90, Public domain 2004 James Van Buskirk.
4 ! Note: The problem solved here was not an infinite loop issue. Middle-end
5 ! could not handle the array constructor unfolded by the front end.
6 ! WARNING: Potential resource hog.
7 ! Jerry DeLisle <jvdelisle@gcc.gnu.org>
10 integer, parameter :: dp
= selected_real_kind(15,300)
11 integer, parameter :: N
= 257
15 complex(dp
), parameter :: ri(2) = (/(1,0),(0,1)/)
24 h1
= (/(sum((/(exp(-2*pi
*(0,1)*mod(k
*L
,N
)/N
)*h2(L
),L
=0,N
-1)/)),k
=0,N
-1)/)