summaryrefslogtreecommitdiff
path: root/src/lib/lapack/dlacpy.f
diff options
context:
space:
mode:
Diffstat (limited to 'src/lib/lapack/dlacpy.f')
-rw-r--r--src/lib/lapack/dlacpy.f87
1 files changed, 0 insertions, 87 deletions
diff --git a/src/lib/lapack/dlacpy.f b/src/lib/lapack/dlacpy.f
deleted file mode 100644
index d72603a5..00000000
--- a/src/lib/lapack/dlacpy.f
+++ /dev/null
@@ -1,87 +0,0 @@
- SUBROUTINE DLACPY( UPLO, M, N, A, LDA, B, LDB )
-*
-* -- LAPACK auxiliary routine (version 3.1) --
-* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd..
-* November 2006
-*
-* .. Scalar Arguments ..
- CHARACTER UPLO
- INTEGER LDA, LDB, M, N
-* ..
-* .. Array Arguments ..
- DOUBLE PRECISION A( LDA, * ), B( LDB, * )
-* ..
-*
-* Purpose
-* =======
-*
-* DLACPY copies all or part of a two-dimensional matrix A to another
-* matrix B.
-*
-* Arguments
-* =========
-*
-* UPLO (input) CHARACTER*1
-* Specifies the part of the matrix A to be copied to B.
-* = 'U': Upper triangular part
-* = 'L': Lower triangular part
-* Otherwise: All of the matrix A
-*
-* M (input) INTEGER
-* The number of rows of the matrix A. M >= 0.
-*
-* N (input) INTEGER
-* The number of columns of the matrix A. N >= 0.
-*
-* A (input) DOUBLE PRECISION array, dimension (LDA,N)
-* The m by n matrix A. If UPLO = 'U', only the upper triangle
-* or trapezoid is accessed; if UPLO = 'L', only the lower
-* triangle or trapezoid is accessed.
-*
-* LDA (input) INTEGER
-* The leading dimension of the array A. LDA >= max(1,M).
-*
-* B (output) DOUBLE PRECISION array, dimension (LDB,N)
-* On exit, B = A in the locations specified by UPLO.
-*
-* LDB (input) INTEGER
-* The leading dimension of the array B. LDB >= max(1,M).
-*
-* =====================================================================
-*
-* .. Local Scalars ..
- INTEGER I, J
-* ..
-* .. External Functions ..
- LOGICAL LSAME
- EXTERNAL LSAME
-* ..
-* .. Intrinsic Functions ..
- INTRINSIC MIN
-* ..
-* .. Executable Statements ..
-*
- IF( LSAME( UPLO, 'U' ) ) THEN
- DO 20 J = 1, N
- DO 10 I = 1, MIN( J, M )
- B( I, J ) = A( I, J )
- 10 CONTINUE
- 20 CONTINUE
- ELSE IF( LSAME( UPLO, 'L' ) ) THEN
- DO 40 J = 1, N
- DO 30 I = J, M
- B( I, J ) = A( I, J )
- 30 CONTINUE
- 40 CONTINUE
- ELSE
- DO 60 J = 1, N
- DO 50 I = 1, M
- B( I, J ) = A( I, J )
- 50 CONTINUE
- 60 CONTINUE
- END IF
- RETURN
-*
-* End of DLACPY
-*
- END