LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
lapack_int LAPACKE_dorbdb_work ( int  matrix_layout,
char  trans,
char  signs,
lapack_int  m,
lapack_int  p,
lapack_int  q,
double *  x11,
lapack_int  ldx11,
double *  x12,
lapack_int  ldx12,
double *  x21,
lapack_int  ldx21,
double *  x22,
lapack_int  ldx22,
double *  theta,
double *  phi,
double *  taup1,
double *  taup2,
double *  tauq1,
double *  tauq2,
double *  work,
lapack_int  lwork 
)

Definition at line 36 of file lapacke_dorbdb_work.c.

44 {
45  lapack_int info = 0;
46  if( matrix_layout == LAPACK_COL_MAJOR ) {
47  /* Call LAPACK function and adjust info */
48  LAPACK_dorbdb( &trans, &signs, &m, &p, &q, x11, &ldx11, x12, &ldx12,
49  x21, &ldx21, x22, &ldx22, theta, phi, taup1, taup2,
50  tauq1, tauq2, work, &lwork, &info );
51  if( info < 0 ) {
52  info = info - 1;
53  }
54  } else if( matrix_layout == LAPACK_ROW_MAJOR ) {
55  lapack_int nrows_x11 = ( LAPACKE_lsame( trans, 'n' ) ? p : q);
56  lapack_int nrows_x12 = ( LAPACKE_lsame( trans, 'n' ) ? p : m-q);
57  lapack_int nrows_x21 = ( LAPACKE_lsame( trans, 'n' ) ? m-p : q);
58  lapack_int nrows_x22 = ( LAPACKE_lsame( trans, 'n' ) ? m-p : m-q);
59  lapack_int ldx11_t = MAX(1,nrows_x11);
60  lapack_int ldx12_t = MAX(1,nrows_x12);
61  lapack_int ldx21_t = MAX(1,nrows_x21);
62  lapack_int ldx22_t = MAX(1,nrows_x22);
63  double* x11_t = NULL;
64  double* x12_t = NULL;
65  double* x21_t = NULL;
66  double* x22_t = NULL;
67  /* Check leading dimension(s) */
68  if( ldx11 < q ) {
69  info = -8;
70  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
71  return info;
72  }
73  if( ldx12 < m-q ) {
74  info = -10;
75  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
76  return info;
77  }
78  if( ldx21 < q ) {
79  info = -12;
80  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
81  return info;
82  }
83  if( ldx22 < m-q ) {
84  info = -14;
85  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
86  return info;
87  }
88  /* Query optimal working array(s) size if requested */
89  if( lwork == -1 ) {
90  LAPACK_dorbdb( &trans, &signs, &m, &p, &q, x11, &ldx11_t, x12,
91  &ldx12_t, x21, &ldx21_t, x22, &ldx22_t, theta, phi,
92  taup1, taup2, tauq1, tauq2, work, &lwork, &info );
93  return (info < 0) ? (info - 1) : info;
94  }
95  /* Allocate memory for temporary array(s) */
96  x11_t = (double*)LAPACKE_malloc( sizeof(double) * ldx11_t * MAX(1,q) );
97  if( x11_t == NULL ) {
98  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
99  goto exit_level_0;
100  }
101  x12_t = (double*)
102  LAPACKE_malloc( sizeof(double) * ldx12_t * MAX(1,m-q) );
103  if( x12_t == NULL ) {
104  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
105  goto exit_level_1;
106  }
107  x21_t = (double*)LAPACKE_malloc( sizeof(double) * ldx21_t * MAX(1,q) );
108  if( x21_t == NULL ) {
109  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
110  goto exit_level_2;
111  }
112  x22_t = (double*)
113  LAPACKE_malloc( sizeof(double) * ldx22_t * MAX(1,m-q) );
114  if( x22_t == NULL ) {
115  info = LAPACK_TRANSPOSE_MEMORY_ERROR;
116  goto exit_level_3;
117  }
118  /* Transpose input matrices */
119  LAPACKE_dge_trans( matrix_layout, nrows_x11, q, x11, ldx11, x11_t,
120  ldx11_t );
121  LAPACKE_dge_trans( matrix_layout, nrows_x12, m-q, x12, ldx12, x12_t,
122  ldx12_t );
123  LAPACKE_dge_trans( matrix_layout, nrows_x21, q, x21, ldx21, x21_t,
124  ldx21_t );
125  LAPACKE_dge_trans( matrix_layout, nrows_x22, m-q, x22, ldx22, x22_t,
126  ldx22_t );
127  /* Call LAPACK function and adjust info */
128  LAPACK_dorbdb( &trans, &signs, &m, &p, &q, x11_t, &ldx11_t, x12_t,
129  &ldx12_t, x21_t, &ldx21_t, x22_t, &ldx22_t, theta, phi,
130  taup1, taup2, tauq1, tauq2, work, &lwork, &info );
131  if( info < 0 ) {
132  info = info - 1;
133  }
134  /* Transpose output matrices */
135  LAPACKE_dge_trans( LAPACK_COL_MAJOR, nrows_x11, q, x11_t, ldx11_t, x11,
136  ldx11 );
137  LAPACKE_dge_trans( LAPACK_COL_MAJOR, nrows_x12, m-q, x12_t, ldx12_t,
138  x12, ldx12 );
139  LAPACKE_dge_trans( LAPACK_COL_MAJOR, nrows_x21, q, x21_t, ldx21_t, x21,
140  ldx21 );
141  LAPACKE_dge_trans( LAPACK_COL_MAJOR, nrows_x22, m-q, x22_t, ldx22_t,
142  x22, ldx22 );
143  /* Release memory and exit */
144  LAPACKE_free( x22_t );
145 exit_level_3:
146  LAPACKE_free( x21_t );
147 exit_level_2:
148  LAPACKE_free( x12_t );
149 exit_level_1:
150  LAPACKE_free( x11_t );
151 exit_level_0:
152  if( info == LAPACK_TRANSPOSE_MEMORY_ERROR ) {
153  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
154  }
155  } else {
156  info = -1;
157  LAPACKE_xerbla( "LAPACKE_dorbdb_work", info );
158  }
159  return info;
160 }
LAPACKE_dge_trans
void LAPACKE_dge_trans(int matrix_layout, lapack_int m, lapack_int n, const double *in, lapack_int ldin, double *out, lapack_int ldout)
Definition: lapacke_dge_trans.c:40
LAPACK_dorbdb
void LAPACK_dorbdb(char *trans, char *signs, lapack_int *m, lapack_int *p, lapack_int *q, double *x11, lapack_int *ldx11, double *x12, lapack_int *ldx12, double *x21, lapack_int *ldx21, double *x22, lapack_int *ldx22, double *theta, double *phi, double *taup1, double *taup2, double *tauq1, double *tauq2, double *work, lapack_int *lwork, lapack_int *info)
LAPACK_ROW_MAJOR
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:119
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

Here is the call graph for this function:

Here is the caller graph for this function: