fn_solve.hpp
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030 template<typename eT, typename T1, typename T2>
00031 inline
00032 bool
00033 solve(Mat<eT>& X, const Base<eT,T1>& A_in, const Base<eT,T2>& B_in)
00034 {
00035 arma_extra_debug_sigprint();
00036
00037 const unwrap<T1> tmp1(A_in.get_ref());
00038 const unwrap<T2> tmp2(B_in.get_ref());
00039
00040 const Mat<eT>& A = tmp1.M;
00041 const Mat<eT>& B = tmp2.M;
00042
00043 arma_debug_check( (A.n_rows != B.n_rows), "solve(): number of rows in A and B must be the same" );
00044
00045 if(A.n_rows == A.n_cols)
00046 {
00047 return auxlib::solve(X, A, B);
00048 }
00049 else
00050 if(A.n_rows > A.n_cols)
00051 {
00052 arma_extra_debug_print("solve(): detected over-determined system");
00053 return auxlib::solve_od(X, A, B);
00054 }
00055 else
00056 {
00057 arma_extra_debug_print("solve(): detected under-determined system");
00058 return auxlib::solve_ud(X, A, B);
00059 }
00060 }
00061
00062
00063
00064 template<typename eT, typename T1, typename T2>
00065 inline
00066 Mat<eT>
00067 solve(const Base<eT,T1>& A_in, const Base<eT,T2>& B_in)
00068 {
00069 arma_extra_debug_sigprint();
00070
00071 Mat<eT> X;
00072 bool info = solve(X, A_in, B_in);
00073
00074 if(info == false)
00075 {
00076 arma_print("solve(): solution not found");
00077 }
00078
00079 return X;
00080 }
00081
00082
00083
00084