If we have the Hilbert space $\mathcal H = L^2(\mathbb R^3, \mathbb C^4)$ and a Hamiltonian:
$$H=\gamma^i p_i + m \gamma^0$$
where $\gamma^i$ are matrices and $\{\gamma^i,\gamma^j\}=\delta^{ij}$. A statement I found in a book is that the Hilbert space decomposes into the orthogonal sum
$$\mathcal H = \mathcal H_+ \oplus \mathcal H_-$$
Where $\mathcal H_\pm$ is given by the positive/negative (including zero) eigenspaces of the Dirac Hamiltonian.
As far as I can tell, the solutions of the Dirac equation are not in $L^2$. Is this correct?
What is meant then with this decomposition?