Suppose $\mathcal{X}$ and $\mathcal{Y}$ are two Hilbert spaces. Let $A:\mathcal{X} \mapsto \mathcal{Y}$ be a bounded linear operator. Consider a linear operator equation $Ax=b$. My question is what would be the set of all solutions of this equation provided that it is solvable.
When $A$ is an $m\times n$ real matrix, hence a linear operator from $\mathbb{R}^m$ to $\mathbb{R}^n$, we know that set of all roots of $Ax=b$, if it was solvable, is $$ \{x = A^+b + (I-A^+A)c,\:\text{for all $c\in\mathbb{R}^n$,}\} $$where $A^+$ is the generalized inverse, e.g. the Moore-Penrose pseudoinverse, of $A$.
Do we have similar conclusion for the general linear operator $A$?