Cpp, Numerical Methods, Programming.
allocation here the cost of this evaluation is certainly blowing up the computational time. In case of Eigen, Blaze and Fastor, since they have their own implementation the call to norm(u - u_old) is fused in to a single fully unrolled AVX loop with further FMA instructions. Here is the full