TBCI Numerical high perf. C++ Library 2.8.0
bdmatsolverlib.cc
Go to the documentation of this file.
1
6/* $Id: bdmatsolverlib.cc,v 1.5.2.4 2019/05/28 11:13:02 garloff Exp $ */
7
8#ifndef NUM
9# define NUM double
10#endif
11
12#include "tbci/band_matrix.h"
13#ifdef PRAGMA_I
14# pragma implementation "bd_lu_solver.h"
15#endif
16#include "tbci/solver/bd_lu_solver.h"
17#ifdef PRAGMA_I
18# pragma implementation "svd_solver.h"
19#endif
20#include "tbci/solver/svd_solver.h"
21
22//NAMESPACE_TBCI
23
24#define T NUM
25#include "bd_lu_solver_inst.h"
26
27#define MAT TBCI__ BdMatrix<NUM >
28#define VEC TBCI__ Vector<NUM >
29#include "svd_solver_inst.h"
30
31#if 0
32#define TEMPLATE template
33
34// explicit instatiation of all those things ...
35
36//bd_lu_solver.h <NUM>
37TEMPLATE int lu_decomp (BdMatrix<NUM>&);
38TEMPLATE TVector<NUM> LU_fwd_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
39TEMPLATE TVector<NUM> LU_bkw_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
40TEMPLATE TVector<NUM> LU_solve (const BdMatrix<NUM>&, const Vector<NUM>&);
41TEMPLATE TVector<NUM> lu_solve (BdMatrix<NUM>&, const Vector<NUM>&);
42TEMPLATE TMatrix<NUM> LU_solve (const BdMatrix<NUM>&, const Matrix<NUM>&);
43TEMPLATE TMatrix<NUM> lu_solve (BdMatrix<NUM>&, const Matrix<NUM>&);
44TEMPLATE NUM LU_det (const BdMatrix<NUM>&);
45TEMPLATE NUM lu_det (BdMatrix<NUM>&);
46TEMPLATE TMatrix<NUM> LU_invert (const BdMatrix<NUM>&);
47TEMPLATE TMatrix<NUM> lu_invert (BdMatrix<NUM>&);
48#endif /* 0 */
49
50//NAMESPACE_END
LU_bkw_subst
TVector< T > LU_bkw_subst(const BdMatrix< T > &lu, const Vector< T > &y)
Definition bd_lu_solver.h:224
lu_det
T lu_det(BdMatrix< T > &mat)
calculates the determinant of a BdMatrix by doing an LU decomposition
Definition bd_lu_solver.h:329
LU_solve
TVector< T > LU_solve(const BdMatrix< T > &lu, const Vector< T > &b)
Solve the equation Ax = b where A IS already LU decomposed.
Definition bd_lu_solver.h:246
LU_det
T LU_det(const BdMatrix< T > &lu)
calculates determinant of an already LU decomposed BdMatrix
Definition bd_lu_solver.h:317
LU_invert
TMatrix< T > LU_invert(const BdMatrix< T > &lu)
return the inverse Matrix for an already LU decomposed BdMatrix
Definition bd_lu_solver.h:338
LU_fwd_subst
TVector< T > LU_fwd_subst(const BdMatrix< T > &lu, const Vector< T > &y)
Definition bd_lu_solver.h:207
lu_invert
TMatrix< T > lu_invert(BdMatrix< T > &mat)
return the inverse Matrix for a BdMatrix by doing an LU decomposition
Definition bd_lu_solver.h:355
NUM
#define NUM
Definition bdmatlib.cc:9
BdMatrix
The class BdMatrix is an implementation to store and do operations on sparse Matrices with a band str...
Definition band_matrix.h:104
TVector
Temporary Base Class Idiom: Class TVector is used for temporary variables.
Definition vector.h:73
TEMPLATE
#define TEMPLATE
Definition cplxinst.cc:22
lu_solve
F_TMatrix< double > lu_solve(const F_Matrix< double > &A, const F_Matrix< double > &B, int overwriteA=0)
Definition lapack.cpp:156
lu_decomp
int lu_decomp(Matrix< T > &) HOT
LU decomposes the TBCI::Matrix mat.
Definition lu_solver.h:94
Generated by doxygen 1.16.1