ROOT logo
AliRoot » STEER » STEER » AliSymMatrix

class AliSymMatrix: public AliMatrixSq

Function Members (Methods)

public:
AliSymMatrix()
AliSymMatrix(Int_t size)
AliSymMatrix(const AliSymMatrix& mat)
virtual~AliSymMatrix()
virtual TMatrixTBase<double>&TMatrixTBase<double>::Abs()
voidTObject::AbstractMethod(const char* method) const
voidAddRows(int nrows = 1)
virtual voidAddToRow(Int_t r, Double_t* valc, Int_t* indc, Int_t n)
virtual voidAliMatrixSq::Allocate(Int_t, Int_t, Int_t, Int_t, Int_t, Int_t)
virtual voidTObject::AppendPad(Option_t* option = "")
virtual TMatrixTBase<double>&TMatrixTBase<double>::Apply(const TElementActionT<double>& action)
virtual TMatrixTBase<double>&TMatrixTBase<double>::Apply(const TElementPosActionT<double>& action)
virtual voidTObject::Browse(TBrowser* b)
static TClass*Class()
virtual const char*TObject::ClassName() const
virtual voidClear(Option_t* option = "")
virtual TObject*TObject::Clone(const char* newname = "") const
virtual doubleTMatrixTBase<double>::ColNorm() const
virtual Int_tTObject::Compare(const TObject* obj) const
virtual voidTObject::Copy(TObject& object) const
AliSymMatrix*DecomposeChol()
virtual voidTObject::Delete(Option_t* option = "")MENU
virtual Double_tTMatrixTBase<double>::Determinant() const
virtual voidTMatrixTBase<double>::Determinant(Double_t& d1, Double_t& d2) const
virtual Double_tDiagElem(Int_t r) const
virtual Double_t&DiagElem(Int_t r)
virtual Int_tTObject::DistancetoPrimitive(Int_t px, Int_t py)
virtual voidTMatrixTBase<double>::Draw(Option_t* option = "")MENU
virtual voidTObject::DrawClass() constMENU
virtual TObject*TObject::DrawClone(Option_t* option = "") constMENU
virtual voidTObject::Dump() constMENU
virtual doubleTMatrixTBase<double>::E2Norm() const
virtual voidTObject::Error(const char* method, const char* msgfmt) const
virtual voidTObject::Execute(const char* method, const char* params, Int_t* error = 0)
virtual voidTObject::Execute(TMethod* method, TObjArray* params, Int_t* error = 0)
virtual voidTObject::ExecuteEvent(Int_t event, Int_t px, Int_t py)
virtual voidTMatrixTBase<double>::ExtractRow(Int_t row, Int_t col, double* v, Int_t n = -1) const
virtual voidTObject::Fatal(const char* method, const char* msgfmt) const
virtual TObject*TObject::FindObject(const char* name) const
virtual TObject*TObject::FindObject(const TObject* obj) const
virtual const Int_t*GetColIndexArray() const
virtual Int_t*GetColIndexArray()
Int_tTMatrixTBase<double>::GetColLwb() const
Int_tTMatrixTBase<double>::GetColUpb() const
virtual Float_tGetDensity() const
virtual Option_t*TObject::GetDrawOption() const
static Long_tTObject::GetDtorOnly()
virtual const char*TObject::GetIconName() const
virtual voidTMatrixTBase<double>::GetMatrix2Array(double* data, Option_t* option = "") const
virtual const Double_t*GetMatrixArray() const
virtual Double_t*GetMatrixArray()
virtual const char*TObject::GetName() const
Int_tTMatrixTBase<double>::GetNcols() const
Int_tTMatrixTBase<double>::GetNoElements() const
Int_tTMatrixTBase<double>::GetNrows() const
virtual char*TObject::GetObjectInfo(Int_t px, Int_t py) const
static Bool_tTObject::GetObjectStat()
virtual Option_t*TObject::GetOption() const
Double_t*GetRow(Int_t r)
virtual const Int_t*GetRowIndexArray() const
virtual Int_t*GetRowIndexArray()
Int_tTMatrixTBase<double>::GetRowLwb() const
Int_tTMatrixTBase<double>::GetRowUpb() const
virtual Int_tGetSize() const
Int_tGetSizeAdded() const
Int_tGetSizeBooked() const
Int_tGetSizeUsed() const
virtual TMatrixDBase&GetSub(Int_t, Int_t, Int_t, Int_t, TMatrixDBase&, Option_t*) const
virtual const char*TObject::GetTitle() const
doubleTMatrixTBase<double>::GetTol() const
virtual UInt_tTObject::GetUniqueID() const
virtual Bool_tTObject::HandleTimer(TTimer* timer)
virtual ULong_tTObject::Hash() const
virtual voidTObject::Info(const char* method, const char* msgfmt) const
virtual Bool_tTObject::InheritsFrom(const char* classname) const
virtual Bool_tTObject::InheritsFrom(const TClass* cl) const
virtual TMatrixTBase<double>&TMatrixTBase<double>::InsertRow(Int_t row, Int_t col, const double* v, Int_t n = -1)
virtual voidTObject::Inspect() constMENU
voidTMatrixTBase<double>::Invalidate()
voidTObject::InvertBit(UInt_t f)
Bool_tInvertChol()
voidInvertChol(AliSymMatrix* mchol)
virtual TClass*IsA() const
virtual Bool_tTObject::IsEqual(const TObject* obj) const
virtual Bool_tTObject::IsFolder() const
Bool_tTObject::IsOnHeap() const
Bool_tTMatrixTBase<double>::IsOwner() const
virtual Bool_tTObject::IsSortable() const
virtual Bool_tAliMatrixSq::IsSymmetric() const
Bool_tTMatrixTBase<double>::IsValid() const
static Bool_tAliMatrixSq::IsZero(Double_t x, Double_t thresh = 1e-64)
Bool_tTObject::IsZombie() const
virtual voidTObject::ls(Option_t* option = "") const
voidTMatrixTBase<double>::MakeValid()
virtual doubleTMatrixTBase<double>::Max() const
voidTObject::MayNotUse(const char* method) const
virtual doubleTMatrixTBase<double>::Min() const
virtual voidMultiplyByVec(const Double_t* vecIn, Double_t* vecOut) const
virtual voidMultiplyByVec(const TVectorD& vecIn, TVectorD& vecOut) const
virtual Int_tTMatrixTBase<double>::NonZeros() const
doubleTMatrixTBase<double>::Norm1() const
virtual TMatrixTBase<double>&TMatrixTBase<double>::NormByDiag(const TVectorT<double>& v, Option_t* option = "D")
doubleTMatrixTBase<double>::NormInf() const
virtual Bool_tTObject::Notify()
voidTObject::Obsolete(const char* method, const char* asOfVers, const char* removedFromVers) const
static voidTObject::operator delete(void* ptr)
static voidTObject::operator delete(void* ptr, void* vp)
static voidTObject::operator delete[](void* ptr)
static voidTObject::operator delete[](void* ptr, void* vp)
void*TObject::operator new(size_t sz)
void*TObject::operator new(size_t sz, void* vp)
void*TObject::operator new[](size_t sz)
void*TObject::operator new[](size_t sz, void* vp)
Bool_tTMatrixTBase<double>::operator!=(double val) const
virtual Double_toperator()(Int_t row, Int_t col) const
virtual Double_t&operator()(Int_t row, Int_t col)
AliSymMatrix&operator+=(const AliSymMatrix& src)
Bool_tTMatrixTBase<double>::operator<(double val) const
Bool_tTMatrixTBase<double>::operator<=(double val) const
AliSymMatrix&operator=(const AliSymMatrix& src)
Bool_tTMatrixTBase<double>::operator==(double val) const
Bool_tTMatrixTBase<double>::operator>(double val) const
Bool_tTMatrixTBase<double>::operator>=(double val) const
virtual voidTObject::Paint(Option_t* option = "")
virtual voidTObject::Pop()
virtual voidPrint(Option_t* option = "") const
virtual voidAliMatrixSq::PrintCOO() const
virtual Double_tAliMatrixSq::Query(Int_t rown, Int_t coln) const
virtual Double_tAliMatrixSq::QueryDiag(Int_t rc) const
virtual TMatrixTBase<double>&TMatrixTBase<double>::Randomize(double alpha, double beta, Double_t& seed)
virtual Int_tTObject::Read(const char* name)
virtual voidTObject::RecursiveRemove(TObject* obj)
virtual voidReset()
voidTObject::ResetBit(UInt_t f)
virtual TMatrixDBase&ResizeTo(Int_t, Int_t, Int_t)
virtual TMatrixDBase&ResizeTo(Int_t, Int_t, Int_t, Int_t, Int_t)
virtual doubleTMatrixTBase<double>::RowNorm() const
virtual voidTObject::SaveAs(const char* filename = "", Option_t* option = "") constMENU
virtual voidTObject::SavePrimitive(basic_ostream<char,char_traits<char> >& out, Option_t* option = "")
voidScale(Double_t coeff)
voidTObject::SetBit(UInt_t f)
voidTObject::SetBit(UInt_t f, Bool_t set)
virtual TMatrixDBase&SetColIndexArray(Int_t*)
virtual voidTObject::SetDrawOption(Option_t* option = "")MENU
static voidTObject::SetDtorOnly(void* obj)
virtual TMatrixTBase<double>&TMatrixTBase<double>::SetMatrixArray(const double* data, Option_t* option = "")
static voidTObject::SetObjectStat(Bool_t stat)
virtual TMatrixDBase&SetRowIndexArray(Int_t*)
voidSetSizeUsed(Int_t sz)
virtual TMatrixDBase&SetSub(Int_t, Int_t, const TMatrixDBase&)
voidAliMatrixSq::SetSymmetric(Bool_t v = kTRUE)
doubleTMatrixTBase<double>::SetTol(double newTol)
virtual voidTObject::SetUniqueID(UInt_t uid)
virtual TMatrixTBase<double>&TMatrixTBase<double>::Shift(Int_t row_shift, Int_t col_shift)
virtual voidShowMembers(TMemberInspector&)
Bool_tSolveChol(Double_t* brhs, Bool_t invert = kFALSE)
Bool_tSolveChol(TVectorD& brhs, Bool_t invert = kFALSE)
Bool_tSolveChol(Double_t* brhs, Double_t* bsol, Bool_t invert = kFALSE)
Bool_tSolveChol(const TVectorD& brhs, TVectorD& bsol, Bool_t invert = kFALSE)
intSolveSpmInv(double* vecB, Bool_t stabilize = kTRUE)
virtual TMatrixTBase<double>&TMatrixTBase<double>::Sqr()
virtual TMatrixTBase<double>&TMatrixTBase<double>::Sqrt()
virtual voidStreamer(TBuffer&)
voidStreamerNVirtual(TBuffer& ClassDef_StreamerNVirtual_b)
virtual doubleTMatrixTBase<double>::Sum() const
virtual voidTObject::SysError(const char* method, const char* msgfmt) const
Bool_tTObject::TestBit(UInt_t f) const
Int_tTObject::TestBits(UInt_t f) const
virtual TMatrixTBase<double>&TMatrixTBase<double>::UnitMatrix()
virtual voidTObject::UseCurrentStyle()
virtual voidTObject::Warning(const char* method, const char* msgfmt) const
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0)
virtual Int_tTObject::Write(const char* name = 0, Int_t option = 0, Int_t bufsize = 0) const
virtual TMatrixTBase<double>&TMatrixTBase<double>::Zero()
protected:
virtual voidTObject::DoError(int level, const char* location, const char* fmt, va_list va) const
static voidTMatrixTBase<double>::DoubleLexSort(Int_t n, Int_t* first, Int_t* second, double* data)
Double_tGetEl(Int_t row, Int_t col) const
virtual Int_tGetIndex(Int_t row, Int_t col) const
static voidTMatrixTBase<double>::IndexedLexSort(Int_t n, Int_t* first, Int_t swapFirst, Int_t* second, Int_t swapSecond, Int_t* index)
voidTObject::MakeZombie()
voidSetEl(Int_t row, Int_t col, Double_t val)
voidAliMatrixSq::Swap(int& r, int& c) const

Data Members

public:
enum TMatrixTBase::[unnamed] { kSizeMax
kWorkMax
};
enum TMatrixTBase::EMatrixStatusBits { kStatus
};
enum TObject::EStatusBits { kCanDelete
kMustCleanup
kObjInCanvas
kIsReferenced
kHasUUID
kCannotPick
kNoContextMenu
kInvalidObject
};
enum TObject::[unnamed] { kIsOnHeap
kNotDeleted
kZombie
kBitMask
kSingleKey
kOverwrite
kWriteDelete
};
protected:
Int_tTMatrixTBase<double>::fColLwblower bound of the col index
Double_t*fElemsElements booked by constructor
Double_t**fElemsAddElements (rows) added dynamicaly
Bool_tTMatrixTBase<double>::fIsOwner!default kTRUE, when Use array kFALSE
Int_tTMatrixTBase<double>::fNcolsnumber of columns
Int_tTMatrixTBase<double>::fNelemsnumber of elements in matrix
Int_tTMatrixTBase<double>::fNrowIndexlength of row index array (= fNrows+1) wich is only used for sparse matrices
Int_tTMatrixTBase<double>::fNrowsnumber of rows
Int_tTMatrixTBase<double>::fRowLwblower bound of the row index
Bool_tAliMatrixSq::fSymmetricis the matrix symmetric? Only lower triangle is filled
doubleTMatrixTBase<double>::fTolsqrt(epsilon); epsilon is smallest number number so that 1+epsilon > 1
static AliSymMatrix*fgBufferbuffer for fast solution
static Int_tfgCopyCntmatrix copy counter

Class Charts

Inheritance Chart:
TMatrixTBase
AliMatrixSq
AliSymMatrix

Function documentation

AliSymMatrix()
 default constructor
AliSymMatrix(Int_t size)
constructor for matrix with defined size
AliSymMatrix(const AliSymMatrix& mat)
 copy constructor
~AliSymMatrix()
AliSymMatrix& operator=(const AliSymMatrix& src)
 assignment operator
AliSymMatrix& operator+=(const AliSymMatrix& src)
 add operator
void Clear(Option_t* option = "")
 clear dynamic part
Float_t GetDensity() const
 get fraction of non-zero elements
void Print(Option_t* option = "") const
 print itself
void MultiplyByVec(const Double_t* vecIn, Double_t* vecOut) const
 fill vecOut by matrix*vecIn
 vector should be of the same size as the matrix
AliSymMatrix* DecomposeChol()
 Return a matrix with Choleski decomposition
 Adopted from Numerical Recipes in C, ch.2-9, http://www.nr.com
 consturcts Cholesky decomposition of SYMMETRIC and
 POSITIVELY-DEFINED matrix a (a=L*Lt)
 Only upper triangle of the matrix has to be filled.
 In opposite to function from the book, the matrix is modified:
 lower triangle and diagonal are refilled.

Bool_t InvertChol()
 Invert matrix using Choleski decomposition

void InvertChol(AliSymMatrix* mchol)
 Invert matrix using Choleski decomposition, provided the Cholseki's L matrix

Bool_t SolveChol(Double_t* brhs, Bool_t invert = kFALSE)
 Adopted from Numerical Recipes in C, ch.2-9, http://www.nr.com
 Solves the set of n linear equations A x = b,
 where a is a positive-definite symmetric matrix.
 a[1..n][1..n] is the output of the routine CholDecomposw.
 Only the lower triangle of a is accessed. b[1..n] is input as the
 right-hand side vector. The solution vector is returned in b[1..n].

Bool_t SolveChol(TVectorD& brhs, Bool_t invert = kFALSE)
Bool_t SolveChol(Double_t* brhs, Double_t* bsol, Bool_t invert = kFALSE)
Bool_t SolveChol(const TVectorD& brhs, TVectorD& bsol, Bool_t invert = kFALSE)
void AddRows(int nrows = 1)
 add empty rows
void Reset()
 if additional rows exist, regularize it
Double_t* GetRow(Int_t r)
 get pointer on the row
int SolveSpmInv(double* vecB, Bool_t stabilize = kTRUE)
   Solution a la MP1: gaussian eliminations
  Obtain solution of a system of linear equations with symmetric matrix
  and the inverse (using 'singular-value friendly' GAUSS pivot)

Int_t GetIndex(Int_t row, Int_t col) const
Double_t operator()(Int_t row, Int_t col) const
Double_t& operator()(Int_t row, Int_t col)
void MultiplyByVec(const TVectorD& vecIn, TVectorD& vecOut) const
void Scale(Double_t coeff)
void AddToRow(Int_t r, Double_t* valc, Int_t* indc, Int_t n)
Int_t GetSize() const
{return fNrowIndex;}
Int_t GetSizeUsed() const
{return fRowLwb;}
Int_t GetSizeBooked() const
{return fNcols;}
Int_t GetSizeAdded() const
{return fNrows;}
Double_t DiagElem(Int_t r) const
{return (*(const AliSymMatrix*)this)(r,r);}
Double_t& DiagElem(Int_t r)
{return (*this)(r,r);}
void SetSizeUsed(Int_t sz)
{fRowLwb = sz;}
const Double_t * GetMatrixArray() const
 ---------------------------------- Dummy methods of MatrixBase
{return fElems;}
Double_t * GetMatrixArray()
{return (Double_t*)fElems;}
const Int_t * GetRowIndexArray() const
{Error("GetRowIndexArray","Dummy"); return 0;}
Int_t * GetRowIndexArray()
{Error("GetRowIndexArray","Dummy"); return 0;}
const Int_t * GetColIndexArray() const
{Error("GetColIndexArray","Dummy"); return 0;}
Int_t * GetColIndexArray()
{Error("GetColIndexArray","Dummy"); return 0;}
TMatrixDBase & SetRowIndexArray(Int_t* )
{Error("SetRowIndexArray","Dummy"); return *this;}
TMatrixDBase & SetColIndexArray(Int_t* )
{Error("SetColIndexArray","Dummy"); return *this;}
TMatrixDBase & GetSub(Int_t , Int_t , Int_t , Int_t , TMatrixDBase& , Option_t* ) const
{Error("GetSub","Dummy"); return *((TMatrixDBase*)this);}
TMatrixDBase & SetSub(Int_t , Int_t , const TMatrixDBase& )
{Error("GetSub","Dummy"); return *this;}
TMatrixDBase & ResizeTo(Int_t , Int_t , Int_t )
{Error("ResizeTo","Dummy"); return *this;}
TMatrixDBase & ResizeTo(Int_t , Int_t , Int_t , Int_t , Int_t )
 ----------------------------- Choleski methods ----------------------------------------
{Error("ResizeTo","Dummy"); return *this;}
Double_t GetEl(Int_t row, Int_t col) const
{return operator()(row,col);}
void SetEl(Int_t row, Int_t col, Double_t val)
{operator()(row,col) = val;}