#ifndef ROOT_TMatrixTBase
#define ROOT_TMatrixTBase
#include <limits>
#ifndef ROOT_TROOT
#include "TROOT.h"
#endif
#ifndef ROOT_TClass
#include "TClass.h"
#endif
#ifndef ROOT_TMatrixFBasefwd
#include "TMatrixFBasefwd.h"
#endif
#ifndef ROOT_TMatrixDBasefwd
#include "TMatrixDBasefwd.h"
#endif
#ifndef ROOT_TVectorFfwd
#include "TVectorFfwd.h"
#endif
#ifndef ROOT_TVectorDfwd
#include "TVectorDfwd.h"
#endif
#ifndef ROOT_TError
#include "TError.h"
#endif
#ifndef ROOT_TMath
#include "TMath.h"
#endif
template<class Element> class TVectorT;
template<class Element> class TElementActionT;
template<class Element> class TElementPosActionT;
R__EXTERN Int_t gMatrixCheck;
template<class Element> class TMatrixTBase : public TObject {
private:
Element *GetElements();
protected:
Int_t fNrows;
Int_t fNcols;
Int_t fRowLwb;
Int_t fColLwb;
Int_t fNelems;
Int_t fNrowIndex;
Element fTol;
Bool_t fIsOwner;
static void DoubleLexSort (Int_t n,Int_t *first,Int_t *second,Element *data);
static void IndexedLexSort(Int_t n,Int_t *first,Int_t swapFirst,
Int_t *second,Int_t swapSecond,Int_t *index);
enum {kSizeMax = 25};
enum {kWorkMax = 100};
enum EMatrixStatusBits {
kStatus = BIT(14)
};
public:
TMatrixTBase():
fNrows(0), fNcols(0), fRowLwb(0), fColLwb(0), fNelems(0), fNrowIndex(0),
fTol(0), fIsOwner(kTRUE) { }
virtual ~TMatrixTBase() {}
inline Int_t GetRowLwb () const { return fRowLwb; }
inline Int_t GetRowUpb () const { return fNrows+fRowLwb-1; }
inline Int_t GetNrows () const { return fNrows; }
inline Int_t GetColLwb () const { return fColLwb; }
inline Int_t GetColUpb () const { return fNcols+fColLwb-1; }
inline Int_t GetNcols () const { return fNcols; }
inline Int_t GetNoElements () const { return fNelems; }
inline Element GetTol () const { return fTol; }
virtual const Element *GetMatrixArray () const = 0;
virtual Element *GetMatrixArray () = 0;
virtual const Int_t *GetRowIndexArray() const = 0;
virtual Int_t *GetRowIndexArray() = 0;
virtual const Int_t *GetColIndexArray() const = 0;
virtual Int_t *GetColIndexArray() = 0;
virtual TMatrixTBase<Element> &SetRowIndexArray(Int_t *data) = 0;
virtual TMatrixTBase<Element> &SetColIndexArray(Int_t *data) = 0;
virtual TMatrixTBase<Element> &SetMatrixArray (const Element *data,Option_t *option="");
inline Element SetTol (Element tol);
virtual void Clear (Option_t *option="") = 0;
inline void Invalidate () { SetBit(kStatus); }
inline void MakeValid () { ResetBit(kStatus); }
inline Bool_t IsValid () const { return !TestBit(kStatus); }
inline Bool_t IsOwner () const { return fIsOwner; }
virtual Bool_t IsSymmetric() const;
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 = 0;
virtual TMatrixTBase<Element> &SetSub(Int_t row_lwb,Int_t col_lwb,const TMatrixTBase<Element> &source) = 0;
virtual void GetMatrix2Array(Element *data,Option_t *option="") const;
virtual TMatrixTBase<Element> &InsertRow (Int_t row,Int_t col,const Element *v,Int_t n = -1);
virtual void ExtractRow (Int_t row,Int_t col, Element *v,Int_t n = -1) const;
virtual TMatrixTBase<Element> &Shift (Int_t row_shift,Int_t col_shift);
virtual TMatrixTBase<Element> &ResizeTo (Int_t nrows,Int_t ncols,Int_t nr_nonzeros=-1) = 0;
virtual TMatrixTBase<Element> &ResizeTo (Int_t row_lwb,Int_t row_upb,Int_t col_lwb,Int_t col_upb,Int_t nr_nonzeros=-1) = 0;
virtual Double_t Determinant() const { AbstractMethod("Determinant()"); return 0.; }
virtual void Determinant(Double_t &d1,Double_t &d2) const { AbstractMethod("Determinant()"); d1 = 0.; d2 = 0.; }
virtual TMatrixTBase<Element> &Zero ();
virtual TMatrixTBase<Element> &Abs ();
virtual TMatrixTBase<Element> &Sqr ();
virtual TMatrixTBase<Element> &Sqrt ();
virtual TMatrixTBase<Element> &UnitMatrix ();
virtual TMatrixTBase<Element> &NormByDiag (const TVectorT<Element> &v,Option_t *option="D");
virtual Element RowNorm () const;
virtual Element ColNorm () const;
virtual Element E2Norm () const;
inline Element NormInf () const { return RowNorm(); }
inline Element Norm1 () const { return ColNorm(); }
virtual Int_t NonZeros () const;
virtual Element Sum () const;
virtual Element Min () const;
virtual Element Max () const;
void Draw (Option_t *option="");
void Print(Option_t *name ="") const;
virtual Element operator()(Int_t rown,Int_t coln) const = 0;
virtual Element &operator()(Int_t rown,Int_t coln) = 0;
Bool_t operator==(Element val) const;
Bool_t operator!=(Element val) const;
Bool_t operator< (Element val) const;
Bool_t operator<=(Element val) const;
Bool_t operator> (Element val) const;
Bool_t operator>=(Element val) const;
virtual TMatrixTBase<Element> &Apply(const TElementActionT<Element> &action);
virtual TMatrixTBase<Element> &Apply(const TElementPosActionT<Element> &action);
virtual TMatrixTBase<Element> &Randomize(Element alpha,Element beta,Double_t &seed);
ClassDef(TMatrixTBase,5)
};
template<class Element> Element TMatrixTBase<Element>::SetTol(Element newTol)
{
const Element oldTol = fTol;
if (newTol >= 0.0)
fTol = newTol;
return oldTol;
}
template<class Element> Bool_t operator== (const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2);
template<class Element> Element E2Norm (const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2);
template<class Element1,class Element2>
Bool_t AreCompatible(const TMatrixTBase<Element1> &m1,const TMatrixTBase<Element2> &m2,Int_t verbose=0);
template<class Element> void Compare (const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2);
template<class Element> Bool_t VerifyMatrixValue (const TMatrixTBase<Element> &m,Element val,
Int_t verbose,Element maxDevAllow);
template<class Element> Bool_t VerifyMatrixValue (const TMatrixTBase<Element> &m,Element val,Int_t verbose)
{ return VerifyMatrixValue(m,val,verbose,Element(0.)); }
template<class Element> Bool_t VerifyMatrixValue (const TMatrixTBase<Element> &m,Element val)
{ return VerifyMatrixValue(m,val,1,Element(0.)); }
template<class Element> Bool_t VerifyMatrixIdentity(const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2,
Int_t verbose,Element maxDevAllow);
template<class Element> Bool_t VerifyMatrixIdentity(const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2,Int_t verbose)
{ return VerifyMatrixIdentity(m1,m2,verbose,Element(0.)); }
template<class Element> Bool_t VerifyMatrixIdentity(const TMatrixTBase<Element> &m1,const TMatrixTBase<Element> &m2)
{ return VerifyMatrixIdentity(m1,m2,1,Element(0.)); }
#endif
ROOT page - Class index - Class Hierarchy - Top of the page
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.