2
$\begingroup$

I would like to know what is the default options of LinearSolve when using LinearSolve[A,b], where A is a sparse matrix. There are several methods available such as Cholesky and Krysol etc., see here.

The reason I ask this is that I have a problem involving instabilities, the matrix becomes almost singular when close to the unstable point. The LinearSolve can find the converged solution with Newton iteration. I try to solve the exact same problem in Fenics, in which I can tune the solver parameters. The available solver parameters in Fenics are listed here. I tried with various solver parameter sets but can't get the converged solution at the unstable point.

I would like to know if I can set the same solver parameters in Fenics with that in Mathematica.

Improve this question
$\endgroup$
1
  • $\begingroup$ Very hard to answer without a concrete example, including code you tried that did not perform well. $\endgroup$
    – Daniel Lichtblau
    Commented Nov 26, 2017 at 16:49

1 Answer 1

Reset to default
1
$\begingroup$

The standard is Method -> "Multifrontal" which employs LU-factorization provided by UMFPACK.

Improve this answer
$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.

Not the answer you're looking for? Browse other questions tagged or ask your own question.