I am trying to solve a large-scale symmetric indefinite linear system:
$$Hx=B$$
Here, the matrix $H$ is the Hessian matrix of large neural networks. So, the system is large-scale, symmetric, and sparse.
My question is: which solver is suitable for my problem?
I did some surveys and found some general answers here and here. But they are too general that I can not find a specific answer.
Any suggestion is welcomed.