const_m2s
scalar_mult(&product, const_m2s(m, i, j, rank), &a[j]);
scalar_mult(&product, const_m2s(m, j, i, rank), &a[j]);