45 template <
class MaTriX,
class VecToR>
52 BIASASSERT(A.num_cols() == N);
53 BIASASSERT(b.dim() == N);
60 typename MaTriX::element_type tmp=0;
63 tmp = tmp + A(i,j)*x(j);
65 x(i) = (b(i) - tmp)/ A(i,i);
72 template <
class MaTriX,
class VecToR>
79 BIASASSERT(A.num_cols() == N);
80 BIASASSERT(b.dim() == N);
88 typename MaTriX::element_type tmp=0;
91 tmp = tmp + A(i,j)*x(j);
100 template <
class MaTriX,
class VecToR>
107 template <
class MaTriX,
class VecToR>
118 template <
class MaTriX,
class VecToR>
125 BIASASSERT(A.num_cols() == N);
126 BIASASSERT(b.dim() == N);
134 typename MaTriX::element_type tmp=0;
137 tmp = tmp + A(i,j)*x(j);
139 x(i) = (b(i) - tmp)/ A(i,i);
146 template <
class MaTriX,
class VecToR>
153 BIASASSERT(A.num_cols() == N);
154 BIASASSERT(b.dim() == N);
162 typename MaTriX::element_type tmp=0;
165 tmp = tmp + A(i,j)*x(j);
174 template <
class MaTriX,
class VecToR>
181 template <
class MaTriX,
class VecToR>
VecToR Lower_triangular_solve(MaTriX &A, VecToR &b)
VecToR Upper_triangular_solve(MaTriX &A, VecToR &b)
TNT_SUBSCRIPT_TYPE Subscript
VecToR Unit_lower_triangular_solve(MaTriX &A, VecToR &b)
VecToR Unit_upper_triangular_solve(MaTriX &A, VecToR &b)
VecToR linear_solve(LowerTriangularView< MaTriX > &A, VecToR &b)