implemented faster way of computing marginal likelihood
found a way to compute the random scaling matrix (at least it seems to work)
fixed a wrong call to the hadamard transform
extended tests
refactored tests
improved documentation
added another test