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>
37 TEMPLATE int lu_decomp (BdMatrix<NUM>&);
38 TEMPLATE TVector<NUM> LU_fwd_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
39 TEMPLATE TVector<NUM> LU_bkw_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
40 TEMPLATE TVector<NUM> LU_solve (const BdMatrix<NUM>&, const Vector<NUM>&);
41 TEMPLATE TVector<NUM> lu_solve (BdMatrix<NUM>&, const Vector<NUM>&);
42 TEMPLATE TMatrix<NUM> LU_solve (const BdMatrix<NUM>&, const Matrix<NUM>&);
43 TEMPLATE TMatrix<NUM> lu_solve (BdMatrix<NUM>&, const Matrix<NUM>&);
44 TEMPLATE NUM LU_det (const BdMatrix<NUM>&);
45 TEMPLATE NUM lu_det (BdMatrix<NUM>&);
46 TEMPLATE TMatrix<NUM> LU_invert (const BdMatrix<NUM>&);
47 TEMPLATE TMatrix<NUM> lu_invert (BdMatrix<NUM>&);
48 #endif /* 0 */
49 
50 //NAMESPACE_END
LU_fwd_subst
TVector< T > LU_fwd_subst(const BdMatrix< T > &lu, const Vector< T > &y)
Definition: bd_lu_solver.h:207
BdMatrix
The class BdMatrix is an implementation to store and do operations on sparse Matrices with a band str...
Definition: band_matrix.h:103
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_bkw_subst
TVector< T > LU_bkw_subst(const BdMatrix< T > &lu, const Vector< T > &y)
Definition: bd_lu_solver.h:224
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_decomp
int lu_decomp(Matrix< T > &) HOT
LU decomposes the TBCI::Matrix mat.
Definition: lu_solver.h:94
LU_det
T LU_det(const BdMatrix< T > &lu)
calculates determinant of an already LU decomposed BdMatrix
Definition: bd_lu_solver.h:317
TEMPLATE
#define TEMPLATE
Definition: cplxinst.cc:22
Matrix
Definition: bvector.h:49
TVector
Temporary Base Class Idiom: Class TVector is used for temporary variables.
Definition: bvector.h:52
Vector
Definition: bvector.h:54
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
NUM
#define NUM
Definition: bdmatsolverlib.cc:9
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
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
Generated by   doxygen 1.8.5