y_eq_x.F
c**********************************************************************
#include "author.inc"
c* $Id: y_eq_x.F,v 1.5 1995/07/30 17:49:38 turner Exp $
c*
c* Copies vector x into vector y (y=x).
c*
c* <PARAMETER LIST>
c*
c* Input:
c* n - number of elements
c* x - vector to be copied
c*
c* Output:
c* y - destination vector
c* status - return status
c*
#include "copyright.inc"
c**********************************************************************
subroutine JT_y_eq_x (n, x, y, status)
implicit none
c
c ... Input:
integer n
real x(n)
c
c ... Output:
integer status
real y(n)
c
#ifdef use_blas
call SCOPY (n, x, 1, y, 1)
#else
c ... Local:
integer i
c
do i=1,n
y(i) = x(i)
enddo
#endif
c
status = 0
return
end