def fwd(self,X): 'forward linear operator' assert np.iscomplexobj(X) == self.cpx,'wrong value for cpx in constructor' return np.einsum('...jk,mj->...mk',X,self.S)