def _matvec(self, x_vec): """ (To be tested): dot operation provided for scipy.sparse.linalg wrapper of self.forward() """ x2 = numpy.reshape(x_vec, self.Nd, order='F') return self.forward(x2)