class ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >


    SMatrix: a generic fixed size D1 x D2 Matrix class.
    The class is template on the scalar type, on the matrix sizes:
    D1 = number of rows and D2 = number of columns
    amd on the representation storage type.
    By default the representation is MatRepStd<T,D1,D2> (standard D1xD2 of type T),
    but it can be of type MatRepSym<T,D> for symmetric matrices DxD, where the storage is only
    D*(D+1)/2.

    See \ref SMatrixDoc.

    Original author is Thorsten Glebe
    HERA-B Collaboration, MPI Heidelberg (Germany)

    @ingroup SMatrixSVector

    @authors T. Glebe, L. Moneta and J. Palacios


 SMatrix: column-wise storage

Function Members (Methods)

public:
~SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >()
floatapply(unsigned int i) const
const float*Array() const
float*Array()
const float&At(unsigned int i, unsigned int j) const
float&At(unsigned int i, unsigned int j)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >::iteratorbegin()
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >::const_iteratorbegin() const
ROOT::Math::SVector<float,2>Col(unsigned int thecol) const
boolDet(float& det)
boolDet2(float& det) const
ROOT::Math::SVector<float,2>Diagonal() const
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >::iteratorend()
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >::const_iteratorend() const
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >Inverse(int& ifail) const
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >InverseFast(int& ifail) const
boolInvert()
boolInvertFast()
boolIsInUse(const float* p) const
ROOT::Math::SVector<float,2*(2+1)/2>LowerBlock() const
booloperator!=(const float& rhs) const
booloperator!=(const ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >& rhs) const
const float&operator()(unsigned int i, unsigned int j) const
float&operator()(unsigned int i, unsigned int j)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator*=(const float& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator+=(const float& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator-=(const float& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator/=(const float& rhs)
booloperator<(const float& rhs) const
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator=(ROOT::Math::SMatrixIdentity)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator=(const float& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&operator=(const ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >&)
booloperator==(const float& rhs) const
booloperator>(const float& rhs) const
ostream&Print(ostream& os) const
ROOT::Math::SVector<float,2>Row(unsigned int therow) const
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >()
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >(ROOT::Math::SMatrixIdentity)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >(const ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >(const float& rhs)
ROOT::Math::SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >SMatrix<float,2,2,ROOT::Math::MatRepSym<float,2> >(const ROOT::Math::SVector<float,2*(2+1)/2>& v, bool lower = true)
floatTrace() const
ROOT::Math::SVector<float,2*(2+1)/2>UpperBlock() const

Data Members

public:
enum { kRows
kCols
kSize
};
public:
ROOT::Math::MatRepSym<float,2>fRep

Class Charts

Inheritance Inherited Members Includes Libraries
Class Charts

Function documentation

SMatrix<T,D1,D2,R>& operator=(const M& rhs)
       Assign from another compatible matrix.
       Possible Symmetirc to general but NOT vice-versa
SMatrix<T,D1,D2,R>& operator=(const Expr<A,T,D1,D2,R2>& rhs)
       Assign from a matrix expression
SMatrix<T,D1,D2,R> & operator=(SMatrixIdentity )
      Assign from an identity matrix
T apply(unsigned int i) const
 @name --- Access functions --- 
 access the parse tree with the index starting from zero and
       following the C convention for the order in accessing
       the matrix elements.
       Same convention for general and symmetric matrices.
const T* Array()
 return read-only pointer to internal array
T* Array()
iterator begin()
 @name --- STL-like interface ---
       The iterators access the matrix element in the order how they are
       stored in memory. The C (row-major) convention is used, and in the
       case of symmetric matrices the iterator spans only the lower diagonal
       block. For example for a symmetric 3x3 matrices the order of the 6
       elements \f${a_0,...a_5}\f$ is:
       \f[
       M = \left( \begin{array}{ccc}
       a_0 & a_1 & a_3  \\
       a_1 & a_2  & a_4  \\
       a_3 & a_4 & a_5   \end{array} \right)
       \f]

 STL iterator interface.
iterator end()
 STL iterator interface.
const_iterator begin()
 STL const_iterator interface.
const_iterator end()
 STL const_iterator interface.
bool operator==(const float& rhs) const
 @name --- Operators --- 
 element wise comparison
bool operator!=(const T& rhs)
 element wise comparison
bool operator!=(const SMatrix<T,D1,D2,R>& rhs)
 element wise comparison
bool operator>(const float& rhs) const
 element wise comparison
bool operator<(const float& rhs) const
 element wise comparison
const T& operator()(unsigned int i, unsigned int j)
      read only access to matrix element, with indices starting from 0
T& operator()(unsigned int i, unsigned int j)
const T& At(unsigned int i, unsigned int j)
      read only access to matrix element, with indices starting from 0.
      Fuction will check index values and it will assert if they are wrong
T& At(unsigned int i, unsigned int j)
SMatrix<T,D1,D2,R>& operator+=(const float& rhs)
      addition with a scalar
SMatrix<T,D1,D2,R>& operator-=(const float& rhs)
      subtraction with a scalar
SMatrix<T,D1,D2,R>& operator*=(const float& rhs)
      multiplication with a scalar
SMatrix<T,D1,D2,R>& operator/=(const float& rhs)
      division with a scalar
bool Invert()
 @name --- Linear Algebra Functions --- 

      Invert a square Matrix ( this method changes the current matrix).
      Return true if inversion is successfull.
      The method used for general square matrices is the LU factorization taken from Dinv routine
      from the CERNLIB (written in C++ from CLHEP authors)
      In case of symmetric matrices Bunch-Kaufman diagonal pivoting method is used
      (The implementation is the one written by the CLHEP authors)
SMatrix<T,D1,D2,R> Inverse(int& ifail) const
      Invert a square Matrix and  returns a new matrix. In case the inversion fails
      the current matrix is returned.
      \param ifail . ifail will be set to 0 when inversion is successfull.
      See ROOT::Math::SMatrix::Invert for the inversion algorithm
bool InvertFast()
      Fast Invertion of a square Matrix ( this method changes the current matrix).
      Return true if inversion is successfull.
      The method used is based on direct inversion using the Cramer rule for
      matrices upto 5x5. Afterwards the same defult algorithm of Invert() is used.
      Note that this method is faster but can suffer from much larger numerical accuracy
      when the condition of the matrix is large
SMatrix<T,D1,D2,R> InverseFast(int& ifail) const
      Invert a square Matrix and  returns a new matrix. In case the inversion fails
      the current matrix is returned.
      \param ifail . ifail will be set to 0 when inversion is successfull.
      See ROOT::Math::SMatrix::InverseFast for the inversion algorithm
bool Det(float& det)
      determinant of square Matrix via Dfact.
      Return true when the calculation is successfull.
      \param det will contain the calculated determinant value
      \b Note: this will destroy the contents of the Matrix!
bool Det2(float& det) const
      determinant of square Matrix via Dfact.
      Return true when the calculation is successfull.
      \param det will contain the calculated determinant value
      \b Note: this will preserve the content of the Matrix!
SVector<T,D2> Row(unsigned int therow) const
      return a full Matrix row as a vector (copy the content in a new vector)
SVector<T,D1> Col(unsigned int thecol) const
      return a full Matrix column as a vector (copy the content in a new vector)
SVector<T,D1> Diagonal() const
      return diagonal elements of a matrix as a Vector.
      It works only for squared matrices D1 == D2, otherwise it will produce a compile error
T Trace() const
      return the trace of a matrix
      Sum of the diagonal elements
SubVector UpperBlock() const
SubVector LowerBlock() const
bool IsInUse(const float* p) const
 @name --- Other Functions --- 

       Function to check if a matrix is sharing same memory location of the passed pointer
       This function is used by the expression templates to avoid the alias problem during
       expression evaluation. When  the matrix is in use, for example in operations
       like A = B * A, a temporary object storing the intermediate result is automatically
       created when evaluating the expression.
std::ostream& Print(ostream& os) const
 submatrices
 Print: used by operator<<()

Author: T. Glebe, L. Moneta, J. Palacios 2005
Last change: root/smatrix:$Id: SMatrix.h 26866 2008-12-12 10:50:07Z moneta $
Last generated: 2008-12-12 12:11
This page has been automatically generated. If you have any comments or suggestions about the page layout send a mail to ROOT support, or contact the developers with any questions or problems regarding ROOT.