ruby-lapack 1.3

data/dev/defs/zlarscl2

@@ -0,0 +1,61 @@
+--- 
+:name: zlarscl2
+:md5sum: 080c93b0555c1909ea0298bd0d4ec2f9
+:category: :subroutine
+:arguments: 
+- m: 
+    :type: integer
+    :intent: input
+- n: 
+    :type: integer
+    :intent: input
+- d: 
+    :type: doublereal
+    :intent: input
+    :dims: 
+    - m
+- x: 
+    :type: doublecomplex
+    :intent: input/output
+    :dims: 
+    - ldx
+    - n
+- ldx: 
+    :type: integer
+    :intent: input
+:substitutions: {}
+
+:fortran_help: "      SUBROUTINE ZLARSCL2 ( M, N, D, X, LDX )\n\n\
+  *  Purpose\n\
+  *  =======\n\
+  *\n\
+  *  ZLARSCL2 performs a reciprocal diagonal scaling on an vector:\n\
+  *    x <-- inv(D) * x\n\
+  *  where the DOUBLE PRECISION diagonal matrix D is stored as a vector.\n\
+  *\n\
+  *  Eventually to be replaced by BLAS_zge_diag_scale in the new BLAS\n\
+  *  standard.\n\
+  *\n\n\
+  *  Arguments\n\
+  *  =========\n\
+  *\n\
+  *     M       (input) INTEGER\n\
+  *     The number of rows of D and X. M >= 0.\n\
+  *\n\
+  *     N       (input) INTEGER\n\
+  *     The number of columns of D and X. N >= 0.\n\
+  *\n\
+  *     D       (input) DOUBLE PRECISION array, length M\n\
+  *     Diagonal matrix D, stored as a vector of length M.\n\
+  *\n\
+  *     X       (input/output) COMPLEX*16 array, dimension (LDX,N)\n\
+  *     On entry, the vector X to be scaled by D.\n\
+  *     On exit, the scaled vector.\n\
+  *\n\
+  *     LDX     (input) INTEGER\n\
+  *     The leading dimension of the vector X. LDX >= 0.\n\
+  *\n\n\
+  *  =====================================================================\n\
+  *\n\
+  *     .. Local Scalars ..\n      INTEGER            I, J\n\
+  *     ..\n"

data/dev/defs/zlartg

@@ -0,0 +1,66 @@
+--- 
+:name: zlartg
+:md5sum: 669790d67331c2cd17991b7c96ecc53d
+:category: :subroutine
+:arguments: 
+- f: 
+    :type: doublecomplex
+    :intent: input
+- g: 
+    :type: doublecomplex
+    :intent: input
+- cs: 
+    :type: doublereal
+    :intent: output
+- sn: 
+    :type: doublecomplex
+    :intent: output
+- r: 
+    :type: doublecomplex
+    :intent: output
+:substitutions: {}
+
+:fortran_help: "      SUBROUTINE ZLARTG( F, G, CS, SN, R )\n\n\
+  *  Purpose\n\
+  *  =======\n\
+  *\n\
+  *  ZLARTG generates a plane rotation so that\n\
+  *\n\
+  *     [  CS  SN  ]     [ F ]     [ R ]\n\
+  *     [  __      ]  .  [   ]  =  [   ]   where CS**2 + |SN|**2 = 1.\n\
+  *     [ -SN  CS  ]     [ G ]     [ 0 ]\n\
+  *\n\
+  *  This is a faster version of the BLAS1 routine ZROTG, except for\n\
+  *  the following differences:\n\
+  *     F and G are unchanged on return.\n\
+  *     If G=0, then CS=1 and SN=0.\n\
+  *     If F=0, then CS=0 and SN is chosen so that R is real.\n\
+  *\n\n\
+  *  Arguments\n\
+  *  =========\n\
+  *\n\
+  *  F       (input) COMPLEX*16\n\
+  *          The first component of vector to be rotated.\n\
+  *\n\
+  *  G       (input) COMPLEX*16\n\
+  *          The second component of vector to be rotated.\n\
+  *\n\
+  *  CS      (output) DOUBLE PRECISION\n\
+  *          The cosine of the rotation.\n\
+  *\n\
+  *  SN      (output) COMPLEX*16\n\
+  *          The sine of the rotation.\n\
+  *\n\
+  *  R       (output) COMPLEX*16\n\
+  *          The nonzero component of the rotated vector.\n\
+  *\n\n\
+  *  Further Details\n\
+  *  ======= =======\n\
+  *\n\
+  *  3-5-96 - Modified with a new algorithm by W. Kahan and J. Demmel\n\
+  *\n\
+  *  This version has a few statements commented out for thread safety\n\
+  *  (machine parameters are computed on each entry). 10 feb 03, SJH.\n\
+  *\n\
+  *  =====================================================================\n\
+  *\n"

data/dev/defs/zlartv

@@ -0,0 +1,82 @@
+--- 
+:name: zlartv
+:md5sum: e5ca2c7119fc89e868f72f604550a7d4
+:category: :subroutine
+:arguments: 
+- n: 
+    :type: integer
+    :intent: input
+- x: 
+    :type: doublecomplex
+    :intent: input/output
+    :dims: 
+    - 1+(n-1)*incx
+- incx: 
+    :type: integer
+    :intent: input
+- y: 
+    :type: doublecomplex
+    :intent: input/output
+    :dims: 
+    - 1+(n-1)*incy
+- incy: 
+    :type: integer
+    :intent: input
+- c: 
+    :type: doublereal
+    :intent: input
+    :dims: 
+    - 1+(n-1)*incc
+- s: 
+    :type: doublecomplex
+    :intent: input
+    :dims: 
+    - 1+(n-1)*incc
+- incc: 
+    :type: integer
+    :intent: input
+:substitutions: {}
+
+:fortran_help: "      SUBROUTINE ZLARTV( N, X, INCX, Y, INCY, C, S, INCC )\n\n\
+  *  Purpose\n\
+  *  =======\n\
+  *\n\
+  *  ZLARTV applies a vector of complex plane rotations with real cosines\n\
+  *  to elements of the complex vectors x and y. For i = 1,2,...,n\n\
+  *\n\
+  *     ( x(i) ) := (        c(i)   s(i) ) ( x(i) )\n\
+  *     ( y(i) )    ( -conjg(s(i))  c(i) ) ( y(i) )\n\
+  *\n\n\
+  *  Arguments\n\
+  *  =========\n\
+  *\n\
+  *  N       (input) INTEGER\n\
+  *          The number of plane rotations to be applied.\n\
+  *\n\
+  *  X       (input/output) COMPLEX*16 array, dimension (1+(N-1)*INCX)\n\
+  *          The vector x.\n\
+  *\n\
+  *  INCX    (input) INTEGER\n\
+  *          The increment between elements of X. INCX > 0.\n\
+  *\n\
+  *  Y       (input/output) COMPLEX*16 array, dimension (1+(N-1)*INCY)\n\
+  *          The vector y.\n\
+  *\n\
+  *  INCY    (input) INTEGER\n\
+  *          The increment between elements of Y. INCY > 0.\n\
+  *\n\
+  *  C       (input) DOUBLE PRECISION array, dimension (1+(N-1)*INCC)\n\
+  *          The cosines of the plane rotations.\n\
+  *\n\
+  *  S       (input) COMPLEX*16 array, dimension (1+(N-1)*INCC)\n\
+  *          The sines of the plane rotations.\n\
+  *\n\
+  *  INCC    (input) INTEGER\n\
+  *          The increment between elements of C and S. INCC > 0.\n\
+  *\n\n\
+  *  =====================================================================\n\
+  *\n\
+  *     .. Local Scalars ..\n      INTEGER            I, IC, IX, IY\n      COMPLEX*16         XI, YI\n\
+  *     ..\n\
+  *     .. Intrinsic Functions ..\n      INTRINSIC          DCONJG\n\
+  *     ..\n"

data/dev/defs/zlarz

@@ -0,0 +1,111 @@
+--- 
+:name: zlarz
+:md5sum: 51817751cd08b24f60b8544807b10bea
+:category: :subroutine
+:arguments: 
+- side: 
+    :type: char
+    :intent: input
+- m: 
+    :type: integer
+    :intent: input
+- n: 
+    :type: integer
+    :intent: input
+- l: 
+    :type: integer
+    :intent: input
+- v: 
+    :type: doublecomplex
+    :intent: input
+    :dims: 
+    - 1+(l-1)*abs(incv)
+- incv: 
+    :type: integer
+    :intent: input
+- tau: 
+    :type: doublecomplex
+    :intent: input
+- c: 
+    :type: doublecomplex
+    :intent: input/output
+    :dims: 
+    - ldc
+    - n
+- ldc: 
+    :type: integer
+    :intent: input
+- work: 
+    :type: doublecomplex
+    :intent: workspace
+    :dims: 
+    - "lsame_(&side,\"L\") ? n : lsame_(&side,\"R\") ? m : 0"
+:substitutions: {}
+
+:fortran_help: "      SUBROUTINE ZLARZ( SIDE, M, N, L, V, INCV, TAU, C, LDC, WORK )\n\n\
+  *  Purpose\n\
+  *  =======\n\
+  *\n\
+  *  ZLARZ applies a complex elementary reflector H to a complex\n\
+  *  M-by-N matrix C, from either the left or the right. H is represented\n\
+  *  in the form\n\
+  *\n\
+  *        H = I - tau * v * v'\n\
+  *\n\
+  *  where tau is a complex scalar and v is a complex vector.\n\
+  *\n\
+  *  If tau = 0, then H is taken to be the unit matrix.\n\
+  *\n\
+  *  To apply H' (the conjugate transpose of H), supply conjg(tau) instead\n\
+  *  tau.\n\
+  *\n\
+  *  H is a product of k elementary reflectors as returned by ZTZRZF.\n\
+  *\n\n\
+  *  Arguments\n\
+  *  =========\n\
+  *\n\
+  *  SIDE    (input) CHARACTER*1\n\
+  *          = 'L': form  H * C\n\
+  *          = 'R': form  C * H\n\
+  *\n\
+  *  M       (input) INTEGER\n\
+  *          The number of rows of the matrix C.\n\
+  *\n\
+  *  N       (input) INTEGER\n\
+  *          The number of columns of the matrix C.\n\
+  *\n\
+  *  L       (input) INTEGER\n\
+  *          The number of entries of the vector V containing\n\
+  *          the meaningful part of the Householder vectors.\n\
+  *          If SIDE = 'L', M >= L >= 0, if SIDE = 'R', N >= L >= 0.\n\
+  *\n\
+  *  V       (input) COMPLEX*16 array, dimension (1+(L-1)*abs(INCV))\n\
+  *          The vector v in the representation of H as returned by\n\
+  *          ZTZRZF. V is not used if TAU = 0.\n\
+  *\n\
+  *  INCV    (input) INTEGER\n\
+  *          The increment between elements of v. INCV <> 0.\n\
+  *\n\
+  *  TAU     (input) COMPLEX*16\n\
+  *          The value tau in the representation of H.\n\
+  *\n\
+  *  C       (input/output) COMPLEX*16 array, dimension (LDC,N)\n\
+  *          On entry, the M-by-N matrix C.\n\
+  *          On exit, C is overwritten by the matrix H * C if SIDE = 'L',\n\
+  *          or C * H if SIDE = 'R'.\n\
+  *\n\
+  *  LDC     (input) INTEGER\n\
+  *          The leading dimension of the array C. LDC >= max(1,M).\n\
+  *\n\
+  *  WORK    (workspace) COMPLEX*16 array, dimension\n\
+  *                         (N) if SIDE = 'L'\n\
+  *                      or (M) if SIDE = 'R'\n\
+  *\n\n\
+  *  Further Details\n\
+  *  ===============\n\
+  *\n\
+  *  Based on contributions by\n\
+  *    A. Petitet, Computer Science Dept., Univ. of Tenn., Knoxville, USA\n\
+  *\n\
+  *  =====================================================================\n\
+  *\n"

data/dev/defs/zlarzb

@@ -0,0 +1,150 @@
+--- 
+:name: zlarzb
+:md5sum: c84f4faa09b773e85cce145ad0d50edf
+:category: :subroutine
+:arguments: 
+- side: 
+    :type: char
+    :intent: input
+- trans: 
+    :type: char
+    :intent: input
+- direct: 
+    :type: char
+    :intent: input
+- storev: 
+    :type: char
+    :intent: input
+- m: 
+    :type: integer
+    :intent: input
+- n: 
+    :type: integer
+    :intent: input
+- k: 
+    :type: integer
+    :intent: input
+- l: 
+    :type: integer
+    :intent: input
+- v: 
+    :type: doublecomplex
+    :intent: input
+    :dims: 
+    - ldv
+    - nv
+- ldv: 
+    :type: integer
+    :intent: input
+- t: 
+    :type: doublecomplex
+    :intent: input
+    :dims: 
+    - ldt
+    - k
+- ldt: 
+    :type: integer
+    :intent: input
+- c: 
+    :type: doublecomplex
+    :intent: input/output
+    :dims: 
+    - ldc
+    - n
+- ldc: 
+    :type: integer
+    :intent: input
+- work: 
+    :type: doublecomplex
+    :intent: workspace
+    :dims: 
+    - ldwork
+    - k
+- ldwork: 
+    :type: integer
+    :intent: input
+:substitutions: 
+  ldwork: "MAX(1,n) ? side = 'l' : MAX(1,m) ? side = 'r' : 0"
+:fortran_help: "      SUBROUTINE ZLARZB( SIDE, TRANS, DIRECT, STOREV, M, N, K, L, V, LDV, T, LDT, C, LDC, WORK, LDWORK )\n\n\
+  *  Purpose\n\
+  *  =======\n\
+  *\n\
+  *  ZLARZB applies a complex block reflector H or its transpose H**H\n\
+  *  to a complex distributed M-by-N  C from the left or the right.\n\
+  *\n\
+  *  Currently, only STOREV = 'R' and DIRECT = 'B' are supported.\n\
+  *\n\n\
+  *  Arguments\n\
+  *  =========\n\
+  *\n\
+  *  SIDE    (input) CHARACTER*1\n\
+  *          = 'L': apply H or H' from the Left\n\
+  *          = 'R': apply H or H' from the Right\n\
+  *\n\
+  *  TRANS   (input) CHARACTER*1\n\
+  *          = 'N': apply H (No transpose)\n\
+  *          = 'C': apply H' (Conjugate transpose)\n\
+  *\n\
+  *  DIRECT  (input) CHARACTER*1\n\
+  *          Indicates how H is formed from a product of elementary\n\
+  *          reflectors\n\
+  *          = 'F': H = H(1) H(2) . . . H(k) (Forward, not supported yet)\n\
+  *          = 'B': H = H(k) . . . H(2) H(1) (Backward)\n\
+  *\n\
+  *  STOREV  (input) CHARACTER*1\n\
+  *          Indicates how the vectors which define the elementary\n\
+  *          reflectors are stored:\n\
+  *          = 'C': Columnwise                        (not supported yet)\n\
+  *          = 'R': Rowwise\n\
+  *\n\
+  *  M       (input) INTEGER\n\
+  *          The number of rows of the matrix C.\n\
+  *\n\
+  *  N       (input) INTEGER\n\
+  *          The number of columns of the matrix C.\n\
+  *\n\
+  *  K       (input) INTEGER\n\
+  *          The order of the matrix T (= the number of elementary\n\
+  *          reflectors whose product defines the block reflector).\n\
+  *\n\
+  *  L       (input) INTEGER\n\
+  *          The number of columns of the matrix V containing the\n\
+  *          meaningful part of the Householder reflectors.\n\
+  *          If SIDE = 'L', M >= L >= 0, if SIDE = 'R', N >= L >= 0.\n\
+  *\n\
+  *  V       (input) COMPLEX*16 array, dimension (LDV,NV).\n\
+  *          If STOREV = 'C', NV = K; if STOREV = 'R', NV = L.\n\
+  *\n\
+  *  LDV     (input) INTEGER\n\
+  *          The leading dimension of the array V.\n\
+  *          If STOREV = 'C', LDV >= L; if STOREV = 'R', LDV >= K.\n\
+  *\n\
+  *  T       (input) COMPLEX*16 array, dimension (LDT,K)\n\
+  *          The triangular K-by-K matrix T in the representation of the\n\
+  *          block reflector.\n\
+  *\n\
+  *  LDT     (input) INTEGER\n\
+  *          The leading dimension of the array T. LDT >= K.\n\
+  *\n\
+  *  C       (input/output) COMPLEX*16 array, dimension (LDC,N)\n\
+  *          On entry, the M-by-N matrix C.\n\
+  *          On exit, C is overwritten by H*C or H'*C or C*H or C*H'.\n\
+  *\n\
+  *  LDC     (input) INTEGER\n\
+  *          The leading dimension of the array C. LDC >= max(1,M).\n\
+  *\n\
+  *  WORK    (workspace) COMPLEX*16 array, dimension (LDWORK,K)\n\
+  *\n\
+  *  LDWORK  (input) INTEGER\n\
+  *          The leading dimension of the array WORK.\n\
+  *          If SIDE = 'L', LDWORK >= max(1,N);\n\
+  *          if SIDE = 'R', LDWORK >= max(1,M).\n\
+  *\n\n\
+  *  Further Details\n\
+  *  ===============\n\
+  *\n\
+  *  Based on contributions by\n\
+  *    A. Petitet, Computer Science Dept., Univ. of Tenn., Knoxville, USA\n\
+  *\n\
+  *  =====================================================================\n\
+  *\n"