Avoid direct inversion of covariance matrices for numerical stability

[coderabbitai]

As discussed in PR #3877 (comment), it was suggested that we should generally avoid direct inversion of covariance matrices in the ACTS codebase and use solvers like ldlt().solve() for better numerical stability. This improvement can be applied throughout the project.

[github-actions]

This issue/PR has been automatically marked as stale because it has not had recent activity. The stale label will be removed if any interaction occurs.