LinearSolverFactory.F90

Go to the documentation of this file.

module linearsolverfactory
  use simmodule, only: ustop
  use linearsolverbasemodule
  use matrixbasemodule
  use sparsematrixmodule
  use vectorbasemodule
 
  use imslinearsolvermodule
#if defined(__WITH_PETSC__)
  use petscsolvermodule, only: create_petsc_solver
#endif
 
  implicit none
  private
 
  public :: create_linear_solver
 
contains
 
  !> @brief Factory method to create the linear solver object
  !<
  function create_linear_solver(solver_mode, sln_name) result(solver)
    character(len=*) :: solver_mode
    character(len=*) :: sln_name
    class(linearsolverbasetype), pointer :: solver
 
    solver => null()
 
    if (solver_mode == 'IMS') then
      solver => create_ims_solver(sln_name)
      return
#if defined(__WITH_PETSC__)
    else if (solver_mode == 'PETSC') then
      solver => create_petsc_solver(sln_name)
#endif
    else
      call ustop('Unsupported solver mode: '//trim(solver_mode))
    end if
 
  end function create_linear_solver
 
end module linearsolverfactory