Should it be lmd * np.linalg.inv(gmat) as in the paper？

[Sunzhiliang]

probreg/probreg/bcpd.py

Line 118 in 2656589

sigma_mat_inv = lmd * gmat + s2s2 * np.diag(nu)

If we follow the paper,it shoud be
sigma_mat_inv = lmd * gmat + s2s2 * np.diag(nu).
However, if it is changed,the result is wrong.
Is it wrong to use "np.linalg.inv(gmat)" in the paper?

[neka-nat]

Thank you!
I've pushed fixes #48 and this issue to master.
Fixing both is working fine for the bunny.