Solution (source code)

= Solution

For $u\in L^2_{\rm loc}(\Omega)$, the correct formulation is the <distributional identity>
$$
\int_\Omega u\,\Delta\varphi\,dx=0
\qquad\text{for every }\varphi\in C_c^\infty(\Omega).
$$
Thus $\Delta u=0$ as a <distribution>. The <Weyl lemma> applies already to locally integrable distributions, so $u$ agrees almost everywhere with a smooth harmonic function.