summaryrefslogtreecommitdiff
path: root/2.3-1/src/fortran/blas/zcopy.f
diff options
context:
space:
mode:
Diffstat (limited to '2.3-1/src/fortran/blas/zcopy.f')
-rw-r--r--2.3-1/src/fortran/blas/zcopy.f33
1 files changed, 33 insertions, 0 deletions
diff --git a/2.3-1/src/fortran/blas/zcopy.f b/2.3-1/src/fortran/blas/zcopy.f
new file mode 100644
index 00000000..9ccfa880
--- /dev/null
+++ b/2.3-1/src/fortran/blas/zcopy.f
@@ -0,0 +1,33 @@
+ subroutine zcopy(n,zx,incx,zy,incy)
+c
+c copies a vector, x, to a vector, y.
+c jack dongarra, linpack, 4/11/78.
+c modified 12/3/93, array(1) declarations changed to array(*)
+c
+ double complex zx(*),zy(*)
+ integer i,incx,incy,ix,iy,n
+c
+ if(n.le.0)return
+ if(incx.eq.1.and.incy.eq.1)go to 20
+c
+c code for unequal increments or equal increments
+c not equal to 1
+c
+ ix = 1
+ iy = 1
+ if(incx.lt.0)ix = (-n+1)*incx + 1
+ if(incy.lt.0)iy = (-n+1)*incy + 1
+ do 10 i = 1,n
+ zy(iy) = zx(ix)
+ ix = ix + incx
+ iy = iy + incy
+ 10 continue
+ return
+c
+c code for both increments equal to 1
+c
+ 20 do 30 i = 1,n
+ zy(i) = zx(i)
+ 30 continue
+ return
+ end