the zca in linear_decoder_exercise.py you writed is wrong

[TorosFanny]

You should add a np.sqrt above s+epsilon like below.

(u, s, vt) = np.linalg.svd(sigma) zca_white = u.dot(np.diag(1 / np.sqrt(s + epsilon))).dot(u.transpose())