JACOUSTICS::TMatrixDS Struct Reference

Auxiliary data structure for compatibility of symmetric matrix. More...

#include <JKatoomba_t.hh>

Inheritance diagram for JACOUSTICS::TMatrixDS:

Public Member Functions
void	resize (const size_t size)
	Resize matrix. More...
void	reset ()
	Set matrix to the null matrix. More...
void	invert ()
	Invert matrix according SVD decomposition. More...
template<class JVectorND_t >
void	solve (JVectorND_t &u)
	Get solution of equation A x = b. More...

Static Public Attributes
static constexpr double	TOLERANCE = 1.0e-8
	Tolerance for SVD. More...
static std::mutex	mtx
	TDecompSVD. More...

Detailed Description

Auxiliary data structure for compatibility of symmetric matrix.

Inversion is based on SVD algorithm.

Definition at line 356 of file JKatoomba_t.hh.

Member Function Documentation

resize()

void JACOUSTICS::TMatrixDS::resize ( const size_t  size )

inline

Resize matrix.

Definition at line 364 of file JKatoomba_t.hh.

    {
      ResizeTo(size, size);
    }

reset()

void JACOUSTICS::TMatrixDS::reset ( )

inline

Set matrix to the null matrix.

Definition at line 372 of file JKatoomba_t.hh.

    {
      Zero();
    }

invert()

void JACOUSTICS::TMatrixDS::invert ( )

inline

Invert matrix according SVD decomposition.

Definition at line 380 of file JKatoomba_t.hh.

    {
#ifdef THREAD_SAFE
      using namespace std;
 
      unique_lock<mutex> lock(mtx);
#endif
      TDecompSVD svd(*this, TOLERANCE);
 
      Bool_t status;
 
      static_cast<TMatrixD&>(*this) = svd.Invert(status);
    }

solve()

template<class JVectorND_t >

void JACOUSTICS::TMatrixDS::solve ( JVectorND_t &  u )

inline

Get solution of equation A x = b.

Parameters

u	column vector; b on input, x on output(I/O)

Definition at line 400 of file JKatoomba_t.hh.

    {
#ifdef THREAD_SAFE
      using namespace std;
 
      unique_lock<mutex> lock(mtx);
#endif
      TDecompSVD svd(*this, TOLERANCE);
 
      TVectorD b(u.size());
 
      for (size_t i = 0; i != u.size(); ++i) {
        b[i] = u[i];
      }
 
      svd.Solve(b);
 
      for (size_t i = 0; i != u.size(); ++i) {
        u[i] = b[i];
      }
    }

Member Data Documentation

TOLERANCE

constexpr double JACOUSTICS::TMatrixDS::TOLERANCE = 1.0e-8

staticconstexpr

Tolerance for SVD.

Definition at line 359 of file JKatoomba_t.hh.

mtx

std::mutex JACOUSTICS::TMatrixDS::mtx

static

TDecompSVD.

Declaration of mutex.

Definition at line 422 of file JKatoomba_t.hh.

---

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

• JKatoomba_t.hh
• JKatoomba_t.cc