------- Comment #4 from pault at gcc dot gnu dot org 2008-01-11 19:12 ------- (In reply to comment #3) Breaking up the expression for h1, thusly:
hh = (/(exp(-2*pi*(0,1)*mod(k*L,N)/N)*h2(L),L=0,N-1)/)
h1 = (/(sum(hh),k=0,N-1)/)
cures the problem.
Paul
--
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=32489
