This should be` np.dot(self.cholsigmainv, X)` as before ?
(we are not using np.matrices, and * is always elementwise.