For a positive definite symmetric matrix A, the linear equation
Ax=b
can be solved either by torch.cholesky_solve or by torch.solve. We might expect torch.cholesky_solve is faster because it uses the positive definite property. However, torch.cholesky_solve is slower than torch.solve. For a matrix of shape 128*128, the time consuming is:
torch.cholesky_solve  10.0ms
torch.solve  7.03ms
(NVIDIA GPU Titan X performance)
Why? So how to achieve the best performance? Thank you!