LAPACK 3.11.0
LAPACK: Linear Algebra PACKage
Loading...
Searching...
No Matches

◆ LAPACKE_ztrevc_work()

lapack_int LAPACKE_ztrevc_work ( int  matrix_layout,
char  side,
char  howmny,
const lapack_logical select,
lapack_int  n,
lapack_complex_double t,
lapack_int  ldt,
lapack_complex_double vl,
lapack_int  ldvl,
lapack_complex_double vr,
lapack_int  ldvr,
lapack_int  mm,
lapack_int m,
lapack_complex_double work,
double *  rwork 
)

Definition at line 35 of file lapacke_ztrevc_work.c.

42{
43 lapack_int info = 0;
44 if( matrix_layout == LAPACK_COL_MAJOR ) {
45 /* Call LAPACK function and adjust info */
46 LAPACK_ztrevc( &side, &howmny, select, &n, t, &ldt, vl, &ldvl, vr,
47 &ldvr, &mm, m, work, rwork, &info );
48 if( info < 0 ) {
49 info = info - 1;
50 }
51 } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
52 lapack_int ldt_t = MAX(1,n);
53 lapack_int ldvl_t = MAX(1,n);
54 lapack_int ldvr_t = MAX(1,n);
55 lapack_complex_double* t_t = NULL;
56 lapack_complex_double* vl_t = NULL;
57 lapack_complex_double* vr_t = NULL;
58 /* Check leading dimension(s) */
59 if( ldt < n ) {
60 info = -7;
61 LAPACKE_xerbla( "LAPACKE_ztrevc_work", info );
62 return info;
63 }
64 if( ldvl < mm ) {
65 info = -9;
66 LAPACKE_xerbla( "LAPACKE_ztrevc_work", info );
67 return info;
68 }
69 if( ldvr < mm ) {
70 info = -11;
71 LAPACKE_xerbla( "LAPACKE_ztrevc_work", info );
72 return info;
73 }
74 /* Allocate memory for temporary array(s) */
76 LAPACKE_malloc( sizeof(lapack_complex_double) * ldt_t * MAX(1,n) );
77 if( t_t == NULL ) {
79 goto exit_level_0;
80 }
81 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'l' ) ) {
82 vl_t = (lapack_complex_double*)
84 ldvl_t * MAX(1,mm) );
85 if( vl_t == NULL ) {
87 goto exit_level_1;
88 }
89 }
90 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'r' ) ) {
91 vr_t = (lapack_complex_double*)
93 ldvr_t * MAX(1,mm) );
94 if( vr_t == NULL ) {
96 goto exit_level_2;
97 }
98 }
99 /* Transpose input matrices */
100 LAPACKE_zge_trans( matrix_layout, n, n, t, ldt, t_t, ldt_t );
101 if( ( LAPACKE_lsame( side, 'l' ) || LAPACKE_lsame( side, 'b' ) ) &&
102 LAPACKE_lsame( howmny, 'b' ) ) {
103 LAPACKE_zge_trans( matrix_layout, n, mm, vl, ldvl, vl_t, ldvl_t );
104 }
105 if( ( LAPACKE_lsame( side, 'r' ) || LAPACKE_lsame( side, 'b' ) ) &&
106 LAPACKE_lsame( howmny, 'b' ) ) {
107 LAPACKE_zge_trans( matrix_layout, n, mm, vr, ldvr, vr_t, ldvr_t );
108 }
109 /* Call LAPACK function and adjust info */
110 LAPACK_ztrevc( &side, &howmny, select, &n, t_t, &ldt_t, vl_t, &ldvl_t,
111 vr_t, &ldvr_t, &mm, m, work, rwork, &info );
112 if( info < 0 ) {
113 info = info - 1;
114 }
115 /* Transpose output matrices */
116 LAPACKE_zge_trans( LAPACK_COL_MAJOR, n, n, t_t, ldt_t, t, ldt );
117 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'l' ) ) {
118 LAPACKE_zge_trans( LAPACK_COL_MAJOR, n, mm, vl_t, ldvl_t, vl,
119 ldvl );
120 }
121 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'r' ) ) {
122 LAPACKE_zge_trans( LAPACK_COL_MAJOR, n, mm, vr_t, ldvr_t, vr,
123 ldvr );
124 }
125 /* Release memory and exit */
126 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'r' ) ) {
127 LAPACKE_free( vr_t );
128 }
129exit_level_2:
130 if( LAPACKE_lsame( side, 'b' ) || LAPACKE_lsame( side, 'l' ) ) {
131 LAPACKE_free( vl_t );
132 }
133exit_level_1:
134 LAPACKE_free( t_t );
135exit_level_0:
136 if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
137 LAPACKE_xerbla( "LAPACKE_ztrevc_work", info );
138 }
139 } else {
140 info = -1;
141 LAPACKE_xerbla( "LAPACKE_ztrevc_work", info );
142 }
143 return info;
144}
#define lapack_int
Definition: lapack.h:87
#define lapack_complex_double
Definition: lapack.h:64
#define LAPACK_ztrevc(...)
Definition: lapack.h:21500
#define LAPACK_COL_MAJOR
Definition: lapacke.h:53
#define LAPACKE_free(p)
Definition: lapacke.h:46
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:52
#define LAPACKE_malloc(size)
Definition: lapacke.h:43
#define LAPACK_TRANSPOSE_MEMORY_ERROR
Definition: lapacke.h:56
lapack_logical LAPACKE_lsame(char ca, char cb)
Definition: lapacke_lsame.c:35
void LAPACKE_xerbla(const char *name, lapack_int info)
void LAPACKE_zge_trans(int matrix_layout, lapack_int m, lapack_int n, const lapack_complex_double *in, lapack_int ldin, lapack_complex_double *out, lapack_int ldout)
#define MAX(x, y)
Definition: lapacke_utils.h:46
Here is the call graph for this function:
Here is the caller graph for this function: