def _default_decomposable_op(self): probs = asarray(self.probs).reshape((1, -1)) # 2D array return (rbf_kernel(probs.T, gamma=self.gamma_quantile) if self.gamma_quantile != npinf else eye(len(self.probs)))