libf/dyn3d/ssum.f90

subroutine scopy(n, sx, incx, sy, incy)

  ! From cray.F, v 1.1.1.1 2004/05/19 12:53:05

  ! This subroutine should not exist in a Fortran 95 program. If the
  ! actual arguments are of rank 1 then replace each call to this
  ! subroutine by the simple statement in this subroutine.

  IMPLICIT NONE

  integer, intent(in):: n, incx, incy
  real, intent(in):: sx((n-1)*incx+1)
  real, intent(inout):: sy((n-1)*incy+1)

  !-------------------------

  sy(::incy) = sx(::incx)

end subroutine scopy

!***********************************

real function ssum(n,sx,incx)

  ! From cray.F, v 1.1.1.1 2004/05/19 12:53:05

  IMPLICIT NONE

  integer, intent(in):: n, incx
  real, intent(in):: sx((n-1)*incx+1)

  !-----------------------

  ssum=sum(sx(::incx))

end function ssum
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