LAPACK  3.6.1 LAPACK: Linear Algebra PACKage
 lapack_logical LAPACKE_dtb_nancheck ( int matrix_layout, char uplo, char diag, lapack_int n, lapack_int kd, const double * ab, lapack_int ldab )

Definition at line 37 of file lapacke_dtb_nancheck.c.

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

Here is the call graph for this function:

Here is the caller graph for this function: