NonlinearSystem.h File Reference

#include <iostream>
#include <sstream>
#include <functional>
#include <cerrno>
#include <iomanip>
#include <limits>
#include <algorithm>
#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/lu.hpp>
#include <boost/numeric/ublas/vector.hpp>
#include <boost/numeric/ublas/io.hpp>
#include "VoxelPools.h"

Include dependency graph for NonlinearSystem.h:

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes
class	NonlinearSystem
class	ReacInfo

Typedefs
typedef ublas::matrix< value_type >	matrix_type
typedef double	value_type
typedef ublas::vector< value_type >	vector_type

Functions
template<class T >
bool	inverse (const ublas::matrix< T > &input, ublas::matrix< T > &inverse)

Typedef Documentation

typedef ublas::matrix<value_type> matrix_type

Definition at line 43 of file NonlinearSystem.h.

typedef double value_type

Definition at line 41 of file NonlinearSystem.h.

typedef ublas::vector<value_type> vector_type

Definition at line 42 of file NonlinearSystem.h.

Function Documentation

template<class T >

bool inverse	(	const ublas::matrix< T > &	input,
		ublas::matrix< T > &	inverse
	)

Definition at line 65 of file NonlinearSystem.h.

Referenced by NonlinearSystem::compute_jacobians().

{
    using namespace boost::numeric::ublas;
    typedef permutation_matrix<std::size_t> pmatrix;
    // create a working copy of the input
    matrix<T> A(input);
    // create a permutation matrix for the LU-factorization
    pmatrix pm(A.size1());

    // perform LU-factorization
    int res = lu_factorize(A,pm);
    if( res != 0 ) return false;

    // create identity matrix of "inverse"
    inverse.assign(ublas::identity_matrix<T>(A.size1()));

    // backsubstitute to get the inverse
    lu_substitute(A, pm, inverse);

    return true;
}

Here is the caller graph for this function: