12 #ifndef ROOT_TDecompBK
13 #define ROOT_TDecompBK
21 #ifndef ROOT_TDecompBase
24 #ifndef ROOT_TMatrixDSym
65 {
MayNotUse(
"Det(Double_t&,Double_t&)"); }
virtual Bool_t TransSolve(TMatrixDColumn &b)
virtual Bool_t Solve(TVectorD &b)
Solve Ax=b assuming the BK form of A is stored in fU . Solution returned in b.
void MayNotUse(const char *method) const
Use this method to signal that a method (defined in a base class) may not be called in a derived clas...
virtual const TMatrixDBase & GetDecompMatrix() const
double inv(double x)
For comparisons.
#define ClassDef(name, id)
void Print(Option_t *opt="") const
Print the class members.
virtual TVectorD TransSolve(const TVectorD &b, Bool_t &ok)
TDecompBK & operator=(const TDecompBK &source)
Assigment operator.
virtual Int_t GetNrows() const
Bool_t TestBit(UInt_t f) const
virtual void Det(Double_t &, Double_t &)
Matrix determinant det = d1*TMath::Power(2.,d2)
virtual void SetMatrix(const TMatrixDSym &a)
Set the matrix to be decomposed, decomposition status is reset.
virtual Int_t GetNcols() const
virtual Bool_t TransSolve(TVectorD &b)
virtual TVectorD Solve(const TVectorD &b, Bool_t &ok)
virtual Bool_t Decompose()
Matrix A is decomposed in components U and D so that A = U*D*U^T If the decomposition succeeds...