In my current understanding, matrix formulation and wave-function formulation of QM are basically the same because $\left|\psi\right>$ and $\psi(x)$ are really the same mathematical object: A vector in the (vector) space of complex functions. My issue with this is the status of $\hat{x}$ and its eigenvectors $\left|x\right>$.
Let’s ask a simple question: What is the dimensionality of the vector space? We can simply write that every vector can be decomposed in a sum of energy eigenstates:
$\left|\psi\right> = \sum E_n\left|n\right>$
Or every vector can be decomposed in a sum of position eigenstates:
$\left|\psi\right> = \int \psi(x)\left|x\right>dx$
In the first case, we have an infinite enumerable basis. In the second case, we have an infinite non-enumerable basis. But they are supposed to be two different orthonormal basis of the same vector space!
Or put another way, how can we write $\psi(x) = \left<\psi|x\right>$ if the set of $\left|x\right>$s is enumerable but the domain of $\psi(x)$ isn’t ?