complete loglik (axis)

[tlienart]

follow the example set for loglik(g::GaussNP, ...)

[tlienart]

Should all accept abstract arrays, apart maybe from gradloglik (?)

[tlienart]

also complete uloglik

[tlienart]

Now in separate file gaussian_loglik, ultimately may want to restore a bit of structure there for potential other families