Went back to Hessian giving back a matrix of doubles. It is twice faster