I read here that if a matrix is positive definite then it's invertible.
If I take some vector nx1 and multiply it by another one that's 1xn then I get an nxn matrix. But It doesn't seem to be invertible according to a short sympy script I wrote.
a1, a2, a3, a4, a5, a6 = sp.symbols('a1 a2 a3 a4 a5 a6')
va = sp.Matrix([a1, a2, a3, a4, a5, a6])
pos_def_matrix = va*va.T
pos_def_matrix.inv()
which gives me the error.
Matrix det == 0; not invertible.
So am I misunderstanding something?
————————————
The operation performed in a paper I’m trying to understand is this
$w_t = E_t[R^e_{t+1} {R^e_{t+1}}^T]^{-1} E_t[R^e_{t+1}]$
But it seems we showed $R^e_{t+1} {R^e_{t+1}}^T$ is not possible.