version 1.7, 2011/07/22 07:38:01
|
version 1.13, 2017/06/17 10:53:42
|
Line 1
|
Line 1
|
|
*> \brief \b DCOPY |
|
* |
|
* =========== DOCUMENTATION =========== |
|
* |
|
* Online html documentation available at |
|
* http://www.netlib.org/lapack/explore-html/ |
|
* |
|
* Definition: |
|
* =========== |
|
* |
|
* SUBROUTINE DCOPY(N,DX,INCX,DY,INCY) |
|
* |
|
* .. Scalar Arguments .. |
|
* INTEGER INCX,INCY,N |
|
* .. |
|
* .. Array Arguments .. |
|
* DOUBLE PRECISION DX(*),DY(*) |
|
* .. |
|
* |
|
* |
|
*> \par Purpose: |
|
* ============= |
|
*> |
|
*> \verbatim |
|
*> |
|
*> DCOPY copies a vector, x, to a vector, y. |
|
*> uses unrolled loops for increments equal to one. |
|
*> \endverbatim |
|
* |
|
* Authors: |
|
* ======== |
|
* |
|
*> \author Univ. of Tennessee |
|
*> \author Univ. of California Berkeley |
|
*> \author Univ. of Colorado Denver |
|
*> \author NAG Ltd. |
|
* |
|
*> \date December 2016 |
|
* |
|
*> \ingroup double_blas_level1 |
|
* |
|
*> \par Further Details: |
|
* ===================== |
|
*> |
|
*> \verbatim |
|
*> |
|
*> jack dongarra, linpack, 3/11/78. |
|
*> modified 12/3/93, array(1) declarations changed to array(*) |
|
*> \endverbatim |
|
*> |
|
* ===================================================================== |
SUBROUTINE DCOPY(N,DX,INCX,DY,INCY) |
SUBROUTINE DCOPY(N,DX,INCX,DY,INCY) |
|
* |
|
* -- Reference BLAS level1 routine (version 3.7.0) -- |
|
* -- Reference BLAS is a software package provided by Univ. of Tennessee, -- |
|
* -- Univ. of California Berkeley, Univ. of Colorado Denver and NAG Ltd..-- |
|
* December 2016 |
|
* |
* .. Scalar Arguments .. |
* .. Scalar Arguments .. |
INTEGER INCX,INCY,N |
INTEGER INCX,INCY,N |
* .. |
* .. |
Line 6
|
Line 63
|
DOUBLE PRECISION DX(*),DY(*) |
DOUBLE PRECISION DX(*),DY(*) |
* .. |
* .. |
* |
* |
* Purpose |
|
* ======= |
|
* |
|
* DCOPY copies a vector, x, to a vector, y. |
|
* uses unrolled loops for increments equal to one. |
|
* |
|
* Further Details |
|
* =============== |
|
* |
|
* jack dongarra, linpack, 3/11/78. |
|
* modified 12/3/93, array(1) declarations changed to array(*) |
|
* |
|
* ===================================================================== |
* ===================================================================== |
* |
* |
* .. Local Scalars .. |
* .. Local Scalars .. |
Line 40
|
Line 85
|
DY(I) = DX(I) |
DY(I) = DX(I) |
END DO |
END DO |
IF (N.LT.7) RETURN |
IF (N.LT.7) RETURN |
END IF |
END IF |
MP1 = M + 1 |
MP1 = M + 1 |
DO I = MP1,N,7 |
DO I = MP1,N,7 |
DY(I) = DX(I) |
DY(I) = DX(I) |
Line 51
|
Line 96
|
DY(I+5) = DX(I+5) |
DY(I+5) = DX(I+5) |
DY(I+6) = DX(I+6) |
DY(I+6) = DX(I+6) |
END DO |
END DO |
ELSE |
ELSE |
* |
* |
* code for unequal increments or equal increments |
* code for unequal increments or equal increments |
* not equal to 1 |
* not equal to 1 |