Ax = b for a lot of different b‘s

2 views (last 30 days)
Hello,
I have a system of equation Ax = b, where A is very large (n = 3e5 - 8e5) and sparse coefficient matrix, b is source term and x is the solution. Of course it's not a problem to solve such a problem in MATLAB, however I have to solve Ax = b for a lot of different b‘s without any change in the coefficient matrix A. The inversion of the matrix A take so long, that it is faster to solve the problem for each new b vector separately using an iterative solver. Unfortunately it is still way too slow, so I'm looking for a way how to speed it up.
Does anybody know, how to improve the performance?
Thank you for help!

Accepted Answer

Matt J
Matt J on 18 Apr 2013
Edited: Matt J on 18 Apr 2013
Neither direct inversion of A, nor iterative methods are required. Just create a matrix B whose columns are the different b. Then do
X=A\B
Each column X(:,i) will be the solution for the corresponding b.
  4 Comments
Matt J
Matt J on 19 Apr 2013
Edited: Matt J on 19 Apr 2013
well it really helps! I have lost some control, e.g. precision, parralelization or preconditioning
Not sure why you feel you've lost precision or why you think pre-conditioning matters when you're using a non-iterative method. It's true you have no control over the parallelization, but parallelization is done for you internally by MLDIVIDE.
David Kusnirak
David Kusnirak on 19 Apr 2013
Edited: David Kusnirak on 19 Apr 2013
It works just fine as I said, thanks for help. The precision is far behind my tolerance treshold I used to use with the iterative solver, and of course preconditioning is not necessary now. I just felt better to have more response about what's going on, as the mldivide is somehow just a 'blackbox' for me.

Sign in to comment.

More Answers (0)

Categories

Find more on Programming in Help Center and File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!