MATRIX::VtSqMatrix Class Reference

#include <VtSqMatrix.hh>

Inheritance diagram for MATRIX::VtSqMatrix:

Collaboration diagram for MATRIX::VtSqMatrix:

Public Member Functions
— Constructors —
	VtSqMatrix (const int row)
	VtSqMatrix (const VtSqMatrix &rhs)
	VtSqMatrix (const VtMatrix &rhs)
virtual	~VtSqMatrix ()
— Matrix operations —
bool	VtDinv (int dim=0)
	transform to inverse More...
const VtSqMatrix	dinv (int dim=0) const
	return inverse More...
virtual double	det () const
	return determinant More...
virtual const VtMatrix	operator* (const VtMatrix &rhs)
const VtVector	operator* (const VtVector &rhs) const
Public Member Functions inherited from MATRIX::VtMatrix
	VtMatrix (const unsigned int row, const unsigned int col)
	VtMatrix (const VtMatrix &rhs)
virtual	~VtMatrix ()
unsigned int	nrow () const
	no of rows $n$ More...
unsigned int	ncol () const
	no of columns $m$ More...
int	size () const
	$m\times n$ More...
VtMatrix_row	operator[] (int row)
VtMatrix_row_const	operator[] (int row) const
virtual double	operator() (unsigned int row, unsigned int col) const
virtual double &	operator() (const unsigned int row, const unsigned int col)
double	get (unsigned int row, unsigned int col) const
double &	get (unsigned int row, unsigned int col)
void	VtT (void)
	transform into transpose matrix More...
const VtMatrix	T (void) const
	return transpose More...
virtual void	place_at (const VtMatrix &rhs, const unsigned int row, const unsigned int col)
	copy a smaller matrix at a certain place More...
virtual void	place_at (const VtVector &rhs, const unsigned int row, const unsigned int col)
	copy a vector at a certain place More...
void	copy (const VtMatrix &rhs)
	to be used if matrix dimensions are not equal More...
void	clear (void)
	set matrix elements to 0 More...
const VtMatrix &	operator= (const VtMatrix &rhs)
	$\textbf{A} = \textbf{B}$ More...
const VtMatrix &	operator= (const VtNegMatrix &rhs)
virtual const VtMatrix &	operator+= (const double rhs)
	$\textbf{A} = (a_{\mu\nu} + \alpha)$ More...
virtual const VtMatrix &	operator-= (const double rhs)
	$\textbf{A} = (a_{\mu\nu} - \alpha)$ More...
virtual const VtMatrix &	operator*= (const double rhs)
	$\textbf{A} = (a_{\mu\nu} \cdot\alpha)$ More...
virtual const VtMatrix &	operator/= (const double rhs)
	$\textbf{A} = (a_{\mu\nu} / \alpha)$ More...
const VtMatrix &	operator+= (const VtMatrix &rhs)
	$\textbf{A} = (a_{\mu\nu} + b_{\mu\nu})$ More...
const VtMatrix &	operator-= (const VtMatrix &rhs)
	$\textbf{A} = (a_{\mu\nu} - b_{\mu\nu})$ More...
const VtMatrix	operator+ (const VtMatrix &rhs) const
	$\textbf{A} + \textbf{B}$ More...
const VtMatrix	operator+ (const VtNegMatrix &rhs) const
const VtMatrix	operator- (const VtMatrix &rhs) const
	$\textbf{A} - \textbf{B}$ More...
const VtMatrix	operator- (const VtNegMatrix &rhs) const
const VtNegMatrix	operator- (void) const
	$-\textbf{A}$ More...
const VtMatrix	operator* (const VtMatrix &rhs) const
	$\textbf{A}\cdot\textbf{B} = \sum_{\nu=1}^n a_{\mu\nu}b_{\nu\lambda}$ More...
const VtVector	operator* (const VtVector &rhs) const
	$\textbf{A}\cdot\vec{v} = (\sum_{\nu=1}^n a_{\mu\nu}v_{\nu})$ More...
double *	array () const
	return pointer to internal array More...
virtual void	print (std::ostream &os) const

Additional Inherited Members
Protected Attributes inherited from MATRIX::VtMatrix
double *	m
double *	work
unsigned int	m_nrow
unsigned int	m_ncol

Detailed Description

Class for quadratic matrices

Constructor & Destructor Documentation

VtSqMatrix() [1/3]

MATRIX::VtSqMatrix::VtSqMatrix ( const int  row )

inline

: VtMatrix(row,row) {}

VtSqMatrix() [2/3]

MATRIX::VtSqMatrix::VtSqMatrix ( const VtSqMatrix &  rhs )

inline

: VtMatrix(rhs) {}

VtSqMatrix() [3/3]

MATRIX::VtSqMatrix::VtSqMatrix ( const VtMatrix &  rhs )

inline

                                    : VtMatrix(rhs) { 
#ifndef VtFAST
      assert(rhs.ncol()==rhs.nrow()); 
#endif
    }

~VtSqMatrix()

MATRIX::VtSqMatrix::~VtSqMatrix ( )

virtual

{}

Member Function Documentation

det()

double MATRIX::VtSqMatrix::det ( void  ) const

virtual

return determinant

Reimplemented in MATRIX::VtSymMatrix.

                                   {
    
    int ifail, jfail;
    int n = nrow();
    double det;
 
    memcpy( work, m, size() * sizeof(double) );
    double *work2 = new double[n];
 
    Dfact1(&n, work, &n, work2, &ifail, &det, &jfail);
    if(jfail != 0) {
      det = 0;
#ifdef VtDEBUG
      cout << "VtMatrix::det() error: Determinant computation failed!" << endl;
#endif
    }
 
    delete[] work2;
   
    return det;
  }

dinv()

const VtSqMatrix MATRIX::VtSqMatrix::dinv

(

int

dim = 0

)

const

return inverse

                                                 {
    VtSqMatrix tmp(*this);
    tmp.VtDinv(dim);
    return tmp;
  }

operator*() [1/2]

const VtMatrix MATRIX::VtSqMatrix::operator* ( const VtMatrix &  rhs )

virtual

                                                          {
    return VtMatrix::operator*(rhs);
  }

operator*() [2/2]

const VtVector MATRIX::VtSqMatrix::operator*

(

const VtVector &

rhs

)

const

                                                                {
    return VtMatrix::operator*(rhs);
  }

VtDinv()

bool MATRIX::VtSqMatrix::VtDinv

(

int

dim = 0

)

transform to inverse

                                 {
    bool valid = false;
 
    if(dim == 0) dim = nrow();
#ifndef VtFAST
    else {
      assert(dim <= static_cast<int>(nrow()));
    }
#endif
 
    const size_t msize = size() * sizeof(double);
 
    memcpy( work, m, msize );
 
    int    ifail = 0;
    double *work2 = new double[dim];
    Dinv1(&dim, work, &dim, work2, &ifail);
    delete[] work2;
 
    if(ifail != -1) {
      memcpy( m, work, msize );
      valid = true;
    }
#ifdef VtDEBUG
    else {
      cout << "VtMatrix::dinv() error: Matrix inversion failed!!!" << endl;
    }
#endif
 
    return valid;
  } 

---

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

• src/libVt++/vt++/include/VtSqMatrix.hh
• src/libVt++/vt++/VtSqMatrix.C