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

Definition at line 36 of file lapacke_zhbgvx_work.c.

47 {
48  lapack_int info = 0;
49  if( matrix_layout == LAPACK_COL_MAJOR ) {
50  /* Call LAPACK function and adjust info */
51  LAPACK_zhbgvx( &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_double* ab_t = NULL;
63  lapack_complex_double* bb_t = NULL;
64  lapack_complex_double* q_t = NULL;
65  lapack_complex_double* z_t = NULL;
66  /* Check leading dimension(s) */
67  if( ldab < n ) {
68  info = -9;
69  LAPACKE_xerbla( "LAPACKE_zhbgvx_work", info );
70  return info;
71  }
72  if( ldbb < n ) {
73  info = -11;
74  LAPACKE_xerbla( "LAPACKE_zhbgvx_work", info );
75  return info;
76  }
77  if( ldq < n ) {
78  info = -13;
79  LAPACKE_xerbla( "LAPACKE_zhbgvx_work", info );
80  return info;
81  }
82  if( ldz < n ) {
83  info = -22;
84  LAPACKE_xerbla( "LAPACKE_zhbgvx_work", info );
85  return info;
86  }
87  /* Allocate memory for temporary array(s) */
88  ab_t = (lapack_complex_double*)
89  LAPACKE_malloc( sizeof(lapack_complex_double) * 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_double*)
95  LAPACKE_malloc( sizeof(lapack_complex_double) * 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_double*)
102  LAPACKE_malloc( sizeof(lapack_complex_double) *
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_double*)
111  LAPACKE_malloc( sizeof(lapack_complex_double) *
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_zhb_trans( matrix_layout, uplo, n, ka, ab, ldab, ab_t, ldab_t );
120  LAPACKE_zhb_trans( matrix_layout, uplo, n, kb, bb, ldbb, bb_t, ldbb_t );
121  /* Call LAPACK function and adjust info */
122  LAPACK_zhbgvx( &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_zhb_trans( LAPACK_COL_MAJOR, uplo, n, ka, ab_t, ldab_t, ab,
130  ldab );
131  LAPACKE_zhb_trans( LAPACK_COL_MAJOR, uplo, n, kb, bb_t, ldbb_t, bb,
132  ldbb );
133  if( LAPACKE_lsame( jobz, 'v' ) ) {
134  LAPACKE_zge_trans( LAPACK_COL_MAJOR, n, n, q_t, ldq_t, q, ldq );
135  }
136  if( LAPACKE_lsame( jobz, 'v' ) ) {
137  LAPACKE_zge_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_zhbgvx_work", info );
154  }
155  } else {
156  info = -1;
157  LAPACKE_xerbla( "LAPACKE_zhbgvx_work", info );
158  }
159  return info;
160 }
LAPACK_ROW_MAJOR
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:119
LAPACK_zhbgvx
void LAPACK_zhbgvx(char *jobz, char *range, char *uplo, lapack_int *n, lapack_int *ka, lapack_int *kb, lapack_complex_double *ab, lapack_int *ldab, lapack_complex_double *bb, lapack_int *ldbb, lapack_complex_double *q, lapack_int *ldq, double *vl, double *vu, lapack_int *il, lapack_int *iu, double *abstol, lapack_int *m, double *w, lapack_complex_double *z, lapack_int *ldz, lapack_complex_double *work, double *rwork, lapack_int *iwork, lapack_int *ifail, lapack_int *info)
lapack_complex_double
#define lapack_complex_double
Definition: lapacke.h:90
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
LAPACKE_zhb_trans
void LAPACKE_zhb_trans(int matrix_layout, char uplo, lapack_int n, lapack_int kd, const lapack_complex_double *in, lapack_int ldin, lapack_complex_double *out, lapack_int ldout)
Definition: lapacke_zhb_trans.c:40
LAPACK_TRANSPOSE_MEMORY_ERROR
#define LAPACK_TRANSPOSE_MEMORY_ERROR
Definition: lapacke.h:123
LAPACKE_zge_trans
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)
Definition: lapacke_zge_trans.c:40

Here is the call graph for this function:

Here is the caller graph for this function: