I would like an example where the maximum principle fails in a bounded smooth domain $\Omega$ where one has a solution which is not $C^2(\Omega)$ to $Lu=0$ where $L$ is elliptic and linear. This obviously must rely on the coefficients being discontiuous for the elliptic operator since otherwise one can do interior regularity estimates. All of the examples I have tried to come up with turn out to not actually be weak solutions so I'm stuck on this. Perhaps maximum principles extend to the non-smooth setting? Any insight on this would be appreciated.