Lapack-etc.h

Go to the documentation of this file.

#ifndef MATRIX_LAPACK_ETC_H
#define MATRIX_LAPACK_ETC_H
 
/* Copy and paste from WRE : */
 
// before any R headers, or define in PKG_CPPFLAGS
#ifndef USE_FC_LEN_T
# define USE_FC_LEN_T
#endif
#include <Rconfig.h>
 
#ifdef PR18534fixed
# define usePR18534fix
#endif
 
#include <R_ext/Lapack.h>
#ifndef FCONE
# define FCONE
#endif
 
#define ERROR_LAPACK_1(_ROUTINE_, _INFO_) \
do { \
    if ((_INFO_) < 0) \
        error(_("LAPACK routine '%s': argument %d had illegal value"), \
              #_ROUTINE_, -(_INFO_)); \
} while (0)
 
#define ERROR_LAPACK_2(_ROUTINE_, _INFO_, _WARN_, _LETTER_) \
do { \
    ERROR_LAPACK_1(_ROUTINE_, _INFO_); \
    if ((_INFO_) > 0 && (_WARN_) > 0) { \
        if (_WARN_ > 1) \
            error  (_("LAPACK routine '%s': matrix is exactly singular, %s[i,i]=0, i=%d"), \
                    #_ROUTINE_, #_LETTER_, (_INFO_)); \
        else \
            warning(_("LAPACK routine '%s': matrix is exactly singular, %s[i,i]=0, i=%d"), \
                    #_ROUTINE_, #_LETTER_, (_INFO_)); \
    } \
} while (0)
 
#define ERROR_LAPACK_3(_ROUTINE_, _INFO_, _WARN_, _NPROTECT_) \
do { \
    ERROR_LAPACK_1(_ROUTINE_, _INFO_); \
    if ((_INFO_) > 0 && (_WARN_) > 0) { \
        if (_WARN_ > 1) \
            error  (_("LAPACK routine '%s': leading principal minor of order %d is not positive"), \
                    #_ROUTINE_, (_INFO_)); \
        else { \
            warning(_("LAPACK routine '%s': leading principal minor of order %d is not positive"), \
                    #_ROUTINE_, (_INFO_)); \
            UNPROTECT(_NPROTECT_); \
            return ScalarInteger(_INFO_); \
        } \
    } \
} while (0)
 
#define ERROR_LAPACK_4(_ROUTINE_, _INFO_, _RANK_, _WARN_) \
    do { \
        ERROR_LAPACK_1(_ROUTINE_, _INFO_); \
        if ((_INFO_) > 0 && (_WARN_) > 0) { \
            if (_WARN_ > 1) \
                error  (_("LAPACK routine '%s': matrix is rank deficient or not positive definite, the _computed_ rank is %d"), \
                        #_ROUTINE_, (_RANK_)); \
            else \
                warning(_("LAPACK routine '%s': matrix is rank deficient or not positive definite, the _computed_ rank is %d"), \
                        #_ROUTINE_, (_RANK_)); \
        } \
    } while (0)
 
#endif /* MATRIX_LAPACK_ETC_H */