Gradient maximum principle for harmonic functions (${\sup }_{\Omega }|Dw|={\sup }_{\partial \Omega }|Dw|$)

For a harmonic function $w$ smooth inside a bounded domain and $C^1$ on its closure, $\Delta |Dw{|}^2=2|D^{2}w{|}^2\ge 0$. The weak maximum principle for elliptic operators and continuity of the gradient imply ${\sup }_{\Omega }|Dw|={\sup }_{\partial \Omega }|Dw|$.