docs: def score(self, p, x, y, sample_weight=None):
          """
          Overloading default regression score method
  
          Parameters
          ----------
          p: ndarray
              (Ns, d) array of predictor variables (Ns samples, d dimensions)
              for regression
          x: ndarray
              ndarray of (x, y) points. Needs to be a (Ns, 2) array
              corresponding to the lon/lat, for example.
              array of Points, (x, y, z) pairs of shape (N, 3) for 3d kriging
          y: ndarray
              array of targets (Ns, )
          """
  
          return r2_score(y_pred=self.predict(p, x),
                          y_true=y,
                          sample_weight=sample_weight)