docs: def rmse(y, p):
      """Root Mean Squared Error (RMSE).
  
      Args:
          y (numpy.array): target
          p (numpy.array): prediction
  
      Returns:
          e (numpy.float64): RMSE
      """
  
      # check and get number of samples
      assert y.shape == p.shape
  
      return np.sqrt(mse(y, p))