Hi,
I thinks there is a mistake in this documentation page
http://www.boost.org/doc/libs/1_35_0/libs/numeric/ublas/doc/products.htm
but not sure (for version 1.36 also).
Anyway the questions are :
1- how axpy_prod
(const vector_expression< E1 > &e1, const
matrix_expression< E2 > &e2, V &v, bool init=true) could
compute v += A x or v = A x and v +=
AT x or v = AT x knowing
that there is no params to indicate whether to transpose or not.
Shouldn't one of the axpy_prod be an opb_prod?
2 - axpy_prod
(const matrix_expression< E1 > &e1, const
matrix_expression< E2 > &e2, M &m, bool init=true) and
opb_prod
(const matrix_expression< E1 > &e1, const
matrix_expression< E2 > &e2, M &m, bool init=true) seem
to compute the same thing (based on the doc) M += A X or
M = A X. Shouldn't one of them compute M += AT
X or M = AT X
Can somebody lighten me?
Thanks
Best regards