author Dmitriy Morozov <>
Sat, 15 Sep 2007 14:31:36 -0400
changeset 34 19fb71f2405e
parent 19 efa14432761a
permissions -rw-r--r--
Counter omits color codes only if connected to a TTY

#ifndef __LINALG_H__
#define __LINALG_H__

#include <vector>

#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/lu.hpp>
#include <boost/numeric/ublas/io.hpp>

template<class ValueType_>
class LinearAlgebra
		typedef 					ValueType_												ValueType;
		typedef 					boost::numeric::ublas::matrix<ValueType>				MatrixType;
		typedef 					boost::numeric::ublas::vector<ValueType>				VectorType;

		/* Currently don't need any of this */
		static ValueType 			determinant(const MatrixType& a);
		static void					solve(const MatrixType& a, const VectorType& b, VectorType& x);

		template<class TriangularType_>
		static ValueType			determinant(const boost::numeric::ublas::triangular_adaptor<MatrixType, TriangularType_>& t);

#include "linalg.hpp"