Lapack interface for TBCI NumLib. More...

#include "tbci/basics.h"
#include "tbci/cplx.h"
#include "tbci/vector.h"
#include "tbci/f_matrix.h"
#include "tbci/f_bandmatrix.h"

Include dependency graph for lapack.h:

Go to the source code of this file.

Functions
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)
TVector< double >	lu_solve_expert (const F_BandMatrix< double > &A, const Vector< double > &B, int equilibrate=1)
	Solution of linear eqution systems, partial pivoting.
TVector< double >	lu_solve (F_BandMatrix< double > &A, const Vector< double > &B, int num_lower, int num_upper)
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)
int	eig (const F_Matrix< double > &A, Vector< double > &EW)
	eigenvalues (and eigenvectors) AEV=EWEV of symmetric double Matrix A
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)
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)
	complex hermitean
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)
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)
	selected eigenvalues -vectors: AEV=EWBEV* of hermitean Matrix A, positive definite Matrix B

Variables
NAMESPACE_TBCI typedef cplx< double >	COMPLEX_DOUBLE

Detailed Description

Lapack interface for TBCI NumLib.

Definition in file lapack.h.

Function Documentation

◆ eig() [1/9]

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 )

selected eigenvalues -vectors: A*EV=EW*B*EV of hermitean Matrix A, positive definite Matrix B

Definition at line 419 of file lapack.cpp.

References doublecomplex, doublecomplex::i, i, imag(), own_ev_(), own_ew_(), doublecomplex::r, real, BVector< T >::resize(), STD__, x, and y.

◆ eig() [2/9]

int eig	(	const F_BandMatrix< double > &	A,
		Vector< double > &	EW,
		F_Matrix< double > &	EV )

Definition at line 296 of file lapack.cpp.

References F_BandMatrix< T >::columns(), dsbev_(), F_BandMatrix< T >::get_fortran_matrix(), F_Matrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), integer, F_BandMatrix< T >::numSub(), and F_BandMatrix< T >::numSuper().

◆ eig() [3/9]

int eig	(	const F_BandMatrix< double > &	A,
		Vector< double > &	EW,
		F_Matrix< double > &	EV,
		double	ew_min,
		double	ew_max )

Definition at line 320 of file lapack.cpp.

References F_BandMatrix< T >::columns(), dsbevx_(), F_Matrix< T >::get_col(), F_BandMatrix< T >::get_fortran_matrix(), F_Matrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), i, integer, F_BandMatrix< T >::numSub(), F_BandMatrix< T >::numSuper(), BVector< T >::resize(), F_Matrix< T >::resize(), and F_Matrix< T >::set_col().

◆ eig() [4/9]

int eig	(	const F_Matrix< COMPLEX_DOUBLE > &	A,
		Vector< COMPLEX_DOUBLE > &	EW )

◆ eig() [5/9]

int eig	(	const F_Matrix< COMPLEX_DOUBLE > &	A,
		Vector< COMPLEX_DOUBLE > &	EW,
		F_Matrix< COMPLEX_DOUBLE > &	EV )

References NAMESPACE_END.

◆ eig() [6/9]

int eig	(	const F_Matrix< cplx< double > > &	A,
		Vector< double > &	EW,
		F_Matrix< cplx< double > > &	EVec )

complex hermitean

Definition at line 539 of file lapack.cpp.

References doublecomplex, BVector< T >::get_fortran_vector(), integer, STD__, and zheev_().

◆ eig() [7/9]

int eig	(	const F_Matrix< double > &	A,
		Vector< double > &	EW )

eigenvalues (and eigenvectors) A*EV=EW*EV of symmetric double Matrix A

Definition at line 263 of file lapack.cpp.

References F_Matrix< T >::columns(), dsyev_(), F_Matrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), and integer.

◆ eig() [8/9]

int eig	(	const F_Matrix< double > &	A,
		Vector< double > &	EW,
		F_Matrix< double > &	EV )

Definition at line 279 of file lapack.cpp.

References F_Matrix< T >::columns(), dsyev_(), F_Matrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), and integer.

◆ eig() [9/9]

int eig	(	F_Matrix< cplx< double > >	A,
		Vector< double > &	EW )

Definition at line 567 of file lapack.cpp.

References doublecomplex, BVector< T >::get_fortran_vector(), integer, STD__, and zheev_().

◆ inv() [1/2]

F_TMatrix< COMPLEX_DOUBLE > inv	(	const F_Matrix< COMPLEX_DOUBLE > &	A,
		int	overwriteA = 0 )

◆ inv() [2/2]

F_TMatrix< double > inv	(	const F_Matrix< double > &	A,
		int	overwriteA = 0 )

Definition at line 181 of file lapack.cpp.

References F_Matrix< T >::columns(), dgesv_(), F_Matrix< T >::get_fortran_matrix(), F_TMatrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), i, integer, F_Matrix< T >::resize(), F_TMatrix< T >::setval(), and STD__.

Referenced by LU_invert(), and LU_invert().

◆ lu_solve() [1/5]

TVector< double > lu_solve	(	const F_BandMatrix< double > &	A,
		const Vector< double > &	B )

Definition at line 74 of file lapack.cpp.

References F_BandMatrix< T >::clear(), F_BandMatrix< T >::columns(), dgbsv_(), F_BandMatrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), i, integer, F_BandMatrix< T >::ldab(), F_BandMatrix< T >::numSub(), F_BandMatrix< T >::numSuper(), and STD__.

◆ lu_solve() [2/5]

F_TMatrix< COMPLEX_DOUBLE > lu_solve	(	const F_Matrix< COMPLEX_DOUBLE > &	A,
		const F_Matrix< COMPLEX_DOUBLE > &	B,
		int	overwriteA = 0 )

◆ lu_solve() [3/5]

F_TMatrix< double > lu_solve	(	const F_Matrix< double > &	A,
		const F_Matrix< double > &	B,
		int	overwriteA = 0 )

Definition at line 156 of file lapack.cpp.

References F_Matrix< T >::columns(), dgesv_(), F_Matrix< T >::get_fortran_matrix(), F_TMatrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), integer, and STD__.

Referenced by expm(), and TMatrix< T >::FGD.

◆ lu_solve() [4/5]

TVector< double > lu_solve	(	const F_Matrix< double > &	A,
		const Vector< double > &	B,
		int	overwriteA = 0 )

Definition at line 129 of file lapack.cpp.

References F_Matrix< T >::columns(), dgesv_(), F_Matrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), integer, and STD__.

◆ lu_solve() [5/5]

TVector< double > lu_solve	(	F_BandMatrix< double > &	A,
		const Vector< double > &	B,
		int	num_lower,
		int	num_upper )

Definition at line 105 of file lapack.cpp.

References F_BandMatrix< T >::columns(), dgbsv_(), F_BandMatrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), integer, F_BandMatrix< T >::ldab(), F_BandMatrix< T >::numSuper(), and STD__.

◆ lu_solve_expert()

TVector< double > lu_solve_expert	(	const F_BandMatrix< double > &	A,
		const Vector< double > &	B,
		int	equilibrate = 1 )

Solution of linear eqution systems, partial pivoting.

Definition at line 24 of file lapack.cpp.

References c, F_BandMatrix< T >::columns(), dgbsvx_(), dlamch_(), fact(), F_BandMatrix< T >::get_fortran_matrix(), BVector< T >::get_fortran_vector(), integer, F_BandMatrix< T >::ldab(), F_BandMatrix< T >::numSub(), F_BandMatrix< T >::numSuper(), TVector< T >::size(), STD__, and x.

Variable Documentation

◆ COMPLEX_DOUBLE

NAMESPACE_TBCI typedef cplx<double> COMPLEX_DOUBLE

Definition at line 18 of file lapack.h.

Functions

Variables