Changeset 616 for trunk/c++_tools/utility/SVD.h
 Timestamp:
 Aug 31, 2006, 10:52:02 AM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/c++_tools/utility/SVD.h
r597 r616 1 // $Id$2 3 1 #ifndef _theplu_utility_svd_ 4 2 #define _theplu_utility_svd_ 5 3 6 #include <c++_tools/gslapi/matrix.h> 7 #include <c++_tools/gslapi/vector.h> 4 // $Id$ 5 6 #include <c++_tools/utility/matrix.h> 7 #include <c++_tools/utility/vector.h> 8 8 9 9 #include <gsl/gsl_linalg.h> … … 52 52 /// input matrix is copied for further use in the object. 53 53 /// 54 inline SVD(const gslapi::matrix& Ain)54 inline SVD(const utility::matrix& Ain) 55 55 : U_(Ain), V_(Ain.columns(),Ain.columns()), s_(Ain.columns()) {} 56 56 … … 73 73 /// is undefined. 74 74 /// 75 inline const gslapi::vector& s(void) const { return s_; }75 inline const utility::vector& s(void) const { return s_; } 76 76 77 77 /// … … 83 83 /// @return Whatever GSL returns. 84 84 /// 85 inline int solve( gslapi::vector b, gslapi::vector x)85 inline int solve(utility::vector b, utility::vector x) 86 86 { return gsl_linalg_SV_solve(U_.gsl_matrix_p(), V_.gsl_matrix_p(), 87 87 s_.gsl_vector_p(), b.gsl_vector_p(), … … 96 96 /// is undefined. 97 97 /// 98 inline const gslapi::matrix& U(void) const { return U_; }98 inline const utility::matrix& U(void) const { return U_; } 99 99 100 100 /// … … 106 106 /// is undefined. 107 107 /// 108 inline const gslapi::matrix& V(void) const { return V_; }108 inline const utility::matrix& V(void) const { return V_; } 109 109 110 110 private: … … 115 115 int modified_golub_reinsch(void); 116 116 117 gslapi::matrix U_, V_;118 gslapi::vector s_;117 utility::matrix U_, V_; 118 utility::vector s_; 119 119 }; 120 120
Note: See TracChangeset
for help on using the changeset viewer.