mbs_lapack.h

Go to the documentation of this file.

 
#ifndef MBS_LAPACK_H_
#define MBS_LAPACK_H_
#include <stdlib.h>
#include "mbs_message.h"
 
#ifdef LAPACK_FOUND
 
// necessary for compatibility with VS
#if __STDC_VERSION__ < 1999901L
    struct complex_float { float x, y; };
    struct complex_double { double x, y; };
    #define lapack_complex_float  struct complex_float
    #define lapack_complex_double  struct complex_double
#endif
// include lapacke header to get lapack function declaration
#include "lapacke.h"
 
#else
 
#define LPK_ERR ">> LAPACK>> "
#define LPK_PRE "         >> "
 
// Type conversion.
#define lapack_int int
#define lapack_logical int
 
#define LAPACKE_malloc( size ) malloc( size )
#define LAPACKE_free( p ) free( p )
 
typedef lapack_logical (*LAPACK_D_SELECT2) ( const double*, const double* );
 
// Declare fake functions to catch error (alphabetical ordering)
static inline void LAPACK_dgees( char* jobvs, char* sort, LAPACK_D_SELECT2 select,
                         lapack_int* n, double* a, lapack_int* lda, lapack_int* sdim,
                         double* wr, double* wi, double* vs, lapack_int* ldvs,
                         double* work, lapack_int* lwork, lapack_logical* bwork,
                         lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dgeev( char* jobvl, char* jobvr, lapack_int* n, double* a,
                         lapack_int* lda, double* wr, double* wi, double* vl,
                         lapack_int* ldvl, double* vr, lapack_int* ldvr, double* work,
                         lapack_int* lwork, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dgelss( lapack_int* m, lapack_int* n, lapack_int* nrhs, double* a,
                          lapack_int* lda, double* b, lapack_int* ldb, double* s,
                          double* rcond, lapack_int* rank, double* work,
                          lapack_int* lwork, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dgesvd( char* jobu, char* jobvt, lapack_int* m, lapack_int* n,
                          double* a, lapack_int* lda, double* s, double* u,
                          lapack_int* ldu, double* vt, lapack_int* ldvt, double* work,
                          lapack_int* lwork, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
 
}
 
static inline void LAPACK_dgetrf(lapack_int* m, lapack_int* n, double* a, lapack_int* lda,
                   lapack_int* ipiv, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dgetri( lapack_int* n, double* a, lapack_int* lda,
                   const lapack_int* ipiv, double* work, lapack_int* lwork,
                   lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dpotrf( char* uplo, lapack_int* n, double* a, lapack_int* lda,
                          lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dpotrs( char* uplo, lapack_int* n, lapack_int* nrhs,
                          const double* a, lapack_int* lda, double* b,
                          lapack_int* ldb, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
 
static inline void LAPACK_dtrsen( char* job, char* compq, const lapack_logical* select,
                   lapack_int* n, double* t, lapack_int* ldt, double* q,
                   lapack_int* ldq, double* wr, double* wi, lapack_int* m,
                   double* s, double* sep, double* work, lapack_int* lwork,
                   lapack_int* iwork, lapack_int* liwork, lapack_int *info )
{
    *info = -1;
    mbs_msg(LPK_ERR"%s()\n"
            LPK_PRE"Code compiled without LAPACK libraries.\n"
            LPK_PRE"To use this module, compile MBsysC with LAPACK libraries.\n",
            __func__);
}
#endif
#endif