lapack.h

Go to the documentation of this file.

/* $Id: lapack.h,v 1.6.2.3 2019/05/28 11:13:02 garloff Exp $ */
 
#ifndef TBCI_LAPACK_H
#define TBCI_LAPACK_H
 
#include "tbci/basics.h" 
#include "tbci/cplx.h" 
#include "tbci/vector.h"
#include "tbci/f_matrix.h"
#include "tbci/f_bandmatrix.h"
 
NAMESPACE_TBCI
 
typedef cplx<double> COMPLEX_DOUBLE;
 
 
//************  Linear Equation systems  ***************************************************
// double
F_TMatrix<double> inv(const F_Matrix<double>& A, int overwriteA=0);
F_TMatrix<double> lu_solve(const F_Matrix<double>& A, const F_Matrix<double>& B, int overwriteA=0);
TVector<double> lu_solve(const F_Matrix<double>& A, const Vector<double>& B, int overwriteA=0);
TVector<double> lu_solve(const F_BandMatrix<double>& A, const Vector<double>& B); //Kopiert A
TVector<double> lu_solve_expert(const F_BandMatrix<double>& A, const Vector<double>& B, int equilibrate=1); //A,B ueberschreiben
TVector<double> lu_solve(F_BandMatrix<double>& A, const Vector<double>& B, int num_lower, int num_upper); //Ueberschreibt A
 
// complex
F_TMatrix<COMPLEX_DOUBLE> inv(const F_Matrix<COMPLEX_DOUBLE>& A, int overwriteA=0);
F_TMatrix<COMPLEX_DOUBLE> lu_solve(const F_Matrix<COMPLEX_DOUBLE>& A, const F_Matrix<COMPLEX_DOUBLE>& B, int overwriteA=0);
 
 
//************  Standard Eigenvalue problems  **********************************************
// eigenvalues (and eigenvectors) A*EV=EW*EV of symmetric double Matrix A
int eig(const F_Matrix<double>& A, Vector<double>& EW);
int eig(const F_Matrix<double>& A, Vector<double>& EW, F_Matrix<double>& EV);
int eig(const F_BandMatrix<double>& A, Vector<double>& EW, F_Matrix<double>& EV);
int eig(const F_BandMatrix<double>& A, Vector<double>& EW, F_Matrix<double>& EV, 
        double ew_min, double ew_max);
 
// eigenvalues (and eigenvectors) A*EV=EW*EV of hermitean complex Matrix A
int eig(F_Matrix<cplx<double> > A, Vector<double>& EW);
int eig(const F_Matrix<cplx<double> >& A, Vector<double>& EW, F_Matrix<cplx<double> >& EVec);
 
// eigenvalues (and eigenvectors) A*EV=EW*EV of general complex Matrix A
int eig(const F_Matrix<COMPLEX_DOUBLE>& A, Vector<COMPLEX_DOUBLE>& EW);
int eig(const F_Matrix<COMPLEX_DOUBLE>& A, Vector<COMPLEX_DOUBLE>& EW, F_Matrix<COMPLEX_DOUBLE>& EV);
 
 
//************  Generalized Eigenvalue problems  ******************************************
int eig(const F_BandMatrix<cplx<double> >& A, const F_BandMatrix<cplx<double> >& B,
        double EW_min, double EW_max,
        Vector<double>& Eigenwerte, F_Matrix<cplx<double> >& Eigenvectoren);
 
NAMESPACE_END
#endif /* TBCI_LAPACK_H */