void Solve(const TVectorD &rhs,TVectorD &sol) {Solve(rhs.GetMatrixArray(),sol.GetMatrixArray());}
//
void Print(Option_t* option="") const;
- void SetDecomposed() {SetBit(kDecomposedBit);}
+ void SetDecomposed(Bool_t v=kTRUE) {SetBit(kDecomposedBit,v);}
Bool_t IsDecomposed() const {return TestBit(kDecomposedBit);}
//
- void MultiplyByVec(Double_t* vecIn, Double_t* vecOut) const;
- void MultiplyByVec(TVectorD &vecIn, TVectorD &vecOut) const;
+ void MultiplyByVec(const Double_t* vecIn, Double_t* vecOut) const;
+ void MultiplyByVec(const TVectorD &vecIn, TVectorD &vecOut) const;
void AddToRow(Int_t r, Double_t *valc,Int_t *indc,Int_t n);
//
// protected:
}
//___________________________________________________________
-inline void AliSymBDMatrix::MultiplyByVec(TVectorD &vecIn, TVectorD &vecOut) const
+inline void AliSymBDMatrix::MultiplyByVec(const TVectorD &vecIn, TVectorD &vecOut) const
{
MultiplyByVec(vecIn.GetMatrixArray(), vecOut.GetMatrixArray());
}