Yeah, you are right there. But we and GPML have different optimization techniques and the difference is big (something about 1e-1). I plan add parameter's comparison (but little bit later), when problem with constrained parameters (i.e. non negative, etc) is fully solved, since for example l-bfgs and some other optimization algorithms has very bad results.