error: user defined reduction not found for ‘norm_b_square’
42 | #pragma omp for reduction(+:norm_b_square)
На что он ругается?
norm_b_square, а в цикле работаешь с *norm_b_square. А если ты имел в виду одноимённый norm_b_square из функции main, то нет, нельзя так сослаться на автоматическую переменную из другой функции.