Cpp, Numerical Methods, Programming.
XTensor are unfortunately not performing really well at this benchmark. Looking closer it seems like the bottleneck in their case is the norm function which requires a LAPACK call on the expression norm(u - u_old) // armadillo xt::linalg::norm(u - u_old) // xtensor I don’t know how Armadillo and