summaryrefslogtreecommitdiff
path: root/src/fortran/blas/zcopy.f
diff options
context:
space:
mode:
authorjofret2010-06-21 06:24:38 +0000
committerjofret2010-06-21 06:24:38 +0000
commit8b44229ef44f0558ce045e46ff833fb44df913c9 (patch)
tree4fdd89b2258c18c905d2f17ef4bea8ea9cf32893 /src/fortran/blas/zcopy.f
parentaf0366230e14cc75d9e9a183375ee9cb69fb46b6 (diff)
downloadscilab2c-8b44229ef44f0558ce045e46ff833fb44df913c9.tar.gz
scilab2c-8b44229ef44f0558ce045e46ff833fb44df913c9.tar.bz2
scilab2c-8b44229ef44f0558ce045e46ff833fb44df913c9.zip
Tagging the 2.0 release of scilab2crelease-2.0
Diffstat (limited to 'src/fortran/blas/zcopy.f')
-rw-r--r--src/fortran/blas/zcopy.f33
1 files changed, 33 insertions, 0 deletions
diff --git a/src/fortran/blas/zcopy.f b/src/fortran/blas/zcopy.f
new file mode 100644
index 00000000..9ccfa880
--- /dev/null
+++ b/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