Уже который десяток раз накалываюсь на том, что ublas::norm_2(x) это вторая норма, а я почему-то жду от неё норму в квадрате, т.е. ublas::inner_prod(x,x).
кстати, не имел ли ты дела с библиотекой Eigen?
кстати, не имел ли ты дела с библиотекой Eigen?