doing ``log(det(...))`` is prone to overflow / underflow. Julia already had an ``logdet`` function that can compute the logarithm of determinant of a positive definite matrix.