lapack_lib.h

Go to the documentation of this file.

#ifndef TBCI_LAPACK_LIB_H
#define TBCI_LAPACK_LIB_H

#ifdef __cplusplus
extern "C" {
#endif
#include "tbci/lapack/f2c.h"

#if 0   
typedef /*long*/ int integer;
typedef char *address;
typedef short int shortint;
/*typedef float real;*/
/*#define real float;*/
typedef double doublereal;
/*typedef struct { real r, i; } complex;*/
typedef struct { doublereal r, i; } doublecomplex;
typedef /*long*/ int logical;
typedef short int shortlogical;
typedef char logical1;
typedef char integer1;
#endif

#ifdef _DOXYGEN
# define integer /*long*/ int
# define doublereal double
# define doublecomplex __complex__ double
# define real float
# define complex __complex__ float
#endif

/* FIXME: should this be dlamch_(char) ? */
double dlamch_(const char*);


// ********************** real linear Systems *********************************************
/* Subroutine */
int dgesv_(integer *n, integer *nrhs, doublereal *a, integer *lda,
           integer *ipiv, doublereal *b, integer *ldb, integer *info);

// Band
/* Subroutine */
int dgbsv_(integer *n, integer *kl, integer *ku, integer *nrhs,
           doublereal *ab, integer *ldab, integer *ipiv, doublereal *b,
           integer *ldb, integer *info);
/* Subroutine */
/* FIXME: Should this be char fact, char trans ? */
int dgbsvx_(const char *fact, const char *trans, integer *n, integer *kl,
            integer *ku, integer *nrhs, doublereal *ab, integer *ldab,
            doublereal *afb, integer *ldafb, integer *ipiv, char *equed,
            doublereal *r, doublereal *c, doublereal *b, integer *ldb,
            doublereal *x, integer *ldx,
            doublereal *rcond, doublereal *ferr, doublereal *berr,
            doublereal *work, integer *iwork, integer *info);
/* Subroutine */
int dgbsvx2_(const char *fact, const char *trans, integer *n, integer *kl,
             integer *ku, integer *nrhs, doublereal *ab, integer *ldab,
             doublereal *afb, integer *ldafb, integer *ipiv, char *equed,
             doublereal *r, doublereal *c, doublereal *b, integer *ldb,
             doublereal *x, integer *ldx,
             doublereal *rcond, doublereal *ferr, doublereal *berr,
             doublereal *work, integer *iwork, integer *info);


// ********************** complex linear Systems ******************************************
/* Subroutine */
int zgesv_(integer *n, integer *nrhs, doublecomplex *a,
           integer *lda, integer *ipiv, doublecomplex *b,
           integer *ldb, integer *info);


// ********************** real Eigenvalues ************************************************
int dsyev_(const char *jobz, const char *uplo, integer *n, doublereal *a,
           integer *lda, doublereal *w, doublereal *work, integer *lwork,
           integer *info);

int dspev_(const char *jobz, const char *uplo, integer *n, doublereal *ap,
           doublereal *w, doublereal *z, integer *ldz, doublereal *work,
           integer *info);

int dgeev_(const char *jobvl, const char *jobvr, integer *n, doublereal *a,
           integer *lda, doublereal *wr, doublereal *wi, doublereal *vl,
           integer *ldvl, doublereal *vr, integer *ldvr, doublereal *work,
           integer *lwork, integer *info);

//symmetric band matrix
/* Subroutine */
int dsbev_(const char *jobz, const char *uplo, integer *n, integer *kd,
           doublereal *ab, integer *ldab, doublereal *w, doublereal *z,
           integer *ldz, doublereal *work, integer *info);

/* Subroutine */
int dsbevx_(const char *jobz, const char *range, const char *uplo, integer *n,
            integer *kd, doublereal *ab, integer *ldab, doublereal *q,
            integer *ldq, doublereal *vl, doublereal *vu, integer *il,
            integer *iu, doublereal *abstol, integer *m, doublereal *w,
            doublereal *z, integer *ldz, doublereal *work, integer *iwork,
            integer *ifail, integer *info);

// ********************** complex Eigenvalues ********************************************
int zgeev_(const char *jobvl, const char *jobvr, integer *n,
           doublecomplex *a, integer *lda, doublecomplex *w, doublecomplex *vl,
           integer *ldvl, doublecomplex *vr, integer *ldvr, doublecomplex *work,
           integer *lwork, doublereal *rwork, integer *info);


int zheev_(const char *jobz, const char *uplo, integer *n, doublecomplex *a,
           integer *lda, doublereal *w, doublecomplex *work, integer *lwork,
           doublereal *rwork, integer *info);


int zgehrd_(integer *n, integer *ilo, integer *ihi, doublecomplex *a,
            integer *lda, doublecomplex *tau, doublecomplex *work,
            integer *lwork, integer *info);

int zunghr_(integer *n, integer *ilo, integer *ihi, doublecomplex *a,
            integer *lda, doublecomplex *tau, doublecomplex *work,
            integer *lwork, integer *info);

int zhseqr_(const char *job, const char *compz, integer *n, integer *ilo, integer *ihi,
            doublecomplex *h, integer *ldh, doublecomplex *w,
            doublecomplex *z, integer *ldz, doublecomplex *work,
            integer *lwork, integer *info);



// Eigene Fortran-Routinen
long own_ew_(long *jobz, long *uplo, long *n, long *m, doublecomplex *ab, long *ldab,
             doublecomplex* bb, long *ldbb, double *vl, double *vu, double *w,
             doublecomplex *q, long *ldq, doublecomplex *work, double *rwork,
             long *iwork);
long own_ev_(long *n, long *m, double *w, doublecomplex *z, long *ldz, doublecomplex *q,
             long *ldq, doublecomplex *work, double *rwork, long *iwork, long *ifail);

#ifdef __cplusplus
}
#endif

#endif