include/geometry/linalg.hpp
author Dmitriy Morozov <dmitriy@mrzv.org>
Sun, 10 Jun 2012 11:26:38 -0700
branchdev
changeset 269 86a02dac0db2
parent 19 efa14432761a
permissions -rw-r--r--
subcomplex in show_complex() + pan/zoom in ComplexViewer2D

template<class ValueType_>
typename LinearAlgebra<ValueType_>::ValueType 					
LinearAlgebra<ValueType_>::
determinant(const MatrixType& a)
{
	using namespace boost::numeric::ublas;
	MatrixType m = a;
	lu_factorize(m);
	return determinant(triangular_adaptor<MatrixType, upper>(m));
}

template<class ValueType_>
void
LinearAlgebra<ValueType_>::
solve(const MatrixType& a, const VectorType& b, VectorType& x)
{
	using namespace boost::numeric::ublas;
	MatrixType m = a;
	x = b;
	lu_factorize(m);
	//lu_substitute(m, x);
	inplace_solve(m,x,unit_lower_tag());
	inplace_solve(m,x,upper_tag());
}


/* Private */
template<class ValueType_>
template<class TriangularType_>
typename LinearAlgebra<ValueType_>::ValueType 					
LinearAlgebra<ValueType_>::
determinant(const boost::numeric::ublas::triangular_adaptor<MatrixType, TriangularType_>& t)
{
	ValueType res = ValueType(1);
	for (typename MatrixType::size_type i = 0; i < t.size1(); ++i)
		res *= t(i,i);
	return res;
}