update documentation (taken from PhD thesis draft)
add hyper-diffusion collisional operator

Collapse of the CvK modes seems to be more sensitive
on higher order diffusivities. This module includes
higher order diffusion operators, so we can study
saturation/non-linear state.