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 */