ImsLinearSolver.f90

Go to the documentation of this file.

module imslinearsolvermodule
  use kindmodule, only: i4b, dp
  use constantsmodule, only: lensolutionname
  use linearsolverbasemodule
  use matrixbasemodule
  use vectorbasemodule
  use sparsematrixmodule
  use imslinearsettingsmodule
  use convergencesummarymodule
 
  implicit none
  private
 
  public :: create_ims_solver
 
  type, public, extends(linearsolverbasetype) :: imslinearsolvertype
  contains
    procedure :: initialize => ims_initialize
    procedure :: print_summary => ims_print_summary
    procedure :: solve => ims_solve
    procedure :: destroy => ims_destroy
 
    procedure :: create_matrix => ims_create_matrix
  end type
 
contains
 
  function create_ims_solver(sln_name) result(solver)
    class(linearsolverbasetype), pointer :: solver
    character(len=LENSOLUTIONNAME) :: sln_name
    ! local
    class(imslinearsolvertype), pointer :: ims_solver
 
    allocate (ims_solver)
    solver => ims_solver
    solver%name = sln_name
 
  end function create_ims_solver
 
  subroutine ims_initialize(this, matrix, linear_settings, convergence_summary)
    class(imslinearsolvertype) :: this
    class(matrixbasetype), pointer :: matrix
    type(imslinearsettingstype), pointer :: linear_settings
    type(convergencesummarytype), pointer :: convergence_summary
  end subroutine ims_initialize
 
  subroutine ims_print_summary(this)
    class(imslinearsolvertype) :: this
  end subroutine ims_print_summary
 
  subroutine ims_solve(this, kiter, rhs, x, cnvg_summary)
    class(imslinearsolvertype) :: this
    integer(I4B) :: kiter
    class(vectorbasetype), pointer :: rhs
    class(vectorbasetype), pointer :: x
    type(convergencesummarytype) :: cnvg_summary
  end subroutine ims_solve
 
  subroutine ims_destroy(this)
    class(imslinearsolvertype) :: this
  end subroutine ims_destroy
 
  function ims_create_matrix(this) result(matrix)
    class(imslinearsolvertype) :: this
    class(matrixbasetype), pointer :: matrix
    ! local
    class(sparsematrixtype), pointer :: ims_matrix
 
    allocate (ims_matrix)
    matrix => ims_matrix
 
  end function ims_create_matrix
 
end module imslinearsolvermodule