49 views (last 30 days)

Show older comments

I have a sparse, banded matrix, A, for which I am solving the linear system Ax=b using A\b in Matlab.

The documentation for mldivide shows a flow chart that explains which solver it will use. However, I would like to automate the process of determining the best algorithm. Therefore, I would really just like to know which algorithm mldivide chooses for different inputs. Is there a way to get this information out?

Grzegorz Knor
on 19 Apr 2017

spparms('spumoni',2)

sparse(rand(3))\rand(3,1)

On my computer it produces:

sp\: bandwidth = 2+1+2.

sp\: is A diagonal? no.

sp\: is band density (1) > bandden (0.5) to try banded solver? yes.

sp\: is LAPACK's banded solver successful? yes.

Jan
on 27 Oct 2015

Edited: Jan
on 27 Oct 2015

I'm not sure if I understand your question. You explain:

The documentation for mldivide shows a flow chart that explains which solver it will use.

Is this not a valid answer of your question already? Or are you looking for the name of the library function instead of the name of the solver? Matlab does not offer the source code for mldivide.

Steven Lord
on 30 Oct 2015

Even if it DID tell you, for what purpose would you use that information?

Perhaps you would be better served by the LINSOLVE function, which allows you to specify various properties of the coefficient matrix for the system you're trying to solve. This can avoid some of the testing for those matrix properties (which for large matrices could be expensive.)

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

Start Hunting!