Hello everyone! Roughly speaking, the question is - what are the possibilities to solve the linear system Ax = b with multiple righthand sides?
I am solving using PETSC a linear system Ax = b. Solution method is BiCGStab with my user-defined preconditioner. Can I somehow exploit the fact that I will need to solve a lot of systems with the same matrix but different righthand sides? I mean, can I simultaneously calculate all the solution vectors for all righthand sides like it is done, e.g., in PARDISO? As far as I understand, the only way is to call KSP_solve in a loop over all righthand sides, each time finding only one solution for the current righthand side (changing inside the loop). Am I right? I would greatly appreciate any help! Thanks in advance! Best regards, Kirill Voronin
