virtualexchangemodule Module Reference

Data Types
type	virtualexchangetype
	The Virtual Exchange is based on two Virtual Models and is therefore not always strictly local or remote. We have to consider three different cases: More...

Functions/Subroutines
subroutine	vx_create (this, name, exg_id, m1_id, m2_id)
	Create the virtual exchange base. More...
subroutine	init_virtual_data (this)
subroutine	vx_prepare_stage (this, stage)
subroutine	vx_get_recv_items (this, stage, rank, virtual_items)
subroutine	vx_get_send_items (this, stage, rank, virtual_items)
logical(lgp) function	vx_has_mover (this)
	Checks if there is an active mover in the exchange. More...
subroutine	vx_destroy (this)
subroutine	allocate_data (this)
subroutine	deallocate_data (this)
class(virtualexchangetype) function, pointer, public	get_virtual_exchange (exg_id)
	Returns a virtual exchange with the specified id. More...
class(virtualexchangetype) function, pointer, public	get_virtual_exchange_from_list (list, idx)
class(virtualexchangetype) function, pointer, private	cast_as_virtual_exchange (obj_ptr)

Function/Subroutine Documentation

allocate_data()

subroutine virtualexchangemodule::allocate_data ( class(virtualexchangetype)  this )

private

Definition at line 263 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
 
    allocate (this%nexg)
    allocate (this%naux)
    allocate (this%ianglex)
    allocate (this%nodem1)
    allocate (this%nodem2)
    allocate (this%ihc)
    allocate (this%cl1)
    allocate (this%cl2)
    allocate (this%hwva)
    allocate (this%auxvar)
 

cast_as_virtual_exchange()

class(virtualexchangetype) function, pointer, private virtualexchangemodule::cast_as_virtual_exchange ( class(*), pointer  obj_ptr )

private

Definition at line 331 of file VirtualExchange.f90.

    class(*), pointer :: obj_ptr
    class(VirtualExchangeType), pointer :: virtual_exg
 
    virtual_exg => null()
    select type (obj_ptr)
    class is (virtualexchangetype)
      virtual_exg => obj_ptr
    end select
 

Here is the caller graph for this function:

deallocate_data()

subroutine virtualexchangemodule::deallocate_data ( class(virtualexchangetype)  this )

private

Definition at line 279 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
 
    deallocate (this%nexg)
    deallocate (this%naux)
    deallocate (this%ianglex)
    deallocate (this%nodem1)
    deallocate (this%nodem2)
    deallocate (this%ihc)
    deallocate (this%cl1)
    deallocate (this%cl2)
    deallocate (this%hwva)
    deallocate (this%auxvar)
 

get_virtual_exchange()

class(virtualexchangetype) function, pointer, public virtualexchangemodule::get_virtual_exchange

(

integer(i4b)

exg_id

)

Definition at line 297 of file VirtualExchange.f90.

    use virtualdatalistsmodule, only: virtual_exchange_list
    integer(I4B) :: exg_id
    class(VirtualExchangeType), pointer :: virtual_exg
    ! local
    integer(I4B) :: i
    class(*), pointer :: ve
 
    virtual_exg => null()
    do i = 1, virtual_exchange_list%Count()
      ve => virtual_exchange_list%GetItem(i)
      select type (ve)
      class is (virtualexchangetype)
        if (ve%id == exg_id) then
          virtual_exg => ve
          return
        end if
      end select
    end do
 

Here is the caller graph for this function:

get_virtual_exchange_from_list()

class(virtualexchangetype) function, pointer, public virtualexchangemodule::get_virtual_exchange_from_list	(	type(listtype)	list,
		integer(i4b)	idx
	)

Definition at line 319 of file VirtualExchange.f90.

    type(ListType) :: list
    integer(I4B) :: idx
    class(VirtualExchangeType), pointer :: virtual_exg
    ! local
    class(*), pointer :: obj_ptr
 
    obj_ptr => list%GetItem(idx)
    virtual_exg => cast_as_virtual_exchange(obj_ptr)
 

Here is the call graph for this function:

Here is the caller graph for this function:

init_virtual_data()

subroutine virtualexchangemodule::init_virtual_data ( class(virtualexchangetype)  this )

private

Definition at line 120 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    ! local
    logical(LGP) :: is_nodem1_local
    logical(LGP) :: is_nodem2_local
 
    ! exchanges can be hybrid with both local and remote
    ! fields, nodem1/2 array only local when corresponding
    ! model sits on the same process
    is_nodem1_local = this%v_model1%is_local
    is_nodem2_local = this%v_model2%is_local
    call this%set(this%nexg%base(), 'NEXG', '', map_all_type)
    call this%set(this%naux%base(), 'NAUX', '', map_all_type)
    call this%set(this%ianglex%base(), 'IANGLEX', '', map_all_type)
    call this%set(this%nodem1%base(), 'NODEM1', '', &
                  map_all_type, is_nodem1_local)
    call this%set(this%nodem2%base(), 'NODEM2', '', &
                  map_all_type, is_nodem2_local)
    call this%set(this%ihc%base(), 'IHC', '', map_all_type)
    call this%set(this%cl1%base(), 'CL1', '', map_all_type)
    call this%set(this%cl2%base(), 'CL2', '', map_all_type)
    call this%set(this%hwva%base(), 'HWVA', '', map_all_type)
    call this%set(this%auxvar%base(), 'AUXVAR', '', map_all_type)
 

vx_create()

subroutine virtualexchangemodule::vx_create	(	class(virtualexchangetype)	this,
		character(len=*)	name,
		integer(i4b)	exg_id,
		integer(i4b)	m1_id,
		integer(i4b)	m2_id
	)

Definition at line 97 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    character(len=*) :: name
    integer(I4B) :: exg_id
    integer(I4B) :: m1_id
    integer(I4B) :: m2_id
    ! local
    logical(LGP) :: is_local
 
    this%v_model1 => get_virtual_model(m1_id)
    this%v_model2 => get_virtual_model(m2_id)
 
    ! 1) both models local:  is_local = true
    ! 2) only one of them:   is_local = true
    ! 3) both models remote: is_local = false
    is_local = this%v_model1%is_local .or. this%v_model2%is_local
    call this%VirtualDataContainerType%vdc_create(name, exg_id, is_local)
 
    call this%allocate_data()
    call this%init_virtual_data()
 

vx_destroy()

subroutine virtualexchangemodule::vx_destroy ( class(virtualexchangetype)  this )

private

Definition at line 255 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
 
    call this%VirtualDataContainerType%destroy()
    call this%deallocate_data()
 

vx_get_recv_items()

subroutine virtualexchangemodule::vx_get_recv_items ( class(virtualexchangetype)  this, integer(i4b)  stage, integer(i4b)  rank, type(stlvecint)  virtual_items  )

private

Definition at line 176 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    integer(I4B) :: stage
    integer(I4B) :: rank
    type(STLVecInt) :: virtual_items
    ! local
    integer(I4B) :: nodem1_idx, nodem2_idx
    class(*), pointer :: vdi
 
    vdi => this%nodem1
    nodem1_idx = this%virtual_data_list%GetIndex(vdi)
    vdi => this%nodem2
    nodem2_idx = this%virtual_data_list%GetIndex(vdi)
 
    if (this%v_model1%is_local .and. &
        this%v_model2%orig_rank == rank) then
      ! this is our dual exchange on the other rank,
      ! only receive nodem2
      if (this%nodem2%check_stage(stage)) then
        call virtual_items%push_back(nodem2_idx)
      end if
    else if (this%v_model2%is_local .and. &
             this%v_model1%orig_rank == rank) then
      ! the reverse case...
      if (this%nodem1%check_stage(stage)) then
        call virtual_items%push_back(nodem1_idx)
      end if
    else
      ! receive all using base
      call this%VirtualDataContainerType%get_recv_items(stage, rank, &
                                                        virtual_items)
    end if
 

vx_get_send_items()

subroutine virtualexchangemodule::vx_get_send_items ( class(virtualexchangetype)  this, integer(i4b)  stage, integer(i4b)  rank, type(stlvecint)  virtual_items  )

private

Definition at line 211 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    integer(I4B) :: stage
    integer(I4B) :: rank
    type(STLVecInt) :: virtual_items
    ! local
    integer(I4B) :: nodem1_idx, nodem2_idx
    class(*), pointer :: vdi
 
    vdi => this%nodem1
    nodem1_idx = this%virtual_data_list%GetIndex(vdi)
    vdi => this%nodem2
    nodem2_idx = this%virtual_data_list%GetIndex(vdi)
    if (this%v_model1%is_local .and. &
        this%v_model2%orig_rank == rank) then
      ! this is our dual exchange on the other rank,
      ! only send nodem1
      if (this%nodem1%check_stage(stage)) then
        call virtual_items%push_back(nodem1_idx)
      end if
    else if (this%v_model2%is_local .and. &
             this%v_model1%orig_rank == rank) then
      ! the reverse case...
      if (this%nodem2%check_stage(stage)) then
        call virtual_items%push_back(nodem2_idx)
      end if
    else
      ! send all of it
      call this%VirtualDataContainerType%get_send_items(stage, rank, &
                                                        virtual_items)
    end if
 

vx_has_mover()

logical(lgp) function virtualexchangemodule::vx_has_mover ( class(virtualexchangetype)  this )

private

Definition at line 247 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    logical(LGP) :: has_mover
 
    has_mover = .false.
 

vx_prepare_stage()

subroutine virtualexchangemodule::vx_prepare_stage ( class(virtualexchangetype)  this, integer(i4b)  stage  )

private

Definition at line 146 of file VirtualExchange.f90.

    class(VirtualExchangeType) :: this
    integer(I4B) :: stage
    ! local
    integer(I4B) :: nexg, naux
 
    if (stage == stg_aft_exg_df) then
 
      call this%map(this%nexg%base(), (/stg_aft_exg_df/))
      call this%map(this%naux%base(), (/stg_aft_exg_df/))
      call this%map(this%ianglex%base(), (/stg_aft_exg_df/))
 
    else if (stage == stg_aft_con_cr) then
 
      nexg = this%nexg%get()
      naux = this%naux%get()
      call this%map(this%nodem1%base(), nexg, (/stg_aft_con_cr, &
                                                stg_bfr_con_df/))
      call this%map(this%nodem2%base(), nexg, (/stg_aft_con_cr, &
                                                stg_bfr_con_df/))
      call this%map(this%ihc%base(), nexg, (/stg_aft_con_cr/))
      call this%map(this%cl1%base(), nexg, (/stg_aft_con_cr/))
      call this%map(this%cl2%base(), nexg, (/stg_aft_con_cr/))
      call this%map(this%hwva%base(), nexg, (/stg_aft_con_cr/))
      call this%map(this%auxvar%base(), naux, nexg, (/stg_aft_con_cr/))
 
    end if