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 415 of file JKatoomba_t.hh.

Member Function Documentation

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

inline

Resize matrix.

Definition at line 423 of file JKatoomba_t.hh.

    {
      ResizeTo(size, size);
    }

void JACOUSTICS::TMatrixDS::reset ( )

inline

Set matrix to the null matrix.

Definition at line 431 of file JKatoomba_t.hh.

    {
      Zero();
    }

void JACOUSTICS::TMatrixDS::invert ( )

inline

Invert matrix according SVD decomposition.

Definition at line 439 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);
    }

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 459 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

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

static

Tolerance for SVD.

Definition at line 418 of file JKatoomba_t.hh.

std::mutex JACOUSTICS::TMatrixDS::mtx

static

TDecompSVD.

Declaration of mutex.

Definition at line 481 of file JKatoomba_t.hh.

---

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

• JKatoomba_t.hh
• JKatoomba_t.cc