TBCI Numerical high perf. C++ Library: bdmatsolverlib.cc Source File

00001 
00006 /* $Id: bdmatsolverlib.cc,v 1.5.2.3 2001/11/16 01:21:46 garloff Exp $ */
00007 
00008 #ifndef NUM
00009 # define NUM double
00010 #endif
00011 
00012 #include "band_matrix.h"
00013 #ifdef PRAGMA_I
00014 # pragma implementation "bd_lu_solver.h"
00015 #endif
00016 #include "solver/bd_lu_solver.h"
00017 #ifdef PRAGMA_I
00018 # pragma implementation "svd_solver.h"
00019 #endif
00020 #include "solver/svd_solver.h"
00021 
00022 //NAMESPACE_TBCI
00023 
00024 #define T NUM
00025 #include "bd_lu_solver_inst.h"
00026 
00027 #define MAT TBCI__ BdMatrix<NUM > 
00028 #define VEC TBCI__ Vector<NUM >
00029 #include "svd_solver_inst.h"
00030 
00031 #if 0
00032 #define TEMPLATE template
00033 
00034 // explicit instatiation of all those things ...
00035 
00036 //bd_lu_solver.h <NUM>
00037 TEMPLATE int lu_decomp (BdMatrix<NUM>&);
00038 TEMPLATE TVector<NUM> LU_fwd_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
00039 TEMPLATE TVector<NUM> LU_bkw_subst (const BdMatrix<NUM>&, const Vector<NUM>&);
00040 TEMPLATE TVector<NUM> LU_solve (const BdMatrix<NUM>&, const Vector<NUM>&);
00041 TEMPLATE TVector<NUM> lu_solve (BdMatrix<NUM>&, const Vector<NUM>&);
00042 TEMPLATE TMatrix<NUM> LU_solve (const BdMatrix<NUM>&, const Matrix<NUM>&);
00043 TEMPLATE TMatrix<NUM> lu_solve (BdMatrix<NUM>&, const Matrix<NUM>&);
00044 TEMPLATE NUM LU_det (const BdMatrix<NUM>&);
00045 TEMPLATE NUM lu_det (BdMatrix<NUM>&);
00046 TEMPLATE TMatrix<NUM> LU_invert (const BdMatrix<NUM>&);
00047 TEMPLATE TMatrix<NUM> lu_invert (BdMatrix<NUM>&);
00048 #endif /* 0 */
00049 
00050 //NAMESPACE_END