docs: def predict(self,Xstar):
          """
          predict on Xstar
          """
          KV = self._update_cache()
          self.covar.setXstar(Xstar)
          Kstar = self.covar.Kcross()
          Ystar = SP.dot(Kstar,KV['alpha'])
          return Ystar