1! { dg-options "-fgraphite-identity -g -O3 -ffast-math" }
2MODULE erf_fn
3CONTAINS
4  SUBROUTINE CALERF(ARG,RESULT,JINT)
5    DIMENSION A(5),B(4),C(9),D(8),P(6),Q(5)
6    IF (Y <= THRESH) THEN
7       DO I = 1, 3
8          XNUM = (XNUM + A(I)) * YSQ
9          XDEN = (XDEN + B(I)) * YSQ
10       END DO
11       RESULT = X * (XNUM + A(4)) / (XDEN + B(4))
12    END IF
13  END SUBROUTINE CALERF
14END MODULE erf_fn
15