- the toolbox la contains implementations for operators; +, -, *, +=, -=
Implementation is for the glas backend only at this stage. The
implementation of the BLAS backend will be the next step for performing
timing comparisons. All is for dense vectors currently.
The organisation of the include files can be improved, as well as the
doc html pages. These are not complete yet. The vector algorithms are
currently subdivided into two classes.
1) At assignment, one loop evaluates the whole expression
e.g. w = x + y + z is evaluated in one loop
2) At assignment, the expression is evaluated term by term, e.g.
w = x + y + z is split up as
w = x ; w += y ; w += z ;
norm_2 can only be computed for vector expressions that can be evaluated
in a single loop, otherwise you need a temporary vector to store
intermediate results, which glas does not create automatically. A
similar restriction holds for dot().