MatrixBase.f90

Go to the documentation of this file.

module matrixbasemodule
  use constantsmodule, only: lenmempath
  use kindmodule, only: i4b, dp
  use sparsemodule, only: sparsematrix
  use vectorbasemodule
  implicit none
  private
 
  type, public, abstract :: matrixbasetype
    character(len=LENMEMPATH) :: memory_path
  contains
    procedure(init_if), deferred :: init
    procedure(destroy_if), deferred :: destroy
    procedure(create_vec_mm_if), deferred :: create_vec_mm
    procedure(create_vec_if), deferred :: create_vec
 
    procedure(get_value_pos_if), deferred :: get_value_pos
    procedure(get_diag_value_if), deferred :: get_diag_value
 
    procedure(set_diag_value_if), deferred :: set_diag_value
    procedure(set_value_pos_if), deferred :: set_value_pos
    procedure(add_value_pos_if), deferred :: add_value_pos
    procedure(add_diag_value_if), deferred :: add_diag_value
    procedure(zero_entries_if), deferred :: zero_entries
    procedure(zero_row_offdiag_if), deferred :: zero_row_offdiag
 
    procedure(get_first_col_pos_if), deferred :: get_first_col_pos
    procedure(get_last_col_pos_if), deferred :: get_last_col_pos
    procedure(get_column_if), deferred :: get_column
    procedure(get_position_if), deferred :: get_position
    procedure(get_position_diag_if), deferred :: get_position_diag
 
    procedure(get_aij_if), deferred :: get_aij
    procedure(get_row_offset_if), deferred :: get_row_offset
 
    procedure(multiply_if), deferred :: multiply
 
  end type matrixbasetype
 
  abstract interface
    subroutine init_if(this, sparse, mem_path)
      import matrixbasetype, sparsematrix
      class(matrixbasetype) :: this
      type(sparsematrix) :: sparse
      character(len=*) :: mem_path
    end subroutine
    subroutine destroy_if(this)
      import matrixbasetype
      class(matrixbasetype) :: this
    end subroutine
    function create_vec_mm_if(this, n, name, mem_path) result(vec)
      import matrixbasetype, vectorbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: n
      character(len=*) :: name
      character(len=*) :: mem_path
      class(vectorbasetype), pointer :: vec
    end function
    function create_vec_if(this, n) result(vec)
      import matrixbasetype, vectorbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: n
      class(vectorbasetype), pointer :: vec
    end function
 
    function get_value_pos_if(this, ipos) result(value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: ipos
      real(dp) :: value
    end function
    function get_diag_value_if(this, irow) result(diag_value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      real(dp) :: diag_value
    end function
 
    subroutine set_diag_value_if(this, irow, diag_value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      real(DP) :: diag_value
    end subroutine
    subroutine set_value_pos_if(this, ipos, value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: ipos
      real(DP) :: value
    end subroutine
    subroutine add_value_pos_if(this, ipos, value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: ipos
      real(DP) :: value
    end subroutine
    subroutine add_diag_value_if(this, irow, value)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      real(DP) :: value
    end subroutine
    subroutine zero_entries_if(this)
      import matrixbasetype
      class(matrixbasetype) :: this
    end subroutine
    subroutine zero_row_offdiag_if(this, irow)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: irow
    end subroutine
    function get_first_col_pos_if(this, irow) result(first_col_pos)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      integer(I4B) :: first_col_pos
    end function
    function get_last_col_pos_if(this, irow) result(last_col_pos)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      integer(I4B) :: last_col_pos
    end function
    function get_column_if(this, ipos) result(icol)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: ipos
      integer(I4B) :: icol
    end function
    function get_position_if(this, irow, icol) result(ipos)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      integer(I4B) :: icol
      integer(I4B) :: ipos
    end function
    function get_position_diag_if(this, irow) result(ipos_diag)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: irow
      integer(I4B) :: ipos_diag
    end function
    subroutine get_aij_if(this, ia, ja, amat)
      import matrixbasetype, i4b, dp
      class(matrixbasetype) :: this
      integer(I4B), dimension(:), pointer, contiguous :: ia
      integer(I4B), dimension(:), pointer, contiguous :: ja
      real(DP), dimension(:), pointer, contiguous :: amat
    end subroutine
    function get_row_offset_if(this) result(offset)
      import matrixbasetype, i4b
      class(matrixbasetype) :: this
      integer(I4B) :: offset
    end function
    subroutine multiply_if(this, vec_x, vec_y)
      import matrixbasetype, vectorbasetype
      class(matrixbasetype) :: this
      class(vectorbasetype), pointer :: vec_x
      class(vectorbasetype), pointer :: vec_y
    end subroutine
  end interface
 
end module matrixbasemodule