Logo ROOT   6.13/01
Reference Guide
List of all members | Public Types | Public Member Functions | Protected Member Functions | Protected Attributes | List of all members
TMatrixT< Element > Class Template Reference

template<class Element>
class TMatrixT< Element >

TMatrixT.

Template class of a general matrix in the linear algebra package

Definition at line 22 of file TMatrixDfwd.h.

Public Types

enum  { kWorkMax = 100 }
 
enum  EMatrixCreatorsOp1 {
  kZero, kUnit, kTransposed, kInverted,
  kAtA
}
 
enum  EMatrixCreatorsOp2 {
  kMult, kTransposeMult, kInvMult, kMultTranspose,
  kPlus, kMinus
}
 

Public Member Functions

 TMatrixT ()
 
 TMatrixT (Int_t nrows, Int_t ncols)
 Constructor for (nrows x ncols) matrix. More...
 
 TMatrixT (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb)
 Constructor for ([row_lwb..row_upb] x [col_lwb..col_upb]) matrix. More...
 
 TMatrixT (Int_t nrows, Int_t ncols, const Element *data, Option_t *option="")
 option="F": array elements contains the matrix stored column-wise like in Fortran, so a[i,j] = elements[i+no_rows*j], else it is supposed that array elements are stored row-wise a[i,j] = elements[i*no_cols+j] More...
 
 TMatrixT (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, const Element *data, Option_t *option="")
 array elements are copied More...
 
 TMatrixT (const TMatrixT< Element > &another)
 Copy constructor. More...
 
 TMatrixT (const TMatrixTSym< Element > &another)
 Copy constructor of a symmetric matrix. More...
 
 TMatrixT (const TMatrixTSparse< Element > &another)
 Copy constructor of a sparse matrix. More...
 
template<class Element2 >
 TMatrixT (const TMatrixT< Element2 > &another)
 
 TMatrixT (EMatrixCreatorsOp1 op, const TMatrixT< Element > &prototype)
 Constructor of matrix applying a specific operation to the prototype. More...
 
 TMatrixT (const TMatrixT< Element > &a, EMatrixCreatorsOp2 op, const TMatrixT< Element > &b)
 Constructor of matrix applying a specific operation to two prototypes. More...
 
 TMatrixT (const TMatrixT< Element > &a, EMatrixCreatorsOp2 op, const TMatrixTSym< Element > &b)
 Constructor of matrix applying a specific operation to two prototypes. More...
 
 TMatrixT (const TMatrixTSym< Element > &a, EMatrixCreatorsOp2 op, const TMatrixT< Element > &b)
 Constructor of matrix applying a specific operation to two prototypes. More...
 
 TMatrixT (const TMatrixTSym< Element > &a, EMatrixCreatorsOp2 op, const TMatrixTSym< Element > &b)
 Constructor of matrix applying a specific operation to two prototypes. More...
 
 TMatrixT (const TMatrixTLazy< Element > &lazy_constructor)
 Constructor using the TMatrixTLazy class. More...
 
virtual ~TMatrixT ()
 
virtual void Clear (Option_t *="")
 
virtual Double_t Determinant () const
 Return the matrix determinant. More...
 
virtual void Determinant (Double_t &d1, Double_t &d2) const
 Return the matrix determinant as d1,d2 where det = d1*TMath::Power(2.0,d2) More...
 
const TMatrixT< Element > EigenVectors (TVectorT< Element > &eigenValues) const
 Return a matrix containing the eigen-vectors ordered by descending values of Re^2+Im^2 of the complex eigen-values . More...
 
virtual const Int_t * GetColIndexArray () const
 
virtual Int_t * GetColIndexArray ()
 
virtual const Element * GetMatrixArray () const
 
virtual Element * GetMatrixArray ()
 
virtual const Int_t * GetRowIndexArray () const
 
virtual Int_t * GetRowIndexArray ()
 
virtual TMatrixTBase< Element > & GetSub (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, TMatrixTBase< Element > &target, Option_t *option="S") const
 Get submatrix [row_lwb..row_upb] x [col_lwb..col_upb]; The indexing range of the returned matrix depends on the argument option: More...
 
TMatrixT< Element > GetSub (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, Option_t *option="S") const
 
TMatrixT< Element > & Invert (Double_t *det=0)
 Invert the matrix and calculate its determinant. More...
 
TMatrixT< Element > & InvertFast (Double_t *det=0)
 Invert the matrix and calculate its determinant, however upto (6x6) a fast Cramer inversion is used . More...
 
void Minus (const TMatrixT< Element > &a, const TMatrixT< Element > &b)
 General matrix summation. Create a matrix C such that C = A - B. More...
 
void Minus (const TMatrixT< Element > &a, const TMatrixTSym< Element > &b)
 General matrix summation. Create a matrix C such that C = A - B. More...
 
void Minus (const TMatrixTSym< Element > &a, const TMatrixT< Element > &b)
 
void Mult (const TMatrixT< Element > &a, const TMatrixT< Element > &b)
 General matrix multiplication. Create a matrix C such that C = A * B. More...
 
void Mult (const TMatrixT< Element > &a, const TMatrixTSym< Element > &b)
 Matrix multiplication, with A general and B symmetric. More...
 
void Mult (const TMatrixTSym< Element > &a, const TMatrixT< Element > &b)
 Matrix multiplication, with A symmetric and B general. More...
 
void Mult (const TMatrixTSym< Element > &a, const TMatrixTSym< Element > &b)
 Matrix multiplication, with A symmetric and B symmetric. More...
 
void MultT (const TMatrixT< Element > &a, const TMatrixT< Element > &b)
 General matrix multiplication. Create a matrix C such that C = A * B^T. More...
 
void MultT (const TMatrixT< Element > &a, const TMatrixTSym< Element > &b)
 
void MultT (const TMatrixTSym< Element > &a, const TMatrixT< Element > &b)
 Matrix multiplication, with A symmetric and B general. More...
 
void MultT (const TMatrixTSym< Element > &a, const TMatrixTSym< Element > &b)
 
TMatrixT< Element > & NormByColumn (const TVectorT< Element > &v, Option_t *option="D")
 Multiply/divide matrix columns by a vector: option: "D" : b(i,j) = a(i,j)/v(i) i = 0,fNrows-1 (default) else : b(i,j) = a(i,j)*v(i) More...
 
TMatrixT< Element > & NormByRow (const TVectorT< Element > &v, Option_t *option="D")
 Multiply/divide matrix rows with a vector: option: "D" : b(i,j) = a(i,j)/v(j) i = 0,fNcols-1 (default) else : b(i,j) = a(i,j)*v(j) More...
 
Element operator() (Int_t rown, Int_t coln) const
 
Element & operator() (Int_t rown, Int_t coln)
 
TMatrixT< Element > & operator*= (Element val)
 Multiply every element of the matrix with val. More...
 
TMatrixT< Element > & operator*= (const TMatrixT< Element > &source)
 Compute target = target * source inplace. More...
 
TMatrixT< Element > & operator*= (const TMatrixTSym< Element > &source)
 Compute target = target * source inplace. More...
 
TMatrixT< Element > & operator*= (const TMatrixTDiag_const< Element > &diag)
 Multiply a matrix row by the diagonal of another matrix matrix(i,j) *= diag(j), j=0,fNcols-1. More...
 
TMatrixT< Element > & operator*= (const TMatrixTRow_const< Element > &row)
 Multiply a matrix by the row of another matrix matrix(i,j) *= another(k,j) for fixed k. More...
 
TMatrixT< Element > & operator*= (const TMatrixTColumn_const< Element > &col)
 Multiply a matrix by the column of another matrix matrix(i,j) *= another(i,k) for fixed k. More...
 
TMatrixT< Element > & operator+= (Element val)
 Add val to every element of the matrix. More...
 
TMatrixT< Element > & operator+= (const TMatrixT< Element > &source)
 Add the source matrix. More...
 
TMatrixT< Element > & operator+= (const TMatrixTSym< Element > &source)
 Add the source matrix. More...
 
TMatrixT< Element > & operator-= (Element val)
 Subtract val from every element of the matrix. More...
 
TMatrixT< Element > & operator-= (const TMatrixT< Element > &source)
 Subtract the source matrix. More...
 
TMatrixT< Element > & operator-= (const TMatrixTSym< Element > &source)
 Subtract the source matrix. More...
 
TMatrixT< Element > & operator/= (const TMatrixTDiag_const< Element > &diag)
 Divide a matrix row by the diagonal of another matrix matrix(i,j) /= diag(j) More...
 
TMatrixT< Element > & operator/= (const TMatrixTRow_const< Element > &row)
 Divide a matrix by the row of another matrix matrix(i,j) /= another(k,j) for fixed k. More...
 
TMatrixT< Element > & operator/= (const TMatrixTColumn_const< Element > &col)
 Divide a matrix by the column of another matrix matrix(i,j) /= another(i,k) for fixed k. More...
 
TMatrixT< Element > & operator= (const TMatrixT< Element > &source)
 Assignment operator. More...
 
TMatrixT< Element > & operator= (const TMatrixTSym< Element > &source)
 Assignment operator. More...
 
TMatrixT< Element > & operator= (const TMatrixTSparse< Element > &source)
 Assignment operator. More...
 
TMatrixT< Element > & operator= (const TMatrixTLazy< Element > &source)
 Assignment operator. More...
 
template<class Element2 >
TMatrixT< Element > & operator= (const TMatrixT< Element2 > &source)
 
TMatrixT< Element > & operator= (Element val)
 Assign val to every element of the matrix. More...
 
const TMatrixTRow_const< Element > operator[] (Int_t rown) const
 
TMatrixTRow< Element > operator[] (Int_t rown)
 
void Plus (const TMatrixT< Element > &a, const TMatrixT< Element > &b)
 General matrix summation. Create a matrix C such that C = A + B. More...
 
void Plus (const TMatrixT< Element > &a, const TMatrixTSym< Element > &b)
 General matrix summation. Create a matrix C such that C = A + B. More...
 
void Plus (const TMatrixTSym< Element > &a, const TMatrixT< Element > &b)
 
TMatrixT< Element > & Rank1Update (const TVectorT< Element > &v, Element alpha=1.0)
 Perform a rank 1 operation on matrix A: A += alpha * v * v^T. More...
 
TMatrixT< Element > & Rank1Update (const TVectorT< Element > &v1, const TVectorT< Element > &v2, Element alpha=1.0)
 Perform a rank 1 operation on matrix A: A += alpha * v1 * v2^T. More...
 
virtual TMatrixTBase< Element > & ResizeTo (Int_t nrows, Int_t ncols, Int_t=-1)
 Set size of the matrix to nrows x ncols New dynamic elements are created, the overlapping part of the old ones are copied to the new structures, then the old elements are deleted. More...
 
virtual TMatrixTBase< Element > & ResizeTo (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, Int_t=-1)
 Set size of the matrix to [row_lwb:row_upb] x [col_lwb:col_upb] New dynamic elemenst are created, the overlapping part of the old ones are copied to the new structures, then the old elements are deleted. More...
 
TMatrixTBase< Element > & ResizeTo (const TMatrixT< Element > &m)
 
virtual TMatrixTBase< Element > & SetColIndexArray (Int_t *)
 
virtual TMatrixTBase< Element > & SetRowIndexArray (Int_t *)
 
virtual TMatrixTBase< Element > & SetSub (Int_t row_lwb, Int_t col_lwb, const TMatrixTBase< Element > &source)
 Insert matrix source starting at [row_lwb][col_lwb], thereby overwriting the part [row_lwb..row_lwb+nrows_source][col_lwb..col_lwb+ncols_source];. More...
 
Element Similarity (const TVectorT< Element > &v) const
 Calculate scalar v * (*this) * v^T. More...
 
TMatrixT< Element > & T ()
 
void TMult (const TMatrixT< Element > &a, const TMatrixT< Element > &b)
 Create a matrix C such that C = A' * B. More...
 
void TMult (const TMatrixT< Element > &a, const TMatrixTSym< Element > &b)
 Create a matrix C such that C = A' * B. More...
 
void TMult (const TMatrixTSym< Element > &a, const TMatrixT< Element > &b)
 
void TMult (const TMatrixTSym< Element > &a, const TMatrixTSym< Element > &b)
 
TMatrixT< Element > & Transpose (const TMatrixT< Element > &source)
 Transpose matrix source. More...
 
TMatrixT< Element > & Use (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, Element *data)
 Use the array data to fill the matrix ([row_lwb..row_upb] x [col_lwb..col_upb]) More...
 
const TMatrixT< Element > & Use (Int_t row_lwb, Int_t row_upb, Int_t col_lwb, Int_t col_upb, const Element *data) const
 
TMatrixT< Element > & Use (Int_t nrows, Int_t ncols, Element *data)
 
const TMatrixT< Element > & Use (Int_t nrows, Int_t ncols, const Element *data) const
 
TMatrixT< Element > & Use (TMatrixT< Element > &a)
 
const TMatrixT< Element > & Use (const TMatrixT< Element > &a) const
 
- Public Member Functions inherited from TMatrixTBase< Element >
 TMatrixTBase ()
 
virtual ~TMatrixTBase ()
 
virtual TMatrixTBase< Element > & Abs ()
 Take an absolute value of a matrix, i.e. apply Abs() to each element. More...
 
virtual TMatrixTBase< Element > & Apply (const TElementActionT< Element > &action)
 Apply action to each matrix element. More...
 
virtual TMatrixTBase< Element > & Apply (const TElementPosActionT< Element > &action)
 Apply action to each element of the matrix. More...
 
virtual Element ColNorm () const
 Column matrix norm, MAX{ SUM{ |M(i,j)|, over i}, over j}. More...
 
void Draw (Option_t *option="")
 Draw this matrix The histogram is named "TMatrixT" by default and no title. More...
 
virtual Element E2Norm () const
 Square of the Euclidian norm, SUM{ m(i,j)^2 }. More...
 
virtual void ExtractRow (Int_t row, Int_t col, Element *v, Int_t n=-1) const
 Store in array v, n matrix elements of row rown starting at column coln. More...
 
Int_t GetColLwb () const
 
Int_t GetColUpb () const
 
virtual void GetMatrix2Array (Element *data, Option_t *option="") const
 Copy matrix data to array . More...
 
Int_t GetNcols () const
 
Int_t GetNoElements () const
 
Int_t GetNrows () const
 
Int_t GetRowLwb () const
 
Int_t GetRowUpb () const
 
Element GetTol () const
 
virtual TMatrixTBase< Element > & InsertRow (Int_t row, Int_t col, const Element *v, Int_t n=-1)
 Copy n elements from array v to row rown starting at column coln. More...
 
void Invalidate ()
 
Bool_t IsOwner () const
 
virtual Bool_t IsSymmetric () const
 Check whether matrix is symmetric. More...
 
Bool_t IsValid () const
 
void MakeValid ()
 
virtual Element Max () const
 return maximum vector element value More...
 
virtual Element Min () const
 return minimum matrix element value More...
 
virtual Int_t NonZeros () const
 Compute the number of elements != 0.0. More...
 
Element Norm1 () const
 
virtual TMatrixTBase< Element > & NormByDiag (const TVectorT< Element > &v, Option_t *option="D")
 option: More...
 
Element NormInf () const
 
Bool_t operator!= (Element val) const
 Are all matrix elements not equal to val? More...
 
Bool_t operator< (Element val) const
 Are all matrix elements < val? More...
 
Bool_t operator<= (Element val) const
 Are all matrix elements <= val? More...
 
Bool_t operator== (Element val) const
 Are all matrix elements equal to val? More...
 
Bool_t operator> (Element val) const
 Are all matrix elements > val? More...
 
Bool_t operator>= (Element val) const
 Are all matrix elements >= val? More...
 
void Print (Option_t *name="") const
 Print the matrix as a table of elements. More...
 
virtual TMatrixTBase< Element > & Randomize (Element alpha, Element beta, Double_t &seed)
 Randomize matrix element values. More...
 
virtual Element RowNorm () const
 Row matrix norm, MAX{ SUM{ |M(i,j)|, over j}, over i}. More...
 
virtual TMatrixTBase< Element > & SetMatrixArray (const Element *data, Option_t *option="")
 Copy array data to matrix . More...
 
Element SetTol (Element tol)
 
virtual TMatrixTBase< Element > & Shift (Int_t row_shift, Int_t col_shift)
 Shift the row index by adding row_shift and the column index by adding col_shift, respectively. More...
 
virtual TMatrixTBase< Element > & Sqr ()
 Square each element of the matrix. More...
 
virtual TMatrixTBase< Element > & Sqrt ()
 Take square root of all elements. More...
 
virtual Element Sum () const
 Compute sum of elements. More...
 
virtual TMatrixTBase< Element > & UnitMatrix ()
 Make a unit matrix (matrix need not be a square one). More...
 
virtual TMatrixTBase< Element > & Zero ()
 Set matrix elements to zero. More...
 

Protected Member Functions

void Allocate (Int_t nrows, Int_t ncols, Int_t row_lwb=0, Int_t col_lwb=0, Int_t init=0, Int_t=-1)
 Allocate new matrix. More...
 
void Delete_m (Int_t size, Element *&)
 Delete data pointer m, if it was assigned on the heap. More...
 
Int_t Memcpy_m (Element *newp, const Element *oldp, Int_t copySize, Int_t newSize, Int_t oldSize)
 Copy copySize doubles from *oldp to *newp . More...
 
Element * New_m (Int_t size)
 Return data pointer . More...
 

Protected Attributes

Element fDataStack [TMatrixTBase< Element >::kSizeMax]
 
Element * fElements
 data container More...
 
- Protected Attributes inherited from TMatrixTBase< Element >
Int_t fColLwb
 
Bool_t fIsOwner
 
Int_t fNcols
 
Int_t fNelems
 
Int_t fNrowIndex
 
Int_t fNrows
 
Int_t fRowLwb
 
Element fTol
 

Additional Inherited Members

- Static Public Member Functions inherited from TMatrixTBase< Element >
static Element & NaNValue ()
 
- Protected Types inherited from TMatrixTBase< Element >
enum  { kSizeMax = 25 }
 
enum  { kWorkMax = 100 }
 
enum  EMatrixStatusBits { kStatus = BIT(14) }
 
- Static Protected Member Functions inherited from TMatrixTBase< Element >
static void DoubleLexSort (Int_t n, Int_t *first, Int_t *second, Element *data)
 default kTRUE, when Use array kFALSE More...
 
static void IndexedLexSort (Int_t n, Int_t *first, Int_t swapFirst, Int_t *second, Int_t swapSecond, Int_t *index)
 Lexical sort on array data using indices first and second. More...
 

#include <TMatrixDfwd.h>

Inheritance diagram for TMatrixT< Element >:
[legend]

Member Enumeration Documentation

◆ anonymous enum

template<class Element>
anonymous enum
Enumerator
kWorkMax 

Definition at line 57 of file TMatrixT.h.

◆ EMatrixCreatorsOp1

template<class Element>
enum TMatrixT::EMatrixCreatorsOp1
Enumerator
kZero 
kUnit 
kTransposed 
kInverted 
kAtA 

Definition at line 58 of file TMatrixT.h.

◆ EMatrixCreatorsOp2

template<class Element>
enum TMatrixT::EMatrixCreatorsOp2
Enumerator
kMult 
kTransposeMult 
kInvMult 
kMultTranspose 
kPlus 
kMinus 

Definition at line 59 of file TMatrixT.h.

Constructor & Destructor Documentation

◆ TMatrixT() [1/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( )
inline

Definition at line 61 of file TMatrixT.h.

◆ TMatrixT() [2/15]

template<class Element >
TMatrixT< Element >::TMatrixT ( Int_t  nrows,
Int_t  ncols 
)

Constructor for (nrows x ncols) matrix.

Definition at line 39 of file TMatrixT.cxx.

◆ TMatrixT() [3/15]

template<class Element >
TMatrixT< Element >::TMatrixT ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb 
)

Constructor for ([row_lwb..row_upb] x [col_lwb..col_upb]) matrix.

Definition at line 48 of file TMatrixT.cxx.

◆ TMatrixT() [4/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( Int_t  no_rows,
Int_t  no_cols,
const Element *  elements,
Option_t *  option = "" 
)

option="F": array elements contains the matrix stored column-wise like in Fortran, so a[i,j] = elements[i+no_rows*j], else it is supposed that array elements are stored row-wise a[i,j] = elements[i*no_cols+j]

array elements are copied

Definition at line 62 of file TMatrixT.cxx.

◆ TMatrixT() [5/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
const Element *  data,
Option_t *  option = "" 
)

array elements are copied

Definition at line 72 of file TMatrixT.cxx.

◆ TMatrixT() [6/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixT< Element > &  another)

Copy constructor.

Definition at line 83 of file TMatrixT.cxx.

◆ TMatrixT() [7/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixTSym< Element > &  another)

Copy constructor of a symmetric matrix.

Definition at line 94 of file TMatrixT.cxx.

◆ TMatrixT() [8/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixTSparse< Element > &  another)

Copy constructor of a sparse matrix.

Definition at line 105 of file TMatrixT.cxx.

◆ TMatrixT() [9/15]

template<class Element>
template<class Element2 >
TMatrixT< Element >::TMatrixT ( const TMatrixT< Element2 > &  another)
inline

Definition at line 69 of file TMatrixT.h.

◆ TMatrixT() [10/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( EMatrixCreatorsOp1  op,
const TMatrixT< Element > &  prototype 
)

Constructor of matrix applying a specific operation to the prototype.

Example: TMatrixT<Element> a(10,12); ...; TMatrixT<Element> b(TMatrixT::kTransposed, a); Supported operations are: kZero, kUnit, kTransposed, kInverted and kAtA.

Definition at line 118 of file TMatrixT.cxx.

◆ TMatrixT() [11/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixT< Element > &  a,
EMatrixCreatorsOp2  op,
const TMatrixT< Element > &  b 
)

Constructor of matrix applying a specific operation to two prototypes.

Example: TMatrixT<Element> a(10,12), b(12,5); ...; TMatrixT<Element> c(a, TMatrixT::kMult, b); Supported operations are: kMult (a*b), kTransposeMult (a'*b), kInvMult (a^(-1)*b)

Definition at line 169 of file TMatrixT.cxx.

◆ TMatrixT() [12/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixT< Element > &  a,
EMatrixCreatorsOp2  op,
const TMatrixTSym< Element > &  b 
)

Constructor of matrix applying a specific operation to two prototypes.

Example: TMatrixT<Element> a(10,12), b(12,5); ...; TMatrixT<Element> c(a, TMatrixT::kMult, b); Supported operations are: kMult (a*b), kTransposeMult (a'*b), kInvMult (a^(-1)*b)

Definition at line 226 of file TMatrixT.cxx.

◆ TMatrixT() [13/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixTSym< Element > &  a,
EMatrixCreatorsOp2  op,
const TMatrixT< Element > &  b 
)

Constructor of matrix applying a specific operation to two prototypes.

Example: TMatrixT<Element> a(10,12), b(12,5); ...; TMatrixT<Element> c(a, TMatrixT::kMult, b); Supported operations are: kMult (a*b), kTransposeMult (a'*b), kInvMult (a^(-1)*b)

Definition at line 283 of file TMatrixT.cxx.

◆ TMatrixT() [14/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixTSym< Element > &  a,
EMatrixCreatorsOp2  op,
const TMatrixTSym< Element > &  b 
)

Constructor of matrix applying a specific operation to two prototypes.

Example: TMatrixT<Element> a(10,12), b(12,5); ...; TMatrixT<Element> c(a, TMatrixT::kMult, b); Supported operations are: kMult (a*b), kTransposeMult (a'*b), kInvMult (a^(-1)*b)

Definition at line 340 of file TMatrixT.cxx.

◆ TMatrixT() [15/15]

template<class Element>
TMatrixT< Element >::TMatrixT ( const TMatrixTLazy< Element > &  lazy_constructor)

Constructor using the TMatrixTLazy class.

Definition at line 395 of file TMatrixT.cxx.

◆ ~TMatrixT()

template<class Element>
virtual TMatrixT< Element >::~TMatrixT ( )
inlinevirtual

Definition at line 83 of file TMatrixT.h.

Member Function Documentation

◆ Allocate()

template<class Element >
void TMatrixT< Element >::Allocate ( Int_t  no_rows,
Int_t  no_cols,
Int_t  row_lwb = 0,
Int_t  col_lwb = 0,
Int_t  init = 0,
Int_t  = -1 
)
protected

Allocate new matrix.

Arguments are number of rows, columns, row lowerbound (0 default) and column lowerbound (0 default).

Definition at line 466 of file TMatrixT.cxx.

◆ Clear()

template<class Element>
virtual void TMatrixT< Element >::Clear ( Option_t *  = "")
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 120 of file TMatrixT.h.

◆ Delete_m()

template<class Element>
void TMatrixT< Element >::Delete_m ( Int_t  size,
Element *&  m 
)
protected

Delete data pointer m, if it was assigned on the heap.

Definition at line 407 of file TMatrixT.cxx.

◆ Determinant() [1/2]

template<class Element >
Double_t TMatrixT< Element >::Determinant ( ) const
virtual

Return the matrix determinant.

Reimplemented from TMatrixTBase< Element >.

Definition at line 1361 of file TMatrixT.cxx.

◆ Determinant() [2/2]

template<class Element >
void TMatrixT< Element >::Determinant ( Double_t &  d1,
Double_t &  d2 
) const
virtual

Return the matrix determinant as d1,d2 where det = d1*TMath::Power(2.0,d2)

Reimplemented from TMatrixTBase< Element >.

Definition at line 1374 of file TMatrixT.cxx.

◆ EigenVectors()

template<class Element>
const TMatrixT< Element > TMatrixT< Element >::EigenVectors ( TVectorT< Element > &  eigenValues) const

Return a matrix containing the eigen-vectors ordered by descending values of Re^2+Im^2 of the complex eigen-values .

If the matrix is asymmetric, only the real part of the eigen-values is returned . For full functionality use TMatrixDEigen .

Definition at line 2307 of file TMatrixT.cxx.

◆ GetColIndexArray() [1/2]

template<class Element>
virtual const Int_t* TMatrixT< Element >::GetColIndexArray ( ) const
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 114 of file TMatrixT.h.

◆ GetColIndexArray() [2/2]

template<class Element>
virtual Int_t* TMatrixT< Element >::GetColIndexArray ( )
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 115 of file TMatrixT.h.

◆ GetMatrixArray() [1/2]

template<class Element >
const Element * TMatrixT< Element >::GetMatrixArray ( ) const
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 222 of file TMatrixT.h.

◆ GetMatrixArray() [2/2]

template<class Element >
Element * TMatrixT< Element >::GetMatrixArray ( )
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 223 of file TMatrixT.h.

◆ GetRowIndexArray() [1/2]

template<class Element>
virtual const Int_t* TMatrixT< Element >::GetRowIndexArray ( ) const
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 112 of file TMatrixT.h.

◆ GetRowIndexArray() [2/2]

template<class Element>
virtual Int_t* TMatrixT< Element >::GetRowIndexArray ( )
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 113 of file TMatrixT.h.

◆ GetSub() [1/2]

template<class Element>
TMatrixTBase< Element > & TMatrixT< Element >::GetSub ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
TMatrixTBase< Element > &  target,
Option_t *  option = "S" 
) const
virtual

Get submatrix [row_lwb..row_upb] x [col_lwb..col_upb]; The indexing range of the returned matrix depends on the argument option:

option == "S" : return [0..row_upb-row_lwb][0..col_upb-col_lwb] (default) else : return [row_lwb..row_upb][col_lwb..col_upb]

Implements TMatrixTBase< Element >.

Definition at line 1089 of file TMatrixT.cxx.

◆ GetSub() [2/2]

template<class Element>
TMatrixT< Element > TMatrixT< Element >::GetSub ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
Option_t *  option = "S" 
) const
inline

Definition at line 242 of file TMatrixT.h.

◆ Invert()

template<class Element >
TMatrixT< Element > & TMatrixT< Element >::Invert ( Double_t *  det = 0)

Invert the matrix and calculate its determinant.

Definition at line 1396 of file TMatrixT.cxx.

◆ InvertFast()

template<class Element >
TMatrixT< Element > & TMatrixT< Element >::InvertFast ( Double_t *  det = 0)

Invert the matrix and calculate its determinant, however upto (6x6) a fast Cramer inversion is used .

Definition at line 1410 of file TMatrixT.cxx.

◆ Memcpy_m()

template<class Element>
Int_t TMatrixT< Element >::Memcpy_m ( Element *  newp,
const Element *  oldp,
Int_t  copySize,
Int_t  newSize,
Int_t  oldSize 
)
protected

Copy copySize doubles from *oldp to *newp .

However take care of the situation where both pointers are assigned to the same stack space

Definition at line 439 of file TMatrixT.cxx.

◆ Minus() [1/3]

template<class Element>
void TMatrixT< Element >::Minus ( const TMatrixT< Element > &  a,
const TMatrixT< Element > &  b 
)

General matrix summation. Create a matrix C such that C = A - B.

Definition at line 580 of file TMatrixT.cxx.

◆ Minus() [2/3]

template<class Element>
void TMatrixT< Element >::Minus ( const TMatrixT< Element > &  a,
const TMatrixTSym< Element > &  b 
)

General matrix summation. Create a matrix C such that C = A - B.

Definition at line 614 of file TMatrixT.cxx.

◆ Minus() [3/3]

template<class Element>
void TMatrixT< Element >::Minus ( const TMatrixTSym< Element > &  a,
const TMatrixT< Element > &  b 
)
inline

Definition at line 93 of file TMatrixT.h.

◆ Mult() [1/4]

template<class Element>
void TMatrixT< Element >::Mult ( const TMatrixT< Element > &  a,
const TMatrixT< Element > &  b 
)

General matrix multiplication. Create a matrix C such that C = A * B.

Definition at line 648 of file TMatrixT.cxx.

◆ Mult() [2/4]

template<class Element>
void TMatrixT< Element >::Mult ( const TMatrixT< Element > &  a,
const TMatrixTSym< Element > &  b 
)

Matrix multiplication, with A general and B symmetric.

Create a matrix C such that C = A * B.

Definition at line 749 of file TMatrixT.cxx.

◆ Mult() [3/4]

template<class Element>
void TMatrixT< Element >::Mult ( const TMatrixTSym< Element > &  a,
const TMatrixT< Element > &  b 
)

Matrix multiplication, with A symmetric and B general.

Create a matrix C such that C = A * B.

Definition at line 697 of file TMatrixT.cxx.

◆ Mult() [4/4]

template<class Element>
void TMatrixT< Element >::Mult ( const TMatrixTSym< Element > &  a,
const TMatrixTSym< Element > &  b 
)

Matrix multiplication, with A symmetric and B symmetric.

(Actually copied for the moment routine for B general) Create a matrix C such that C = A * B.

Definition at line 801 of file TMatrixT.cxx.

◆ MultT() [1/4]

template<class Element>
void TMatrixT< Element >::MultT ( const TMatrixT< Element > &  a,
const TMatrixT< Element > &  b 
)

General matrix multiplication. Create a matrix C such that C = A * B^T.

Definition at line 951 of file TMatrixT.cxx.

◆ MultT() [2/4]

template<class Element>
void TMatrixT< Element >::MultT ( const TMatrixT< Element > &  a,
const TMatrixTSym< Element > &  b 
)
inline

Definition at line 106 of file TMatrixT.h.

◆ MultT() [3/4]

template<class Element>
void TMatrixT< Element >::MultT ( const TMatrixTSym< Element > &  a,
const TMatrixT< Element > &  b 
)

Matrix multiplication, with A symmetric and B general.

Create a matrix C such that C = A * B^T.

Definition at line 1003 of file TMatrixT.cxx.

◆ MultT() [4/4]

template<class Element>
void TMatrixT< Element >::MultT ( const TMatrixTSym< Element > &  a,
const TMatrixTSym< Element > &  b 
)
inline

Definition at line 108 of file TMatrixT.h.

◆ New_m()

template<class Element >
Element * TMatrixT< Element >::New_m ( Int_t  size)
protected

Return data pointer .

if requested size <= kSizeMax, assign pointer to the stack space

Definition at line 421 of file TMatrixT.cxx.

◆ NormByColumn()

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::NormByColumn ( const TVectorT< Element > &  v,
Option_t *  option = "D" 
)

Multiply/divide matrix columns by a vector: option: "D" : b(i,j) = a(i,j)/v(i) i = 0,fNrows-1 (default) else : b(i,j) = a(i,j)*v(i)

Definition at line 1642 of file TMatrixT.cxx.

◆ NormByRow()

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::NormByRow ( const TVectorT< Element > &  v,
Option_t *  option = "D" 
)

Multiply/divide matrix rows with a vector: option: "D" : b(i,j) = a(i,j)/v(j) i = 0,fNcols-1 (default) else : b(i,j) = a(i,j)*v(j)

Definition at line 1689 of file TMatrixT.cxx.

◆ operator()() [1/2]

template<class Element >
Element TMatrixT< Element >::operator() ( Int_t  rown,
Int_t  coln 
) const
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 250 of file TMatrixT.h.

◆ operator()() [2/2]

template<class Element >
Element & TMatrixT< Element >::operator() ( Int_t  rown,
Int_t  coln 
)
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 267 of file TMatrixT.h.

◆ operator*=() [1/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( Element  val)

Multiply every element of the matrix with val.

Definition at line 1874 of file TMatrixT.cxx.

◆ operator*=() [2/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( const TMatrixT< Element > &  source)

Compute target = target * source inplace.

Strictly speaking, it can't be done inplace, though only the row of the target matrix needs to be saved. "Inplace" multiplication is only allowed when the 'source' matrix is square.

Definition at line 1972 of file TMatrixT.cxx.

◆ operator*=() [3/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( const TMatrixTSym< Element > &  source)

Compute target = target * source inplace.

Strictly speaking, it can't be done inplace, though only the row of the target matrix needs to be saved.

Definition at line 2035 of file TMatrixT.cxx.

◆ operator*=() [4/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( const TMatrixTDiag_const< Element > &  diag)

Multiply a matrix row by the diagonal of another matrix matrix(i,j) *= diag(j), j=0,fNcols-1.

Definition at line 2097 of file TMatrixT.cxx.

◆ operator*=() [5/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( const TMatrixTRow_const< Element > &  row)

Multiply a matrix by the row of another matrix matrix(i,j) *= another(k,j) for fixed k.

Definition at line 2234 of file TMatrixT.cxx.

◆ operator*=() [6/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator*= ( const TMatrixTColumn_const< Element > &  col)

Multiply a matrix by the column of another matrix matrix(i,j) *= another(i,k) for fixed k.

Definition at line 2162 of file TMatrixT.cxx.

◆ operator+=() [1/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator+= ( Element  val)

Add val to every element of the matrix.

Definition at line 1842 of file TMatrixT.cxx.

◆ operator+=() [2/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator+= ( const TMatrixT< Element > &  source)

Add the source matrix.

Definition at line 1890 of file TMatrixT.cxx.

◆ operator+=() [3/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator+= ( const TMatrixTSym< Element > &  source)

Add the source matrix.

Definition at line 1910 of file TMatrixT.cxx.

◆ operator-=() [1/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator-= ( Element  val)

Subtract val from every element of the matrix.

Definition at line 1858 of file TMatrixT.cxx.

◆ operator-=() [2/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator-= ( const TMatrixT< Element > &  source)

Subtract the source matrix.

Definition at line 1930 of file TMatrixT.cxx.

◆ operator-=() [3/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator-= ( const TMatrixTSym< Element > &  source)

Subtract the source matrix.

Definition at line 1950 of file TMatrixT.cxx.

◆ operator/=() [1/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator/= ( const TMatrixTDiag_const< Element > &  diag)

Divide a matrix row by the diagonal of another matrix matrix(i,j) /= diag(j)

Definition at line 2127 of file TMatrixT.cxx.

◆ operator/=() [2/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator/= ( const TMatrixTRow_const< Element > &  row)

Divide a matrix by the row of another matrix matrix(i,j) /= another(k,j) for fixed k.

Definition at line 2268 of file TMatrixT.cxx.

◆ operator/=() [3/3]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator/= ( const TMatrixTColumn_const< Element > &  col)

Divide a matrix by the column of another matrix matrix(i,j) /= another(i,k) for fixed k.

Definition at line 2195 of file TMatrixT.cxx.

◆ operator=() [1/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator= ( const TMatrixT< Element > &  source)

Assignment operator.

Definition at line 1732 of file TMatrixT.cxx.

◆ operator=() [2/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator= ( const TMatrixTSym< Element > &  source)

Assignment operator.

Definition at line 1751 of file TMatrixT.cxx.

◆ operator=() [3/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator= ( const TMatrixTSparse< Element > &  source)

Assignment operator.

Definition at line 1770 of file TMatrixT.cxx.

◆ operator=() [4/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator= ( const TMatrixTLazy< Element > &  source)

Assignment operator.

Definition at line 1805 of file TMatrixT.cxx.

◆ operator=() [5/6]

template<class Element>
template<class Element2 >
TMatrixT<Element>& TMatrixT< Element >::operator= ( const TMatrixT< Element2 > &  source)
inline

Definition at line 170 of file TMatrixT.h.

◆ operator=() [6/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::operator= ( Element  val)

Assign val to every element of the matrix.

Definition at line 1826 of file TMatrixT.cxx.

◆ operator[]() [1/2]

template<class Element>
const TMatrixTRow_const<Element> TMatrixT< Element >::operator[] ( Int_t  rown) const
inline

Definition at line 163 of file TMatrixT.h.

◆ operator[]() [2/2]

template<class Element>
TMatrixTRow<Element> TMatrixT< Element >::operator[] ( Int_t  rown)
inline

Definition at line 164 of file TMatrixT.h.

◆ Plus() [1/3]

template<class Element>
void TMatrixT< Element >::Plus ( const TMatrixT< Element > &  a,
const TMatrixT< Element > &  b 
)

General matrix summation. Create a matrix C such that C = A + B.

Definition at line 512 of file TMatrixT.cxx.

◆ Plus() [2/3]

template<class Element>
void TMatrixT< Element >::Plus ( const TMatrixT< Element > &  a,
const TMatrixTSym< Element > &  b 
)

General matrix summation. Create a matrix C such that C = A + B.

Definition at line 546 of file TMatrixT.cxx.

◆ Plus() [3/3]

template<class Element>
void TMatrixT< Element >::Plus ( const TMatrixTSym< Element > &  a,
const TMatrixT< Element > &  b 
)
inline

Definition at line 89 of file TMatrixT.h.

◆ Rank1Update() [1/2]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Rank1Update ( const TVectorT< Element > &  v,
Element  alpha = 1.0 
)

Perform a rank 1 operation on matrix A: A += alpha * v * v^T.

Definition at line 1539 of file TMatrixT.cxx.

◆ Rank1Update() [2/2]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Rank1Update ( const TVectorT< Element > &  v1,
const TVectorT< Element > &  v2,
Element  alpha = 1.0 
)

Perform a rank 1 operation on matrix A: A += alpha * v1 * v2^T.

Definition at line 1567 of file TMatrixT.cxx.

◆ ResizeTo() [1/3]

template<class Element >
TMatrixTBase< Element > & TMatrixT< Element >::ResizeTo ( Int_t  nrows,
Int_t  ncols,
Int_t  = -1 
)
virtual

Set size of the matrix to nrows x ncols New dynamic elements are created, the overlapping part of the old ones are copied to the new structures, then the old elements are deleted.

Implements TMatrixTBase< Element >.

Definition at line 1210 of file TMatrixT.cxx.

◆ ResizeTo() [2/3]

template<class Element >
TMatrixTBase< Element > & TMatrixT< Element >::ResizeTo ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
Int_t  = -1 
)
virtual

Set size of the matrix to [row_lwb:row_upb] x [col_lwb:col_upb] New dynamic elemenst are created, the overlapping part of the old ones are copied to the new structures, then the old elements are deleted.

Implements TMatrixTBase< Element >.

Definition at line 1275 of file TMatrixT.cxx.

◆ ResizeTo() [3/3]

template<class Element>
TMatrixTBase<Element>& TMatrixT< Element >::ResizeTo ( const TMatrixT< Element > &  m)
inline

Definition at line 140 of file TMatrixT.h.

◆ SetColIndexArray()

template<class Element>
virtual TMatrixTBase<Element>& TMatrixT< Element >::SetColIndexArray ( Int_t *  )
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 118 of file TMatrixT.h.

◆ SetRowIndexArray()

template<class Element>
virtual TMatrixTBase<Element>& TMatrixT< Element >::SetRowIndexArray ( Int_t *  )
inlinevirtual

Implements TMatrixTBase< Element >.

Definition at line 117 of file TMatrixT.h.

◆ SetSub()

template<class Element>
TMatrixTBase< Element > & TMatrixT< Element >::SetSub ( Int_t  row_lwb,
Int_t  col_lwb,
const TMatrixTBase< Element > &  source 
)
virtual

Insert matrix source starting at [row_lwb][col_lwb], thereby overwriting the part [row_lwb..row_lwb+nrows_source][col_lwb..col_lwb+ncols_source];.

Implements TMatrixTBase< Element >.

Definition at line 1156 of file TMatrixT.cxx.

◆ Similarity()

template<class Element>
Element TMatrixT< Element >::Similarity ( const TVectorT< Element > &  v) const

Calculate scalar v * (*this) * v^T.

Definition at line 1601 of file TMatrixT.cxx.

◆ T()

template<class Element>
TMatrixT<Element>& TMatrixT< Element >::T ( )
inline

Definition at line 150 of file TMatrixT.h.

◆ TMult() [1/4]

template<class Element>
void TMatrixT< Element >::TMult ( const TMatrixT< Element > &  a,
const TMatrixT< Element > &  b 
)

Create a matrix C such that C = A' * B.

In other words, c[i,j] = SUM{ a[k,i] * b[k,j] }.

Definition at line 852 of file TMatrixT.cxx.

◆ TMult() [2/4]

template<class Element>
void TMatrixT< Element >::TMult ( const TMatrixT< Element > &  a,
const TMatrixTSym< Element > &  b 
)

Create a matrix C such that C = A' * B.

In other words, c[i,j] = SUM{ a[k,i] * b[k,j] }.

Definition at line 902 of file TMatrixT.cxx.

◆ TMult() [3/4]

template<class Element>
void TMatrixT< Element >::TMult ( const TMatrixTSym< Element > &  a,
const TMatrixT< Element > &  b 
)
inline

Definition at line 102 of file TMatrixT.h.

◆ TMult() [4/4]

template<class Element>
void TMatrixT< Element >::TMult ( const TMatrixTSym< Element > &  a,
const TMatrixTSym< Element > &  b 
)
inline

Definition at line 103 of file TMatrixT.h.

◆ Transpose()

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Transpose ( const TMatrixT< Element > &  source)

Transpose matrix source.

Definition at line 1469 of file TMatrixT.cxx.

◆ Use() [1/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Use ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
Element *  data 
)

Use the array data to fill the matrix ([row_lwb..row_upb] x [col_lwb..col_upb])

Definition at line 1053 of file TMatrixT.cxx.

◆ Use() [2/6]

template<class Element>
const TMatrixT<Element>& TMatrixT< Element >::Use ( Int_t  row_lwb,
Int_t  row_upb,
Int_t  col_lwb,
Int_t  col_upb,
const Element *  data 
) const
inline

Definition at line 125 of file TMatrixT.h.

◆ Use() [3/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Use ( Int_t  nrows,
Int_t  ncols,
Element *  data 
)
inline

Definition at line 225 of file TMatrixT.h.

◆ Use() [4/6]

template<class Element>
const TMatrixT< Element > & TMatrixT< Element >::Use ( Int_t  nrows,
Int_t  ncols,
const Element *  data 
) const
inline

Definition at line 227 of file TMatrixT.h.

◆ Use() [5/6]

template<class Element>
TMatrixT< Element > & TMatrixT< Element >::Use ( TMatrixT< Element > &  a)
inline

Definition at line 229 of file TMatrixT.h.

◆ Use() [6/6]

template<class Element>
const TMatrixT< Element > & TMatrixT< Element >::Use ( const TMatrixT< Element > &  a) const
inline

Definition at line 235 of file TMatrixT.h.

Member Data Documentation

◆ fDataStack

template<class Element>
Element TMatrixT< Element >::fDataStack[TMatrixTBase< Element >::kSizeMax]
protected

Definition at line 43 of file TMatrixT.h.

◆ fElements

template<class Element>
Element* TMatrixT< Element >::fElements
protected

data container

Definition at line 44 of file TMatrixT.h.

Libraries for TMatrixT< Element >:
[legend]

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