58 TVectorT(
const TVectorT <Element> &another);
59 TVectorT(
const TMatrixTRow_const <Element> &mr);
61 TVectorT(
const TMatrixTDiag_const <Element> &md);
111 Element
Norm1 ()
const;
115 Element
Sum ()
const;
116 Element
Min ()
const;
117 Element
Max ()
const;
133 Error(
"operator=(const TVectorT2 &)",
"vectors not compatible");
199 return Use(
v.GetLwb(),
v.GetUpb(),
v.GetMatrixArray());
204 return Use(
v.GetLwb(),
v.GetUpb(),
v.GetMatrixArray());
209 this->GetSub(row_lwb,row_upb,tmp,option);
218 const Int_t aind = ind-fRowLwb;
219 if (aind >= fNrows || aind < 0) {
220 Error(
"operator()",
"Request index(%d) outside vector range of %d - %d",ind,fRowLwb,fRowLwb+fNrows);
224 return fElements[aind];
231 const Int_t aind = ind-fRowLwb;
232 if (aind >= fNrows || aind < 0) {
233 Error(
"operator()",
"Request index(%d) outside vector range of %d - %d",ind,fRowLwb,fRowLwb+fNrows);
237 return fElements[aind];
241template<
class Element>
Bool_t operator== (
const TVectorT <Element> &source1,
const TVectorT <Element> &source2);
244template<
class Element> Element
operator* (
const TVectorT <Element> &source1,
const TVectorT <Element> &source2);
249template<
class Element>
253template<
class Element> Element
Dot (
const TVectorT <Element> &source1,
const TVectorT <Element> &source2);
254template <
class Element1,
class Element2>
256template <
class Element1,
class Element2,
class Element3>
258template <
class Element1,
class Element2,
class Element3>
259 Element1
Mult (
const TVectorT <Element1> &
v1,
const TMatrixT <Element2> &
m,
const TVectorT <Element3> &
v2);
261template<
class Element>
TVectorT<Element> &
Add ( TVectorT <Element> &target, Element scalar,
const TVectorT <Element> &source);
262template<
class Element>
TVectorT<Element> &
Add ( TVectorT <Element> &target, Element scalar,
const TMatrixT <Element> &
a,
264template<
class Element>
TVectorT<Element> &
Add ( TVectorT <Element> &target, Element scalar,
const TMatrixTSym <Element> &
a,
269 const TVectorT <Element> &source2);
271 const TVectorT <Element> &source2,
const TVectorT <Element> &select);
272template<
class Element>
TVectorT<Element> &
AddElemDiv ( TVectorT <Element> &target, Element scalar,
const TVectorT <Element> &source1,
273 const TVectorT <Element> &source2);
274template<
class Element>
TVectorT<Element> &
AddElemDiv ( TVectorT <Element> &target, Element scalar,
const TVectorT <Element> &source1,
275 const TVectorT <Element> &source2,
const TVectorT <Element> &select);
277template<
class Element>
TVectorT<Element> &
ElementMult ( TVectorT <Element> &target,
const TVectorT <Element> &source,
const TVectorT <Element> &select);
279template<
class Element>
TVectorT<Element> &
ElementDiv ( TVectorT <Element> &target,
const TVectorT <Element> &source,
const TVectorT <Element> &select);
286template<
class Element>
void Compare (
const TVectorT <Element> &source1,
const TVectorT <Element> &source2);
size_t size(const MatrixT &matrix)
retrieve the size of a square matrix
#define ClassDef(name, id)
void Error(const char *location, const char *msgfmt,...)
Use this function in case an error occurred.
TClass instances represent classes, structs and namespaces in the ROOT type system.
static Element & NaNValue()
Mother of all ROOT objects.
TObject & operator=(const TObject &rhs)
TObject assignment operator.
R__ALWAYS_INLINE Bool_t TestBit(UInt_t f) const
void SetBit(UInt_t f, Bool_t set)
Set or unset the user status bits as specified in f.
virtual void Error(const char *method, const char *msgfmt,...) const
Issue error message.
TVectorT< Element > & Zero()
Set vector elements to zero.
Bool_t operator>=(Element val) const
Are all vector elements >= val?
Element * New_m(Int_t size)
default kTRUE, when Use array kFALSE
TVectorT< Element > & Apply(const TElementActionT< Element > &action)
Apply action to each element of the vector.
void SetElements(const Element *elements)
void Allocate(Int_t nrows, Int_t row_lwb=0, Int_t init=0)
Allocate new vector.
TVectorT< Element > & Sqr()
Square each element of the vector.
Element & operator[](Int_t index)
Element Min() const
return minimum vector element value
Int_t Memcpy_m(Element *newp, const Element *oldp, Int_t copySize, Int_t newSize, Int_t oldSize)
Copy copySize doubles from *oldp to *newp .
void Clear(Option_t *="")
Element & operator()(Int_t index)
TVectorT< Element > & Use(TVectorT< Element > &v)
Element Norm1() const
Compute the 1-norm of the vector SUM{ |v[i]| }.
void Delete_m(Int_t size, Element *&)
Delete data pointer m, if it was assigned on the heap.
TVectorT(const TVectorT< Element2 > &another)
Element NormInf() const
Compute the infinity-norm of the vector MAX{ |v[i]| }.
Element Sum() const
Compute sum of elements.
void AddSomeConstant(Element val, const TVectorT< Element > &select)
Add to vector elements as selected through array select the value val.
Bool_t operator!=(Element val) const
Are all vector elements not equal to val?
Bool_t MatchesNonZeroPattern(const TVectorT< Element > &select)
Check if vector elements as selected through array select are non-zero.
TVectorT< Element > & GetSub(Int_t row_lwb, Int_t row_upb, TVectorT< Element > &target, Option_t *option="S") const
Get subvector [row_lwb..row_upb]; The indexing range of the returned vector depends on the argument o...
TVectorT< Element > & Shift(Int_t row_shift)
TVectorT< Element > & Abs()
Take an absolute value of a vector, i.e. apply Abs() to each element.
Bool_t operator==(Element val) const
Are all vector elements equal to val?
Bool_t operator>(Element val) const
Are all vector elements > val?
TVectorT< Element > & ResizeTo(Int_t n)
TVectorT< Element > & ResizeTo(Int_t lwb, Int_t upb)
Resize the vector to [lwb:upb] .
TVectorT< Element > & operator=(const TVectorT< Element > &source)
Notice that this assignment does NOT change the ownership : if the storage space was adopted,...
TVectorT< Element > & SetSub(Int_t row_lwb, const TVectorT< Element > &source)
Insert vector source starting at [row_lwb], thereby overwriting the part [row_lwb....
void Randomize(Element alpha, Element beta, Double_t &seed)
randomize vector elements value
TVectorT< Element > & Use(Int_t n, Element *data)
void Add(const TVectorT< Element > &v)
Add vector v to this vector.
const Element & operator()(Int_t index) const
TVectorT< Element > & operator+=(Element val)
Add val to every element of the vector.
TVectorT< Element > & Use(Int_t lwb, Int_t upb, Element *data)
Use the array data to fill the vector lwb..upb].
const TVectorT< Element > & Use(Int_t lwb, Int_t upb, const Element *data) const
TVectorT< Element > & SelectNonZeros(const TVectorT< Element > &select)
Keep only element as selected through array select non-zero.
Element Max() const
return maximum vector element value
Bool_t operator<(Element val) const
Are all vector elements < val?
Bool_t operator<=(Element val) const
Are all vector elements <= val?
TVectorT< Element > & operator*=(Element val)
Multiply every element of the vector with val.
void Print(Option_t *option="") const
Print the vector as a list of elements.
const TVectorT< Element > & Use(const TVectorT< Element > &v) const
Bool_t SomePositive(const TVectorT< Element > &select)
Check if vector elements as selected through array select are all positive.
const Element & operator[](Int_t index) const
Int_t NonZeros() const
Compute the number of elements != 0.0.
Bool_t fIsOwner
data container
TVectorT< Element > & ResizeTo(const TVectorT< Element > &v)
const TVectorT< Element > & Use(Int_t n, const Element *data) const
TVectorT< Element > & Invert()
v[i] = 1/v[i]
Int_t GetNoElements() const
Element Norm2Sqr() const
Compute the square of the 2-norm SUM{ v[i]^2 }.
Element * GetMatrixArray()
const Element * GetMatrixArray() const
TVectorT< Element > & operator-=(Element val)
Subtract val from every element of the vector.
Element fDataStack[kSizeMax]
TVectorT< Element > GetSub(Int_t row_lwb, Int_t row_upb, Option_t *option="S") const
TVectorT< Element > & Sqrt()
Take square root of all elements.
Bool_t VerifyVectorValue(const TVectorT< Element > &m, Element val, Int_t verbose, Element maxDevAllow)
Validate that all elements of vector have value val within maxDevAllow .
TVectorT< Element > & AddElemDiv(TVectorT< Element > &target, Element scalar, const TVectorT< Element > &source1, const TVectorT< Element > &source2)
Modify addition: target += scalar * ElementDiv(source1,source2) .
TMatrixT< Element > operator+(const TMatrixT< Element > &source1, const TMatrixT< Element > &source2)
operation this = source1+source2
void Compare(const TMatrixTBase< Element > &m1, const TMatrixTBase< Element > &m2)
Compare two matrices and print out the result of the comparison.
TMatrixT< Element > & ElementMult(TMatrixT< Element > &target, const TMatrixT< Element > &source)
Multiply target by the source, element-by-element.
TVectorT< Element > & AddElemMult(TVectorT< Element > &target, Element scalar, const TVectorT< Element > &source1, const TVectorT< Element > &source2)
Modify addition: target += scalar * ElementMult(source1,source2) .
TMatrixT< Element > & Add(TMatrixT< Element > &target, Element scalar, const TMatrixT< Element > &source)
Modify addition: target += scalar * source.
TMatrixT< Element1 > OuterProduct(const TVectorT< Element1 > &v1, const TVectorT< Element2 > &v2)
Return the matrix M = v1 * v2'.
Bool_t VerifyVectorIdentity(const TVectorT< Element > &m1, const TVectorT< Element > &m2, Int_t verbose, Element maxDevAllow)
Verify that elements of the two vectors are equal within maxDevAllow .
TMatrixT< Element > & ElementDiv(TMatrixT< Element > &target, const TMatrixT< Element > &source)
Divide target by the source, element-by-element.
Element1 Mult(const TVectorT< Element1 > &v1, const TMatrixT< Element2 > &m, const TVectorT< Element3 > &v2)
Perform v1 * M * v2, a scalar result.
TMatrixT< Element > operator-(const TMatrixT< Element > &source1, const TMatrixT< Element > &source2)
operation this = source1-source2
Bool_t AreCompatible(const TVectorT< Element1 > &v1, const TVectorT< Element2 > &v2, Int_t verbose=0)
Check if v1 and v2 are both valid and have the same shape.
Bool_t operator==(const TMatrixTBase< Element > &m1, const TMatrixTBase< Element > &m2)
Check to see if two matrices are identical.
TMatrixT< Element > operator*(Element val, const TMatrixT< Element > &source)
operation this = val*source
Bool_t AreCompatible(const TMatrixTBase< Element1 > &m1, const TMatrixTBase< Element2 > &m2, Int_t verbose=0)
Check that matrice sm1 and m2 areboth valid and have identical shapes .