JMATH::JMatrixND::JMatrixHelper Struct Reference

Auxiliary class for matrix decomposition. More...

#include <JMatrixND.hh>

Inheritance diagram for JMATH::JMatrixND::JMatrixHelper:

Public Types
typedef std::vector < std::vector< double > >	matrix_type
typedef matrix_type::const_iterator	const_row_type
typedef std::vector< double > ::const_iterator	const_col_type
typedef matrix_type::iterator	row_type
typedef std::vector< double > ::iterator	col_type
typedef T	data_type

Public Member Functions
const JMatrixHelper &	set (const JMatrixND &A)
	Set matrix helper. More...
void	decompose (const JMatrixND &A)
	LU decomposition. More...
int	getSign () const
	Get sign corresponding to permutation matrix. More...

Static Public Member Functions
static data_type &	getInstance ()
	Get unique instance of template class. More...

Private Attributes
int	sign

Detailed Description

Auxiliary class for matrix decomposition.

Definition at line 547 of file JMatrixND.hh.

Member Typedef Documentation

typedef std::vector< std::vector<double> > JMATH::JMatrixND_t::matrix_type

inherited

Definition at line 44 of file JMatrixND.hh.

typedef matrix_type::const_iterator JMATH::JMatrixND_t::const_row_type

inherited

Definition at line 46 of file JMatrixND.hh.

typedef std::vector<double>::const_iterator JMATH::JMatrixND_t::const_col_type

inherited

Definition at line 47 of file JMatrixND.hh.

typedef matrix_type::iterator JMATH::JMatrixND_t::row_type

inherited

Definition at line 49 of file JMatrixND.hh.

typedef std::vector<double>::iterator JMATH::JMatrixND_t::col_type

inherited

Definition at line 50 of file JMatrixND.hh.

template<class T>

typedef T JLANG::JSingleton< T >::data_type

inherited

Definition at line 20 of file JSingleton.hh.

Member Function Documentation

const JMatrixHelper& JMATH::JMatrixND::JMatrixHelper::set ( const JMatrixND &  A )

inline

Set matrix helper.

Parameters

A

matrix

Returns

this matrix helper

Definition at line 557 of file JMatrixND.hh.

      {
        static_cast<JMatrixND_t&>(*this) = A;

        return *this;
      }

void JMATH::JMatrixND::JMatrixHelper::decompose ( const JMatrixND &  A )

inline

LU decomposition.

Parameters

A

matrix

Definition at line 570 of file JMatrixND.hh.

      {
        set(A);

        sign = +1;
        
        size_t   k;
        double   val, del;
        row_type row, root;
        col_type col, cursor;
      
        // LU decomposition

        for (root = this->begin(), k = 0; root != this->end(); ++root, ++k) {

          row_type pivot;

          val   = (*root)[k];
          pivot = root;

          for (row = root; ++row != this->end(); ) {
            if (fabs((*row)[k]) > fabs(val)) {
              val   = (*row)[k];
              pivot = row;
            }
          }

          if (val == 0) {
            throw JDivisionByZero("LU decomposition zero pivot");
          }

          if (pivot != root) {

            iter_swap(root, pivot);

            sign = -sign;
          }


          // decompose in LU

          for (row = root; ++row != this->end(); ) {
          
            del = (*row)[k] /= val;
          
            for (col = row->begin() + k, cursor = root->begin() + k; ++col != row->end(); )
              (*col) -= del * (*(++cursor));
          }
        }
      }

int JMATH::JMatrixND::JMatrixHelper::getSign ( ) const

inline

Get sign corresponding to permutation matrix.

Definition at line 625 of file JMatrixND.hh.

      {
        return sign;
      }

template<class T>

static data_type& JLANG::JSingleton< T >::getInstance ( )

inlinestaticinherited

Get unique instance of template class.

Returns

object

Definition at line 27 of file JSingleton.hh.

    {
      static data_type value;

      return value;
    }

Member Data Documentation

int JMATH::JMatrixND::JMatrixHelper::sign

private

Definition at line 631 of file JMatrixND.hh.

---

The documentation for this struct was generated from the following file:

• JMatrixND.hh