diff options
Diffstat (limited to 'bench/btl/libs/STL/STL_interface.hh')
| -rw-r--r-- | bench/btl/libs/STL/STL_interface.hh | 9 |
1 files changed, 9 insertions, 0 deletions
diff --git a/bench/btl/libs/STL/STL_interface.hh b/bench/btl/libs/STL/STL_interface.hh index 9dd9b8ed4..3958d4af5 100644 --- a/bench/btl/libs/STL/STL_interface.hh +++ b/bench/btl/libs/STL/STL_interface.hh @@ -146,6 +146,15 @@ public : X[j] += t2; } } + + static inline void syr2(gene_matrix & A, gene_vector & B, gene_vector & X, int N) + { + for (int j=0; j<N; ++j) + { + for (int i=j; i<N; ++i) + A[j][i] += B[i]*X[j] + B[j]*X[i]; + } + } static inline void atv_product(gene_matrix & A, gene_vector & B, gene_vector & X, int N) { |