Implemented Gibbs kernel for arbitrary l up to first derivatives. Expression is now hard-coded (and hence much faster).