What does the term "regularity" mean? When I was an undergraduate, I took a course on regularity theory for nonlinear elliptic systems. This included topics such as the direct method of calculus of variations, mollifiers, integration over distributions, and the Calderón-Zygmund decomposition. 
In the introduction, there was a brief mention of regularity theory allowing us to change the space and conditions you are working in to make an equation work, or work with weaker conditions.  
Set $\Delta u=0$ in a domain $\Omega$. Integrating over the domain and introducing a term $\varphi$ we concluded that
$$
\int_{\Omega}u\Delta \varphi=0~\forall\varphi\in C^{\infty}_{c}(\Omega),~u\in L^1_{loc}(\Omega).
$$
This has changed the condition of $u$ being harmonic to $L^1$ integrable.
(Note: I am aware that what I have written is not completely mathematically rigorous. My notes are old and a bit rough).
I have also seen the term "regularity" related to differentiability and "regularity conditions" on initial data and "regularity control" on the solution, when referring to work on Strichartz and energy estimates. 
$\textbf{Question}$: Can someone shed some light on the term regularity?
I cannot find a formal definition of the term "regularity". I have pieced together information I have to get a rough idea of what it is.
I appreciate any thoughts or comments. Thank you in advance.
 A: It seems to me that a definition of "regularity" really depends on where it is being applied, for instance the regularity theory of weak solutions refers to the extra differentiability a function may have, due to it solving a PDE. The regularity a solution can inherit depends on the properties of the problem, i.e., the smoothness of a domain boundary, the smoothness of boundary conditions and initial data. Notice here that the smoothness of the initial data falls under the category of "regularity of initial data".
Let's consider the example you gave, let $u$ be a harmonic function, which I assume you mean in the classical sense, i.e., $u\in C^2(\Omega)\cap C(\overline{\Omega})$, and $$-\Delta u=0\quad\mbox{in}\quad\Omega,$$
even in this classical setting we can gain "regularity" by plugging a derivative of $u$ into the equation and noticing that $-\Delta$ commutes with the derivative. We then find that $u\in C^\infty(\Omega)-$it is in fact more regular than a general function in $u\in C^2(\Omega)\cap C(\overline{\Omega})$. If we wanted to extend this smoothness to the boundary, we would need $\partial\Omega$ to satisfy some smoothness conditions.
You mentioned that if we multiply the equation by a test function $\varphi\in C^\infty_c(\Omega)$ and integrate by parts, we obtain
$$\int_\Omega u\Delta\varphi=0\quad\forall\varphi\in C^\infty_c(\Omega),$$
which indeed tells us that $u\in L^1_{loc}(\Omega)$. But we already knew this, since  $u\in C^2(\Omega)\cap C(\overline{\Omega})$.
Regularity theory works the other way: (to make my job easier) take $u\in H^1(\Omega)$, and assume
$$\int_\Omega u\Delta\varphi=0\quad\forall\varphi\in C^\infty_c(\Omega),$$
this now tells us more about the nature of $u$. We can integrate by parts once to find
$$\int_\Omega \nabla u\cdot\nabla\varphi=0\quad\forall\varphi\in C^\infty_c(\Omega),$$
now by the choice of a suitable test function (using difference quotients), we can actually deduce that $u\in H^2_{loc}(\Omega)$, in which case we can justify integrating by parts once more to find
$$\int_\Omega \varphi\Delta u=0\quad\forall\varphi\in C^\infty_c(\Omega).$$
This holds a.e. and so $-\Delta u=0$ a.e.
We can bootstrap the difference quotients - test function method to deduce that $u\in H^k_{loc}(\Omega)$ for all $k\in\Bbb N$, i.e., $u\in C^\infty(\Omega)$. Again if the boundary is more regular, we can extend this notion of improved differentiability to the boundary.
Overall we have used properties of the Laplace equation to deduce that weak solutions behave in a nicer way than they "need to" to be considered a solution, i.e., we go from $H^1(\Omega)$ to $C^\infty(\overline{\Omega})$.
We see that the solution is more "regular", since really the word "regular" could be coined with the notion of "well behaved", and indeed in general a $C^\infty(\overline{\Omega})$ function is more well behaved than a $H^1(\Omega)$ function.
