different methods for log determinants

log_det

import numpy as np

# generate data
X = np.random.randn(1000, 10000)
K = X.dot(X.T)

# np function
np.linalg.slogdet(K)

# using eigh
S1, U1 = np.linalg.eigh(K)
np.log(S1[S1 > 1e-9]).sum()

# using svd
U,S,V = np.linalg.svd(X, full_matrices=False)
S2 = S*S
np.log(S2[S2 > 1e-9]).sum()