Specifying the KSP method to use in solve(a == L)
Hi,
I just spent an hour debugging why
solve(a == L, u, solver_
was using an LU decomposition.
It appears that the interface has changed. Now you're supposed to do solver_
If so,
a) how am I supposed to use another option PETSc provides (like bicg, tfqmr, etc)? It appears this logic only allows me to choose between cg and gmres, which is very restrictive.
b) the snippet of code above used to work in 1.0.0. Shouldn't there be an error message, rather than silently using LU instead of doing what I told it to? Not everyone will be willing to gdb in to see why it's not doing what they expect.
c) is this redesign discussed anywhere? I strongly prefer the old interface, and would like to see why it was changed.
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- DOLFIN Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Patrick Farrell for more information if necessary.