Skip to content

simplify and cleanup the SolverFactory

Nils-Arne Dreier requested to merge simplify_solverfactory into master

This MR simplifies the SolverFactory. This includes the following changes:

  • Unification of the creator interface. All creators (solvers and preconditioners) now get the new type OperatorTraits and a shared_ptr<LinearOperator<X,Y>>. This has the advantage, that all available information is passed to the creator and can be accessed easily via the OperatorTraits. In particular, it is also possible for an iterative solver creator to access the matrix. This was not possible before as the actual operator type was not passed to the creator.
  • The OperatorTraits replace the TypeList that was passed to the creators before. This makes the code far more readable as the type can be accessed by their names instead of a list index.
  • Direct and iterative solvers are no longer kept in separate factories. The reason for keeping them separate was the need of the direct solver to access the matrix. This can now be achieved via the OperatorTraits.

Open questions regarding backward-compatibility:

  • Do we need backward compatibility to the old creator interface?
    • For the direct and iterative solvers, we could add easily the old macros DUNE_REGISTER_ITERARIVE_SOLVER and DUNE_REGISTER_DIRECT_SOLVER the wrap a provided old interface creator to the new interface.
    • I'm not sure how to deal with DUNE_REGISTER_PRECONDITIONER.
    • I added the old macros again and make OperatorTraits compatible to TypeList<matrix_tpye, domain_type, range_type> so that the old interface can be used. For the DUNE_REGISTER_ITERATIVE_SOLVER i added a wrapper that translates the creator to the new interface.
  • The method getSolverFromFactory had a third parameter before to inject a precreated preconditioner. This is now a bit more difficult as the preconditioner is created the the creator of the iterative solvers.
    • Possible solution: Add a setPrecondititioner method to the IterativeSolver class, this could then be used to replace the preconditioner of the created iterative solver. With that we could implement a function with the same functionality as before.
    • I solved that by defining the passed preconditioner in the ProconditionerRegistry with a dummy tag (__passed at runtime__)
Edited by Nils-Arne Dreier

Merge request reports

Loading