LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
lapack_int LAPACKE_chbgvx_work ( int  matrix_layout,
char  jobz,
char  range,
char  uplo,
lapack_int  n,
lapack_int  ka,
lapack_int  kb,
lapack_complex_float *  ab,
lapack_int  ldab,
lapack_complex_float *  bb,
lapack_int  ldbb,
lapack_complex_float *  q,
lapack_int  ldq,
float  vl,
float  vu,
lapack_int  il,
lapack_int  iu,
float  abstol,
lapack_int *  m,
float *  w,
lapack_complex_float *  z,
lapack_int  ldz,
lapack_complex_float *  work,
float *  rwork,
lapack_int *  iwork,
lapack_int *  ifail 
)

Definition at line 36 of file lapacke_chbgvx_work.c.

47 {
48  lapack_int info = 0;
49  if( matrix_layout == LAPACK_COL_MAJOR ) {
50  /* Call LAPACK function and adjust info */
51  LAPACK_chbgvx( &jobz, &range, &uplo, &n, &ka, &kb, ab, &ldab, bb, &ldbb,
52  q, &ldq, &vl, &vu, &il, &iu, &abstol, m, w, z, &ldz,
53  work, rwork, iwork, ifail, &info );
54  if( info < 0 ) {
55  info = info - 1;
56  }
57  } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
58  lapack_int ldab_t = MAX(1,ka+1);
59  lapack_int ldbb_t = MAX(1,kb+1);
60  lapack_int ldq_t = MAX(1,n);
61  lapack_int ldz_t = MAX(1,n);
62  lapack_complex_float* ab_t = NULL;
63  lapack_complex_float* bb_t = NULL;
64  lapack_complex_float* q_t = NULL;
65  lapack_complex_float* z_t = NULL;
66  /* Check leading dimension(s) */
67  if( ldab < n ) {
68  info = -9;
69  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
70  return info;
71  }
72  if( ldbb < n ) {
73  info = -11;
74  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
75  return info;
76  }
77  if( ldq < n ) {
78  info = -13;
79  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
80  return info;
81  }
82  if( ldz < n ) {
83  info = -22;
84  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
85  return info;
86  }
87  /* Allocate memory for temporary array(s) */
88  ab_t = (lapack_complex_float*)
89  LAPACKE_malloc( sizeof(lapack_complex_float) * ldab_t * MAX(1,n) );
90  if( ab_t == NULL ) {
91  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
92  goto exit_level_0;
93  }
94  bb_t = (lapack_complex_float*)
95  LAPACKE_malloc( sizeof(lapack_complex_float) * ldbb_t * MAX(1,n) );
96  if( bb_t == NULL ) {
97  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
98  goto exit_level_1;
99  }
100  if( LAPACKE_lsame( jobz, 'v' ) ) {
101  q_t = (lapack_complex_float*)
102  LAPACKE_malloc( sizeof(lapack_complex_float) *
103  ldq_t * MAX(1,n) );
104  if( q_t == NULL ) {
105  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
106  goto exit_level_2;
107  }
108  }
109  if( LAPACKE_lsame( jobz, 'v' ) ) {
110  z_t = (lapack_complex_float*)
111  LAPACKE_malloc( sizeof(lapack_complex_float) *
112  ldz_t * MAX(1,n) );
113  if( z_t == NULL ) {
114  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
115  goto exit_level_3;
116  }
117  }
118  /* Transpose input matrices */
119  LAPACKE_chb_trans( matrix_layout, uplo, n, ka, ab, ldab, ab_t, ldab_t );
120  LAPACKE_chb_trans( matrix_layout, uplo, n, kb, bb, ldbb, bb_t, ldbb_t );
121  /* Call LAPACK function and adjust info */
122  LAPACK_chbgvx( &jobz, &range, &uplo, &n, &ka, &kb, ab_t, &ldab_t, bb_t,
123  &ldbb_t, q_t, &ldq_t, &vl, &vu, &il, &iu, &abstol, m, w,
124  z_t, &ldz_t, work, rwork, iwork, ifail, &info );
125  if( info < 0 ) {
126  info = info - 1;
127  }
128  /* Transpose output matrices */
129  LAPACKE_chb_trans( LAPACK_COL_MAJOR, uplo, n, ka, ab_t, ldab_t, ab,
130  ldab );
131  LAPACKE_chb_trans( LAPACK_COL_MAJOR, uplo, n, kb, bb_t, ldbb_t, bb,
132  ldbb );
133  if( LAPACKE_lsame( jobz, 'v' ) ) {
134  LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );
135  }
136  if( LAPACKE_lsame( jobz, 'v' ) ) {
137  LAPACKE_cge_trans( LAPACK_COL_MAJOR, n, n, z_t, ldz_t, z, ldz );
138  }
139  /* Release memory and exit */
140  if( LAPACKE_lsame( jobz, 'v' ) ) {
141  LAPACKE_free( z_t );
142  }
143 exit_level_3:
144  if( LAPACKE_lsame( jobz, 'v' ) ) {
145  LAPACKE_free( q_t );
146  }
147 exit_level_2:
148  LAPACKE_free( bb_t );
149 exit_level_1:
150  LAPACKE_free( ab_t );
151 exit_level_0:
152  if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
153  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
154  }
155  } else {
156  info = -1;
157  LAPACKE_xerbla( "LAPACKE_chbgvx_work", info );
158  }
159  return info;
160 }
LAPACKE_cge_trans
void LAPACKE_cge_trans(int matrix_layout, lapack_int m, lapack_int n, const lapack_complex_float *in, lapack_int ldin, lapack_complex_float *out, lapack_int ldout)
Definition: lapacke_cge_trans.c:40
LAPACK_chbgvx
void LAPACK_chbgvx(char *jobz, char *range, char *uplo, lapack_int *n, lapack_int *ka, lapack_int *kb, lapack_complex_float *ab, lapack_int *ldab, lapack_complex_float *bb, lapack_int *ldbb, lapack_complex_float *q, lapack_int *ldq, float *vl, float *vu, lapack_int *il, lapack_int *iu, float *abstol, lapack_int *m, float *w, lapack_complex_float *z, lapack_int *ldz, lapack_complex_float *work, float *rwork, lapack_int *iwork, lapack_int *ifail, lapack_int *info)
LAPACK_ROW_MAJOR
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:119
lapack_complex_float
#define lapack_complex_float
Definition: lapacke.h:74
MAX
#define MAX(x, y)
Definition: lapacke_utils.h:47
LAPACKE_free
#define LAPACKE_free(p)
Definition: lapacke.h:113
LAPACKE_malloc
#define LAPACKE_malloc(size)
Definition: lapacke.h:110
LAPACKE_lsame
lapack_logical LAPACKE_lsame(char ca, char cb)
Definition: lapacke_lsame.c:36
LAPACK_COL_MAJOR
#define LAPACK_COL_MAJOR
Definition: lapacke.h:120
LAPACKE_xerbla
void LAPACKE_xerbla(const char *name, lapack_int info)
Definition: lapacke_xerbla.c:37
lapack_int
#define lapack_int
Definition: lapacke.h:47
LAPACK_TRANSPOSE_MEMORY_ERROR
#define LAPACK_TRANSPOSE_MEMORY_ERROR
Definition: lapacke.h:123
LAPACKE_chb_trans
void LAPACKE_chb_trans(int matrix_layout, char uplo, lapack_int n, lapack_int kd, const lapack_complex_float *in, lapack_int ldin, lapack_complex_float *out, lapack_int ldout)
Definition: lapacke_chb_trans.c:40

Here is the call graph for this function:

Here is the caller graph for this function: