//___________________________________________________
inline void AliMatrixSparse::MultiplyByVec(const TVectorD &vecIn, TVectorD &vecOut) const
{
+ // multiplication
MultiplyByVec((Double_t*)vecIn.GetMatrixArray(),(Double_t*)vecOut.GetMatrixArray());
}
//___________________________________________________
inline Double_t AliMatrixSparse::DiagElem(Int_t row) const
{
+ // get diag elem
AliVectorSparse* rowv = GetRow(row);
if (!rowv) return 0;
if (IsSymmetric()) return (rowv->GetNElems()>0 && rowv->GetLastIndex()==row) ? rowv->GetLastElem() : 0.;
//___________________________________________________
inline Double_t &AliMatrixSparse::DiagElem(Int_t row)
{
+ // get diag elem
AliVectorSparse* rowv = GetRowAdd(row);
if (row>=fNcols) fNcols = row+1;
if (IsSymmetric()) {