SVector Class Properties

The template ROOT::Math::SVector class has 2 template parameters which define, at compile time, its properties. These are:

- type of the contained elements, for example
*float*or double. - size of the vector.

The following constructors are available to create a vector:

- Default constructor for a zero vector (all elements equal to zero)
- Constructor (and assignment) from a vector expression, like v = p*q + w. Due to the expression template technique, no temporary objects are created in this operation.
- Construct a vector passing directly the elements. This is possible only for vector up to size 10.
- Constructor from an iterator copying the data refered by the iterator. It is possible to specify the
*begin*and*end*of the iterator or the*begin*and the size. Note that for the Vector the iterator is not generic and must be of type*T*,*where T is the type of the contained elements.

Here are some examples on how to create a vector. In the following we assume that we are using the namespace ROOT::Math.

SVector>double,N> v; _// create a vector of size N, v[i]=0_

SVector>double,3> v(1,2,3); _// create a vector of size 3, v[0]=1,v[1]=2,v[2]=3_

double a[9] = {1,2,3,4,5,6,7,8,9}; _// input data_

SVector>double,9> v(a,9); _// create a vector using the a[] data_

The single vector elements can be set or retrieved using the *operator[i]* , *operator(i)* or the iterator interface. Notice that the index starts from zero and not from one as in FORTRAN. Also no check is performed on the passed index. Furthermore, all the matrix elements can be set also by using the ROOT::SVector::SetElements function passing a generic iterator. The elements can be accessed also by using the ROOT::Math::SVector::apply(i) function.

v[0] = 1; _ // set the first element _

v(1) = 2; _ // set the second element _

*(v.**begin**()+3) = 3; _// set the third element_

_// set vector elements from a std::vector<double>::iterator</double>_

std::vector <double>w(3);

v.SetElements(w.begin(),w.end());

double x = m(i); _// return the i-th element_

x = m.**apply**(i); _// return the i-th element_

x = *(m.**begin**()+i); _// return the i-th element

In addition there are methods to place a sub-vector in a vector. If the size of the the sub-vector is larger than the vector size a static assert ( a compilation error) is produced.

SVector>double,N> v;

SVector>double,M> w; _// M <= N otherwise a compilation error is obtained later _

_// place a vector of size M starting from element ioff, v[ioff + i] = w[i]_

v.**Place_at**(w,ioff);

_// return a sub-vector of size M starting from v[ioff]: w[i] = v[ioff + i]_

w = v.Sub < SVector>double,M> > (ioff);

For additional Vector functionality see the MatVecFunctions page