Lecture 7: Thread distributions MVM

Scaffold Head: #include "include/hpc_helpers.hpp" #include <iostream> #include <cstdint> #include <vector> #include <thread> template < typename value_t, typename index_t> void init( std::vector<value_t>& A, std::vector<value_t>& x, index_t m, index_t n) { for (index_t row = 0; row < m; row++) for (index_t col = 0; col < n; col++) A[row*n+col] = row >= col ? 1 : 0; for (index_t col = 0; col < m; col++) x[col] = col; } template < typename value_t, typename index_t> void sequential_mult( std::vector<value_t>& A, std::vector<value_t>& x, std::vector<value_t>& b, index_t m, index_t n) { for (index_t row = 0; row < m; row++) { value_t accum = value_t(0); for (index_t col = 0; col < n; col++) accum += A[row*n+col]*x[col]; b[row] = accum; } }

Tests