def METHOD_NAME(self, x):
    """Compute the unit normal at x for Neumann or Robin boundary conditions."""
    raise NotImplementedError(
        "{}.boundary_normal to be implemented".format(self.idstr)
    )