1 2 3 4 5 6 7 8 9 10 11
! { dg-do compile } ! { dg-options "-Ofast" } subroutine h3dall(z,hvec,hder,nterms) complex *16 hvec(0:1),hder(0:1) complex *16 z,zinv,ztmp/1.0/ zinv=1.0/z do i=1,nterms ztmp=zinv*i hder(i)=hvec(i-1)-ztmp*hvec(i) enddo end