Pseudo should implement calculation for Hessian (as is written in the paper), which would also make it faster.