|
ROOT
Reference Guide |
|
Go to the documentation of this file.
12 #ifndef ROOT_TMatrixTUtils
13 #define ROOT_TMatrixTUtils
38 #include <initializer_list>
40 template<
class Element>
class TVectorT;
41 template<
class Element>
class TMatrixT;
140 if (acoln < fMatrix->GetNcols() && acoln >= 0)
143 Error(
"operator()",
"Request col(%d) outside matrix range of %d - %d",
161 inline Element *
GetPtr()
const {
return const_cast<Element *
>(this->fPtr); }
166 const Int_t acoln = i-this->fMatrix->GetColLwb();
167 if (acoln < this->fMatrix->GetNcols() || acoln >= 0)
168 return (this->fPtr)[acoln];
170 Error(
"operator()",
"Request col(%d) outside matrix range of %d - %d",
171 i,this->fMatrix->GetColLwb(),this->fMatrix->GetColLwb()+this->fMatrix->GetNcols());
178 const Int_t acoln = i-this->fMatrix->GetColLwb();
179 if (acoln < this->fMatrix->GetNcols() && acoln >= 0)
180 return (
const_cast<Element *
>(this->fPtr))[acoln];
182 Error(
"operator()",
"Request col(%d) outside matrix range of %d - %d",
183 i,this->fMatrix->GetColLwb(),this->fMatrix->GetColLwb()+this->fMatrix->GetNcols());
191 void Assign (Element val);
192 void operator= (std::initializer_list<Element>
l);
240 if (arown < fMatrix->GetNrows() && arown >= 0)
243 Error(
"operator()",
"Request row(%d) outside matrix range of %d - %d",
261 inline Element *
GetPtr()
const {
return const_cast<Element *
>(this->fPtr); }
266 const Int_t arown = i-this->fMatrix->GetRowLwb();
267 if (arown < this->fMatrix->GetNrows() && arown >= 0)
268 return (this->fPtr)[arown*this->fInc];
270 Error(
"operator()",
"Request row(%d) outside matrix range of %d - %d",
271 i,this->fMatrix->GetRowLwb(),this->fMatrix->GetRowLwb()+this->fMatrix->GetNrows());
278 const Int_t arown = i-this->fMatrix->GetRowLwb();
280 if (arown < this->fMatrix->GetNrows() && arown >= 0)
281 return (
const_cast<Element *
>(this->fPtr))[arown*this->fInc];
283 Error(
"operator()",
"Request row(%d) outside matrix range of %d - %d",
284 i,this->fMatrix->GetRowLwb(),this->fMatrix->GetRowLwb()+this->fMatrix->GetNrows());
291 void Assign (Element val);
294 void operator= (std::initializer_list<Element>
l);
339 if (i < fNdiag && i >= 0)
342 Error(
"operator()",
"Request diagonal(%d) outside matrix range of 0 - %d",i,
fNdiag);
361 inline Element *
GetPtr()
const {
return const_cast<Element *
>(this->fPtr); }
365 if (i < this->fNdiag && i >= 0)
366 return (this->fPtr)[i*this->fInc];
368 Error(
"operator()",
"Request diagonal(%d) outside matrix range of 0 - %d",i,this->fNdiag);
374 if (i < this->fNdiag && i >= 0)
375 return (
const_cast<Element *
>(this->fPtr))[i*this->fInc];
377 Error(
"operator()",
"Request diagonal(%d) outside matrix range of 0 - %d",i,this->fNdiag);
378 return (
const_cast<Element *
>(this->fPtr))[0];
427 if (i < fNelems && i >= 0)
430 Error(
"operator()",
"Request element(%d) outside matrix range of 0 - %d",i,
fNelems);
447 inline Element *
GetPtr()
const {
return const_cast<Element *
>(this->fPtr); }
451 if (i < this->fNelems && i >= 0)
452 return (this->fPtr)[i];
454 Error(
"operator()",
"Request element(%d) outside matrix range of 0 - %d",i,this->fNelems);
460 if (i < this->fNelems && i >= 0)
461 return (
const_cast<Element *
>(this->fPtr))[i];
463 Error(
"operator()",
"Request element(%d) outside matrix range of 0 - %d",i,this->fNelems);
515 const Element *ptr =
fMatrix->GetMatrixArray();
517 Error(
"operator()",
"Request row(%d) outside matrix range of 0 - %d",rown,
fNrowsSub);
521 Error(
"operator()",
"Request column(%d) outside matrix range of 0 - %d",coln,
fNcolsSub);
545 const Element *ptr = this->fMatrix->GetMatrixArray();
546 if (rown >= this->fNrowsSub || rown < 0) {
547 Error(
"operator()",
"Request row(%d) outside matrix range of 0 - %d",rown,this->fNrowsSub);
550 if (coln >= this->fNcolsSub || coln < 0) {
551 Error(
"operator()",
"Request column(%d) outside matrix range of 0 - %d",coln,this->fNcolsSub);
554 const Int_t index = (rown+this->fRowOff)*this->fMatrix->GetNcols()+coln+this->fColOff;
555 return (
const_cast<Element *
>(ptr))[index];
622 inline Element *
GetDataPtr()
const {
return const_cast<Element *
>(this->fDataPtr); }
684 inline Element *
GetDataPtr()
const {
return const_cast<Element *
>(this->fDataPtr); }
TMatrixTColumn_const< Element > & operator=(const TMatrixTColumn_const< Element > &trc)
const Element & operator()(Int_t i) const
virtual ~TElementActionT()
TMatrixTRow_const(const TMatrixTRow_const< Element > &trc)
const Element * GetPtr() const
TMatrixTSparseRow_const(const TMatrixTSparseRow_const< Element > &trc)
TMatrixTRow_const< Element > & operator=(const TMatrixTRow_const< Element > &trc)
static Element & NaNValue()
void Rank1Update(const TVectorT< Element > &vec, Element alpha=1.0)
Perform a rank 1 operation on the matrix: A += alpha * v * v^T.
void Assign(Element val)
Assign val to every element of the matrix column.
Element * GetDataPtr() const
TMatrixTDiag< Element > & operator=(const TMatrixTDiag< Element > &d)
Element & operator()(Int_t i)
Int_t GetRowIndex() const
Element operator()(Int_t i) const
void Assign(Element val)
Assign val to every element of the matrix row.
void operator+=(Element val)
Add val to every element of the matrix.
virtual ~TMatrixTColumn_const()
void operator*=(Element val)
Multiply every element of the sub matrix by val .
const Element & operator()(Int_t i) const
Element & operator()(Int_t i)
const Element & operator[](Int_t i) const
virtual ~TMatrixTSparseDiag_const()
TMatrixTFlat< Element > & operator=(const TMatrixTFlat< Element > &f)
virtual ~TElementPosActionT()
TMatrixTSparseRow< Element > & operator=(const TMatrixTSparseRow< Element > &r)
const TMatrixTBase< Element > * GetMatrix() const
TMatrixTColumn_const(const TMatrixTColumn_const< Element > &trc)
virtual ~TMatrixTFlat_const()
void operator*=(Element val)
Multiply every element of the matrix diagonal by val.
Element operator[](Int_t i) const
Element & operator[](Int_t i)
TElementPosActionT< Element > & operator=(const TElementPosActionT< Element > &)
static constexpr double s
Element operator[](Int_t i) const
const Element * GetPtr() const
const TMatrixTSparse< Element > * fMatrix
const Element & operator()(Int_t i) const
const Element & operator()(Int_t i) const
void operator*=(Element val)
Multiply every element of the matrix row with val.
Element & operator[](Int_t i)
Element operator()(Int_t i) const
void operator*=(Element val)
Multiply every element of the matrix with val.
void operator=(std::initializer_list< Element > l)
Element & operator()(Int_t rown, Int_t coln)
void operator+=(Element val)
Add val to every non-zero (!) element of the matrix row.
TMatrixTSparseDiag_const()
const TMatrixTBase< Element > * GetMatrix() const
void operator=(Element val)
Assign val to every element of the matrix diagonal.
const Int_t * GetColPtr() const
Element operator()(Int_t i) const
TMatrixTSparseDiag_const< Element > & operator=(const TMatrixTSparseDiag_const< Element > &trc)
const TMatrixTBase< Element > * fMatrix
const TMatrixTBase< Element > * GetMatrix() const
TMatrixTSub< Element > & operator=(const TMatrixTSub< Element > &s)
Element * GetDataPtr() const
Element operator[](Int_t i) const
const TMatrixTBase< Element > * GetMatrix() const
Element & operator[](Int_t i)
const Element & operator()(Int_t rown, Int_t coln) const
const Element * GetPtr() const
virtual void Operation(Element &element) const =0
Element & operator[](Int_t i)
TMatrixTSparseDiag_const(const TMatrixTSparseDiag_const< Element > &trc)
const TMatrixTBase< Element > * fMatrix
void operator=(Element val)
Assign val to every element of the sub matrix.
const Element & operator()(Int_t i) const
Element operator()(Int_t i) const
Element & operator[](Int_t i)
void operator*=(Element val)
Assign val to every element of the matrix diagonal.
const Element & operator[](Int_t i) const
void operator+=(Element val)
Add val to every element of the sub matrix.
virtual ~TMatrixTSparseRow_const()
const Element & operator[](Int_t i) const
virtual void Operation(Element &element) const =0
const TMatrixTBase< Element > * fMatrix
const Element & operator[](Int_t i) const
TMatrixTSparseRow_const()
const Element * GetDataPtr() const
TMatrixTColumn< Element > & operator=(const TMatrixTColumn< Element > &c)
virtual ~TMatrixTSub_const()
Element & operator()(Int_t i)
const Element & operator[](Int_t i) const
TMatrixTDiag_const(const TMatrixTDiag_const< Element > &trc)
Element & operator[](Int_t i)
Element operator[](Int_t i) const
TMatrixTFlat_const< Element > & operator=(const TMatrixTFlat_const< Element > &trc)
virtual ~TMatrixTDiag_const()
void operator=(Element val)
Assign val to every non-zero (!) element of the matrix row.
TMatrixTRow< Element > & operator=(const TMatrixTRow< Element > &r)
#define ClassDef(name, id)
const TMatrixTSparse< Element > * fMatrix
void operator*=(Element val)
Multiply every element of the matrix column with val.
void operator+=(Element val)
Add val to every element of the matrix diagonal.
Int_t GetColIndex() const
void operator=(Element val)
TMatrixTSparseRow_const< Element > & operator=(const TMatrixTSparseRow_const< Element > &trc)
const Element * GetDataPtr() const
const Element & operator()(Int_t i) const
const Element & operator[](Int_t i) const
const TMatrixTBase< Element > * GetMatrix() const
static constexpr double ms
TMatrixTSparseDiag< Element > & operator=(const TMatrixTSparseDiag< Element > &d)
Int_t GetRowIndex() const
virtual ~TMatrixTRow_const()
const TMatrixTBase< Element > * GetMatrix() const
void operator=(Element val)
Assign val to every element of the matrix diagonal.
void operator+=(Element val)
Assign val to every element of the matrix diagonal.
void operator+=(Element val)
Add val to every element of the matrix column.
TMatrixTFlat_const(const TMatrixTFlat_const< Element > &trc)
TElementActionT & operator=(const TElementActionT< Element > &)
const Element & operator[](Int_t i) const
void operator=(Element val)
Assign val to every element of the matrix.
const TMatrixTBase< Element > * GetMatrix() const
void operator*=(Element val)
Multiply every element of the matrix row by val.
void operator+=(Element val)
Add val to every element of the matrix row.
const TMatrixTBase< Element > * fMatrix
TMatrixTDiag_const< Element > & operator=(const TMatrixTDiag_const< Element > &trc)
const Element * GetPtr() const
const Element & operator()(Int_t i) const
const TMatrixTBase< Element > * fMatrix
Element & operator()(Int_t i)
const Element & operator[](Int_t i) const
const Element & operator()(Int_t i) const
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
Double_t Drand(Double_t &ix)
Random number generator [0....1] with seed ix.