lapacke_stb_nancheck.c File Reference

#include "lapacke_utils.h"

Go to the source code of this file.

Functions
lapack_logical	LAPACKE_stb_nancheck (int matrix_layout, char uplo, char diag, lapack_int n, lapack_int kd, const float *ab, lapack_int ldab)

Function Documentation

LAPACKE_stb_nancheck()

lapack_logical LAPACKE_stb_nancheck	(	int	matrix_layout,
		char	uplo,
		char	diag,
		lapack_int	n,
		lapack_int	kd,
		const float *	ab,
		lapack_int	ldab )

Definition at line 36 of file lapacke_stb_nancheck.c.

{
    lapack_logical colmaj, upper, unit;
 
    if( ab == NULL ) return (lapack_logical) 0;
 
    colmaj = ( matrix_layout == LAPACK_COL_MAJOR );
    upper  = LAPACKE_lsame( uplo, 'u' );
    unit   = LAPACKE_lsame( diag, 'u' );
 
    if( ( !colmaj && ( matrix_layout != LAPACK_ROW_MAJOR ) ) ||
        ( !upper  && !LAPACKE_lsame( uplo, 'l' ) ) ||
        ( !unit   && !LAPACKE_lsame( diag, 'n' ) ) ) {
        /* Just exit if any of input parameters are wrong */
        return (lapack_logical) 0;
    }
 
    if( unit ) {
        /* Unit case, diagonal should be excluded from the check for NaN. */
        if( colmaj ) {
            if( upper ) {
                return LAPACKE_sgb_nancheck( matrix_layout, n-1, n-1, 0, kd-1,
                                             &ab[ldab], ldab );
            } else {
                return LAPACKE_sgb_nancheck( matrix_layout, n-1, n-1, kd-1, 0,
                                             &ab[1], ldab );
            }
        } else {
            if( upper ) {
                return LAPACKE_sgb_nancheck( matrix_layout, n-1, n-1, 0, kd-1,
                                             &ab[1], ldab );
            } else {
                return LAPACKE_sgb_nancheck( matrix_layout, n-1, n-1, kd-1, 0,
                                             &ab[ldab], ldab );
            }
        }
    } else {
        /* Non-unit case */
        if( upper ) {
            return LAPACKE_sgb_nancheck( matrix_layout, n, n, 0, kd, ab, ldab );
        } else {
            return LAPACKE_sgb_nancheck( matrix_layout, n, n, kd, 0, ab, ldab );
        }
    }
}