code: def logpdf_sum(self, f, y, Y_metadata=None):
          
          return np.sum(self.logpdf(f, y, Y_metadata=Y_metadata))
docs: Convenience function that can overridden for functions where this could
          be computed more efficiently