VectorBase.f90

Go to the documentation of this file.

module vectorbasemodule
  use kindmodule, only: i4b, dp, lgp
  implicit none
  private
 
  type, public, abstract :: vectorbasetype
    logical(LGP) :: is_mem_managed
  contains
    procedure(create_mm_if), deferred :: create_mm
    procedure(destroy_if), deferred :: destroy
    procedure(get_array_if), deferred :: get_array
    procedure(get_ownership_range_if), deferred :: get_ownership_range
    procedure(get_size_if), deferred :: get_size
    procedure(get_value_local_if), deferred :: get_value_local
    procedure(zero_entries_if), deferred :: zero_entries
    procedure(set_value_local_if), deferred :: set_value_local
    procedure(axpy_if), deferred :: axpy
    procedure(norm2_if), deferred :: norm2
    procedure(print_if), deferred :: print
  end type vectorbasetype
 
  abstract interface
    subroutine create_mm_if(this, n, name, mem_path)
      import vectorbasetype, i4b
      class(vectorbasetype) :: this !< this vector
      integer(I4B) :: n !< the nr. of elements in the (local) vector
      character(len=*) :: name !< the variable name (for access through memory manager)
      character(len=*) :: mem_path !< memory path for storing the underlying memory items
    end subroutine create_mm_if
    subroutine destroy_if(this)
      import vectorbasetype
      class(vectorbasetype) :: this !< this vector
    end subroutine destroy_if
    function get_array_if(this) result(array)
      import vectorbasetype, dp
      class(vectorbasetype) :: this !< this vector
      real(dp), dimension(:), pointer, contiguous :: array
    end function get_array_if
    subroutine get_ownership_range_if(this, start, end)
      import vectorbasetype, i4b
      class(vectorbasetype) :: this !< this vector
      integer(I4B) :: start, end !< local range in global numbering
    end subroutine
    function get_size_if(this) result(size)
      import vectorbasetype, i4b
      class(vectorbasetype) :: this !< this vector
      integer(I4B) :: size !< the global vector size
    end function
    function get_value_local_if(this, idx) result(val)
      import vectorbasetype, i4b, dp
      class(vectorbasetype) :: this !< this vector
      integer(I4B) :: idx !< index in local numbering
      real(dp) :: val !< the value at the index
    end function
    subroutine zero_entries_if(this)
      import vectorbasetype
      class(vectorbasetype) :: this
    end subroutine
    subroutine set_value_local_if(this, idx, val)
      import vectorbasetype, i4b, dp
      class(vectorbasetype) :: this !< this vector
      integer(I4B) :: idx !< index in local numbering
      real(DP) :: val !< the value to set
    end subroutine
    subroutine axpy_if(this, alpha, vec_x)
      import vectorbasetype, dp
      class(vectorbasetype) :: this !< this vector y => alpha*x + y
      real(DP) :: alpha !< the factor
      class(vectorbasetype), pointer :: vec_x !< the vector to add
    end subroutine
    function norm2_if(this) result(n2)
      import vectorbasetype, dp
      class(vectorbasetype) :: this !< this vector
      real(dp) :: n2 !< the calculated 2-norm
    end function
    subroutine print_if(this)
      import vectorbasetype
      class(vectorbasetype) :: this
    end subroutine print_if
  end interface
 
end module vectorbasemodule