|
13749 | 13749 | matrix_vector_product(scaled(3.0, A), x, y);
|
13750 | 13750 | }
|
13751 | 13751 |
|
13752 |
| -// y = 3.0 * A * x + 2.0 * y |
13753 |
| -void scaled_matvec_2(mdspan<double, extents<size_t, num_rows, num_cols>> A, |
13754 |
| - mdspan<double, extents<size_t, num_cols>> x, mdspan<double, extents<size_t, num_rows>> y) { |
13755 |
| - matrix_vector_product(scaled(3.0, A), x, scaled(2.0, y), y); |
13756 |
| -} |
13757 |
| - |
13758 | 13752 | // z = 7.0 times the transpose of A, times y
|
13759 | 13753 | void scaled_transposed_matvec(mdspan<double, extents<size_t, num_rows, num_cols>> A,
|
13760 | 13754 | mdspan<double, extents<size_t, num_rows>> y, mdspan<double, extents<size_t, num_cols>> z) {
|
|
13785 | 13779 | \tcode{z} may alias \tcode{y}.
|
13786 | 13780 | \end{itemdescr}
|
13787 | 13781 |
|
| 13782 | +\begin{example} |
| 13783 | +\begin{codeblock} |
| 13784 | +// y = 3.0 * A * x + 2.0 * y |
| 13785 | +void scaled_matvec_2(mdspan<double, extents<size_t, num_rows, num_cols>> A, |
| 13786 | + mdspan<double, extents<size_t, num_cols>> x, mdspan<double, extents<size_t, num_rows>> y) { |
| 13787 | + matrix_vector_product(scaled(3.0, A), x, scaled(2.0, y), y); |
| 13788 | +} |
| 13789 | +\end{codeblock} |
| 13790 | +\end{example} |
| 13791 | + |
13788 | 13792 | \rSec3[linalg.algs.blas2.symv]{Symmetric matrix-vector product}
|
13789 | 13793 |
|
13790 | 13794 | \pnum
|
|
0 commit comments