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