LAPACK  3.6.1
LAPACK: Linear Algebra PACKage
void LAPACKE_stb_trans ( int  matrix_layout,
char  uplo,
char  diag,
lapack_int  n,
lapack_int  kd,
const float *  in,
lapack_int  ldin,
float *  out,
lapack_int  ldout 
)

Definition at line 40 of file lapacke_stb_trans.c.

44 {
45  lapack_logical colmaj, upper, unit;
46 
47  if( in == NULL || out == NULL ) return;
48 
49  colmaj = ( matrix_layout == LAPACK_COL_MAJOR );
50  upper = LAPACKE_lsame( uplo, 'u' );
51  unit = LAPACKE_lsame( diag, 'u' );
52 
53  if( ( !colmaj && ( matrix_layout != LAPACK_ROW_MAJOR ) ) ||
54  ( !upper && !LAPACKE_lsame( uplo, 'l' ) ) ||
55  ( !unit && !LAPACKE_lsame( diag, 'n' ) ) ) {
56  /* Just exit if any of input parameters are wrong */
57  return;
58  }
59 
60  if( unit ) {
61  /* Unit case, diagonal excluded from transposition */
62  if( colmaj ) {
63  if( upper ) {
64  LAPACKE_sgb_trans( matrix_layout, n-1, n-1, 0, kd-1,
65  &in[ldin], ldin, &out[1], ldout );
66  } else {
67  LAPACKE_sgb_trans( matrix_layout, n-1, n-1, kd-1, 0,
68  &in[1], ldin, &out[ldout], ldout );
69  }
70  } else {
71  if( upper ) {
72  LAPACKE_sgb_trans( matrix_layout, n-1, n-1, 0, kd-1,
73  &in[1], ldin, &out[ldout], ldout );
74  } else {
75  LAPACKE_sgb_trans( matrix_layout, n-1, n-1, kd-1, 0,
76  &in[ldin], ldin, &out[1], ldout );
77  }
78  }
79  } else {
80  /* Non-unit case */
81  if( upper ) {
82  LAPACKE_sgb_trans( matrix_layout, n, n, 0, kd, in, ldin, out,
83  ldout );
84  } else {
85  LAPACKE_sgb_trans( matrix_layout, n, n, kd, 0, in, ldin, out,
86  ldout );
87  }
88  }
89 }
#define lapack_logical
Definition: lapacke.h:51
#define LAPACK_ROW_MAJOR
Definition: lapacke.h:119
lapack_logical LAPACKE_lsame(char ca, char cb)
Definition: lapacke_lsame.c:36
#define LAPACK_COL_MAJOR
Definition: lapacke.h:120
void LAPACKE_sgb_trans(int matrix_layout, lapack_int m, lapack_int n, lapack_int kl, lapack_int ku, const float *in, lapack_int ldin, float *out, lapack_int ldout)

Here is the call graph for this function:

Here is the caller graph for this function: