← Back to branch summary

~ubuntu-branches/debian/wheezy/arpack/wheezy

« back to all changes in this revision

Viewing changes to LAPACK/dlacpy.f

  • Committer: Package Import Robot
  • Author(s): Sylvestre Ledru
  • Date: 2011-12-10 20:32:45 UTC
  • mfrom: (1.2.2)
  • Revision ID: package-import@ubuntu.com-20111210203245-g0fo30pqvuo92fqh
Tags: 3.0-1
http://bugs.debian.org/651351
* Switch to arpack-ng since upstream is dead.
* New upstream release
* Daniel Leidert removed from the uploaders (Closes: #651351)

Show diffs side-by-side

added added

removed removed

Lines of Context:
LAPACK/dlacpy.f
1
 
      SUBROUTINE DLACPY( UPLO, M, N, A, LDA, B, LDB )
2
 
*
3
 
*  -- LAPACK auxiliary routine (version 2.0) --
4
 
*     Univ. of Tennessee, Univ. of California Berkeley, NAG Ltd.,
5
 
*     Courant Institute, Argonne National Lab, and Rice University
6
 
*     February 29, 1992
7
 
*
8
 
*     .. Scalar Arguments ..
9
 
      CHARACTER          UPLO
10
 
      INTEGER            LDA, LDB, M, N
11
 
*     ..
12
 
*     .. Array Arguments ..
13
 
      DOUBLE PRECISION   A( LDA, * ), B( LDB, * )
14
 
*     ..
15
 
*
16
 
*  Purpose
17
 
*  =======
18
 
*
19
 
*  DLACPY copies all or part of a two-dimensional matrix A to another
20
 
*  matrix B.
21
 
*
22
 
*  Arguments
23
 
*  =========
24
 
*
25
 
*  UPLO    (input) CHARACTER*1
26
 
*          Specifies the part of the matrix A to be copied to B.
27
 
*          = 'U':      Upper triangular part
28
 
*          = 'L':      Lower triangular part
29
 
*          Otherwise:  All of the matrix A
30
 
*
31
 
*  M       (input) INTEGER
32
 
*          The number of rows of the matrix A.  M >= 0.
33
 
*
34
 
*  N       (input) INTEGER
35
 
*          The number of columns of the matrix A.  N >= 0.
36
 
*
37
 
*  A       (input) DOUBLE PRECISION array, dimension (LDA,N)
38
 
*          The m by n matrix A.  If UPLO = 'U', only the upper triangle
39
 
*          or trapezoid is accessed; if UPLO = 'L', only the lower
40
 
*          triangle or trapezoid is accessed.
41
 
*
42
 
*  LDA     (input) INTEGER
43
 
*          The leading dimension of the array A.  LDA >= max(1,M).
44
 
*
45
 
*  B       (output) DOUBLE PRECISION array, dimension (LDB,N)
46
 
*          On exit, B = A in the locations specified by UPLO.
47
 
*
48
 
*  LDB     (input) INTEGER
49
 
*          The leading dimension of the array B.  LDB >= max(1,M).
50
 
*
51
 
*  =====================================================================
52
 
*
53
 
*     .. Local Scalars ..
54
 
      INTEGER            I, J
55
 
*     ..
56
 
*     .. External Functions ..
57
 
      LOGICAL            LSAME
58
 
      EXTERNAL           LSAME
59
 
*     ..
60
 
*     .. Intrinsic Functions ..
61
 
      INTRINSIC          MIN
62
 
*     ..
63
 
*     .. Executable Statements ..
64
 
*
65
 
      IF( LSAME( UPLO, 'U' ) ) THEN
66
 
         DO 20 J = 1, N
67
 
            DO 10 I = 1, MIN( J, M )
68
 
               B( I, J ) = A( I, J )
69
 
   10       CONTINUE
70
 
   20    CONTINUE
71
 
      ELSE IF( LSAME( UPLO, 'L' ) ) THEN
72
 
         DO 40 J = 1, N
73
 
            DO 30 I = J, M
74
 
               B( I, J ) = A( I, J )
75
 
   30       CONTINUE
76
 
   40    CONTINUE
77
 
      ELSE
78
 
         DO 60 J = 1, N
79
 
            DO 50 I = 1, M
80
 
               B( I, J ) = A( I, J )
81
 
   50       CONTINUE
82
 
   60    CONTINUE
83
 
      END IF
84
 
      RETURN
85
 
*
86
 
*     End of DLACPY
87
 
*
88
 
      END