5 DOUBLE PRECISION DX
(9),DY
(9)
10 DATA DX
/-2, -1, -3, -4, 1, 2, 10, 15, 14/
11 DATA DY
/0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0/
12 CALL DAXPY
(N
,DA
,DX
,DY
)
14 if (DX
(I
).ne
.DY
(I
)) call abort
19 SUBROUTINE DAXPY
(N
,DA
,DX
,DY
)
22 DOUBLE PRECISION DX
(*),DY
(*)
29 DY
(I
) = DY
(I
) + DA*DX
(I
)
34 DY
(I
) = DY
(I
) + DA*DX
(I
)
35 DY
(I
+1) = DY
(I
+1) + DA*DX
(I
+1)
36 DY
(I
+2) = DY
(I
+2) + DA*DX
(I
+2)
37 DY
(I
+3) = DY
(I
+3) + DA*DX
(I
+3)