removing obsolete component from libAliHLTTrigger

[u/mrichter/AliRoot.git] / STEER / AliMatrixSparse.h

diff --git a/STEER/AliMatrixSparse.h b/STEER/AliMatrixSparse.h
index 9dcf5195925a3402ea1480ab15543193ce85ff66..12898b4a4e4951fa2682627356c390a6501d2678 100644 (file)
--- a/STEER/AliMatrixSparse.h
+++ b/STEER/AliMatrixSparse.h
@@ -24,7 +24,11 @@ class AliMatrixSparse : public AliMatrixSq
   //
   AliVectorSparse* GetRow(Int_t ir)        const {return (ir<fNcols) ? fVecs[ir] : 0;}
   AliVectorSparse* GetRowAdd(Int_t ir);
-
+  //
+  virtual Int_t   GetSize()                            const {return fNrows;}
+  virtual Int_t   GetNRows()                           const {return fNrows;}
+  virtual Int_t   GetNCols()                           const {return fNcols;}
+  //
   void Clear(Option_t* option="");
   void Reset()                            {for (int i=fNcols;i--;) GetRow(i)->Reset();}
   void Print(Option_t* option="")                      const;
@@ -38,8 +42,8 @@ class AliMatrixSparse : public AliMatrixSq
   Double_t&        DiagElem(Int_t r);
   void             SortIndices(Bool_t valuesToo=kFALSE);
   //
-  void MultiplyByVec(TVectorD &vecIn, TVectorD &vecOut) const; 
-  void MultiplyByVec(Double_t* vecIn, Double_t* vecOut) const;
+  void MultiplyByVec(const TVectorD &vecIn, TVectorD &vecOut) const; 
+  void MultiplyByVec(const Double_t* vecIn, Double_t* vecOut) const;
   //
   void AddToRow(Int_t r, Double_t *valc,Int_t *indc,Int_t n);
   //
@@ -51,7 +55,7 @@ class AliMatrixSparse : public AliMatrixSq
 };
 
 //___________________________________________________
-inline void AliMatrixSparse::MultiplyByVec(TVectorD &vecIn, TVectorD &vecOut) const 
+inline void AliMatrixSparse::MultiplyByVec(const TVectorD &vecIn, TVectorD &vecOut) const 
 {
   MultiplyByVec((Double_t*)vecIn.GetMatrixArray(),(Double_t*)vecOut.GetMatrixArray());
 }
@@ -81,6 +85,7 @@ inline Double_t& AliMatrixSparse::operator()(Int_t row,Int_t col)
   //  printf("M: findindexAdd\n");
   if (IsSymmetric() && col>row) Swap(row,col); 
   AliVectorSparse* rowv = GetRowAdd(row);
+  if (col>=fNcols) fNcols = col+1;
   return rowv->FindIndexAdd(col);
 }
 
@@ -98,8 +103,11 @@ inline Double_t AliMatrixSparse::DiagElem(Int_t row) const
 inline Double_t &AliMatrixSparse::DiagElem(Int_t row)
 {
   AliVectorSparse* rowv = GetRowAdd(row);
-  if (IsSymmetric()) return (rowv->GetNElems()>0 && rowv->GetLastIndex()==row) ? 
+  if (row>=fNcols) fNcols = row+1;
+  if (IsSymmetric()) {
+    return (rowv->GetNElems()>0 && rowv->GetLastIndex()==row) ? 
                       rowv->GetLastElem() : rowv->FindIndexAdd(row);
+  }
   else return rowv->FindIndexAdd(row);
   //
 }