cunml2 (3p)

Name

cunml2 - torization determined by cgelqf (unblocked algorithm)

Synopsis

SUBROUTINE CUNML2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK,
INFO)

CHARACTER*1 SIDE, TRANS
COMPLEX A(LDA,*), TAU(*), C(LDC,*), WORK(*)
INTEGER M, N, K, LDA, LDC, INFO

SUBROUTINE CUNML2_64(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK,
INFO)

CHARACTER*1 SIDE, TRANS
COMPLEX A(LDA,*), TAU(*), C(LDC,*), WORK(*)
INTEGER*8 M, N, K, LDA, LDC, INFO




F95 INTERFACE
SUBROUTINE UNML2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
WORK, INFO)

CHARACTER(LEN=1) :: SIDE, TRANS
COMPLEX, DIMENSION(:) :: TAU, WORK
COMPLEX, DIMENSION(:,:) :: A, C
INTEGER :: M, N, K, LDA, LDC, INFO

SUBROUTINE UNML2_64(SIDE, TRANS, M, N, K, A, LDA, TAU, C,
LDC, WORK, INFO)

CHARACTER(LEN=1) :: SIDE, TRANS
COMPLEX, DIMENSION(:) :: TAU, WORK
COMPLEX, DIMENSION(:,:) :: A, C
INTEGER(8) :: M, N, K, LDA, LDC, INFO




C INTERFACE
#include <sunperf.h>

void cunml2(char side, char trans, int m, int n, int k, complex *a, int
lda, complex *tau, complex *c, int ldc, int *info);

void  cunml2_64(char  side, char trans, long m, long n, long k, complex
*a, long lda,  complex  *tau,  complex  *c,  long  ldc,  long
*info);

Description

Oracle Solaris Studio Performance Library                           cunml2(3P)



NAME
       cunml2 - multiply a general matrix by the unitary matrix from a LQ fac-
       torization determined by cgelqf (unblocked algorithm)


SYNOPSIS
       SUBROUTINE CUNML2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK,
             INFO)

       CHARACTER*1 SIDE, TRANS
       COMPLEX A(LDA,*), TAU(*), C(LDC,*), WORK(*)
       INTEGER M, N, K, LDA, LDC, INFO

       SUBROUTINE CUNML2_64(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC, WORK,
             INFO)

       CHARACTER*1 SIDE, TRANS
       COMPLEX A(LDA,*), TAU(*), C(LDC,*), WORK(*)
       INTEGER*8 M, N, K, LDA, LDC, INFO




   F95 INTERFACE
       SUBROUTINE UNML2(SIDE, TRANS, M, N, K, A, LDA, TAU, C, LDC,
              WORK, INFO)

       CHARACTER(LEN=1) :: SIDE, TRANS
       COMPLEX, DIMENSION(:) :: TAU, WORK
       COMPLEX, DIMENSION(:,:) :: A, C
       INTEGER :: M, N, K, LDA, LDC, INFO

       SUBROUTINE UNML2_64(SIDE, TRANS, M, N, K, A, LDA, TAU, C,
              LDC, WORK, INFO)

       CHARACTER(LEN=1) :: SIDE, TRANS
       COMPLEX, DIMENSION(:) :: TAU, WORK
       COMPLEX, DIMENSION(:,:) :: A, C
       INTEGER(8) :: M, N, K, LDA, LDC, INFO




   C INTERFACE
       #include <sunperf.h>

       void cunml2(char side, char trans, int m, int n, int k, complex *a, int
                 lda, complex *tau, complex *c, int ldc, int *info);

       void  cunml2_64(char  side, char trans, long m, long n, long k, complex
                 *a, long lda,  complex  *tau,  complex  *c,  long  ldc,  long
                 *info);



PURPOSE
       cunml2 overwrites the general complex M-by-N matrix C with
              Q * C  if SIDE = 'L' and TRANS = 'N', or

              Q**H* C  if SIDE = 'L' and TRANS = 'C', or

              C * Q  if SIDE = 'R' and TRANS = 'N', or

              C * Q**H if SIDE = 'R' and TRANS = 'C',

       where Q is a complex unitary matrix defined as the product of K elemen-
       tary reflectors

              Q = H(K)**H . . . H(2)**H * H(1)**H
       as returned by CGELQF. Q is of order M if SIDE = 'L' and of order N  if
       SIDE = 'R'.


ARGUMENTS
       SIDE (input)
                 = 'L': apply Q or Q' from the Left
                 = 'R': apply Q or Q' from the Right


       TRANS (input)
                 = 'N': apply Q  (No transpose)
                 = 'C': apply Q' (Conjugate transpose)


       M (input) The number of rows of the matrix C. M >= 0.


       N (input) The number of columns of the matrix C. N >= 0.


       K (input) The number of elementary reflectors whose product defines the
                 matrix Q.  If SIDE = 'L', M >= K >= 0; if SIDE = 'R', N >=  K
                 >= 0.


       A (input) (LDA,M)  if  SIDE  =  'L', (LDA,N) if SIDE = 'R' The i-th row
                 must contain the vector which defines the elementary  reflec-
                 tor  H(i),  for  i  = 1,2,...,k, as returned by CGELQF in the
                 first k rows of its array argument A.  A is modified  by  the
                 routine but restored on exit.


       LDA (input)
                 The leading dimension of the array A. LDA >= max(1,K).


       TAU (input)
                 TAU(i)  must  contain  the  scalar  factor  of the elementary
                 reflector H(i), as returned by CGELQF.


       C (input/output)
                 On entry, the m-by-n matrix C.  On exit, C is overwritten  by
                 Q*C or Q'*C or C*Q' or C*Q.


       LDC (input)
                 The leading dimension of the array C. LDC >= max(1,M).


       WORK (workspace)
                 (N) if SIDE = 'L', (M) if SIDE = 'R'


       INFO (output)
                 = 0: successful exit
                 < 0: if INFO = -i, the i-th argument had an illegal value




                                  7 Nov 2015                        cunml2(3P)