If the problem is definite, maybe you should be using shift-invert? (The Julia libraries won’t switch to shift-invert automatically.)
It’s hard to give more concrete advice without knowing your matrices.
It’s also not clear what is “not working” — the residual norms are small, so what is incorrect about the answers?