|
13664 | 13664 | matrix_vector_product(scaled(3.0, A), x, y);
|
13665 | 13665 | }
|
13666 | 13666 |
|
13667 |
| -// y = 3.0 * A * x + 2.0 * y |
13668 |
| -void scaled_matvec_2(mdspan<double, extents<size_t, num_rows, num_cols>> A, |
13669 |
| - mdspan<double, extents<size_t, num_cols>> x, mdspan<double, extents<size_t, num_rows>> y) { |
13670 |
| - matrix_vector_product(scaled(3.0, A), x, scaled(2.0, y), y); |
13671 |
| -} |
13672 |
| - |
13673 | 13667 | // z = 7.0 times the transpose of A, times y
|
13674 | 13668 | void scaled_transposed_matvec(mdspan<double, extents<size_t, num_rows, num_cols>> A,
|
13675 | 13669 | mdspan<double, extents<size_t, num_rows>> y, mdspan<double, extents<size_t, num_cols>> z) {
|
|
13700 | 13694 | \tcode{z} may alias \tcode{y}.
|
13701 | 13695 | \end{itemdescr}
|
13702 | 13696 |
|
| 13697 | +\begin{example} |
| 13698 | +\begin{codeblock} |
| 13699 | +// y = 3.0 * A * x + 2.0 * y |
| 13700 | +void scaled_matvec_2(mdspan<double, extents<size_t, num_rows, num_cols>> A, |
| 13701 | + mdspan<double, extents<size_t, num_cols>> x, mdspan<double, extents<size_t, num_rows>> y) { |
| 13702 | + matrix_vector_product(scaled(3.0, A), x, scaled(2.0, y), y); |
| 13703 | +} |
| 13704 | +\end{codeblock} |
| 13705 | +\end{example} |
| 13706 | + |
13703 | 13707 | \rSec3[linalg.algs.blas2.symv]{Symmetric matrix-vector product}
|
13704 | 13708 |
|
13705 | 13709 | \begin{note}
|
|
0 commit comments