1 |
guez |
3 |
subroutine scopy(n, sx, incx, sy, incy) |
2 |
|
|
|
3 |
|
|
! From cray.F, v 1.1.1.1 2004/05/19 12:53:05 |
4 |
|
|
|
5 |
|
|
! This subroutine should not exist in a Fortran 95 program. If the |
6 |
|
|
! actual arguments are of rank 1 then replace each call to this |
7 |
|
|
! subroutine by the simple statement in this subroutine. |
8 |
|
|
|
9 |
|
|
IMPLICIT NONE |
10 |
|
|
|
11 |
|
|
integer, intent(in):: n, incx, incy |
12 |
|
|
real, intent(in):: sx((n-1)*incx+1) |
13 |
|
|
real, intent(inout):: sy((n-1)*incy+1) |
14 |
|
|
|
15 |
|
|
!------------------------- |
16 |
|
|
|
17 |
|
|
sy(::incy) = sx(::incx) |
18 |
|
|
|
19 |
|
|
end subroutine scopy |
20 |
|
|
|
21 |
|
|
!*********************************** |
22 |
|
|
|
23 |
|
|
real function ssum(n,sx,incx) |
24 |
|
|
|
25 |
|
|
! From cray.F, v 1.1.1.1 2004/05/19 12:53:05 |
26 |
|
|
|
27 |
|
|
IMPLICIT NONE |
28 |
|
|
|
29 |
|
|
integer, intent(in):: n, incx |
30 |
|
|
real, intent(in):: sx((n-1)*incx+1) |
31 |
|
|
|
32 |
|
|
!----------------------- |
33 |
|
|
|
34 |
|
|
ssum=sum(sx(::incx)) |
35 |
|
|
|
36 |
|
|
end function ssum |