抛弃 numpy.linalg.det ，换用 slogdet

numpy.linalg.det 的实现是先进行 LU 分解，再把 U 的对角线绝对值取 log ，加起来。得到 slogdet。

而我们在用 np.log(np.linalg.det(np.dot(beta, fsig2s_inv))) 终究还是取了 log ，不如直接调用 slogdet。

The determinant is computed via LU factorization using the LAPACK
routine ``z/dgetrf``.