m_GlobalToLocal Module

private subroutine GlobalSegMapToIndices_(GSMap, comm, start, length)

Arguments

private subroutine GlobalSegMapToNavigator_(GSMap, comm, oNav)

Arguments

private subroutine GlobalSegMapToIndexArr_(GSMap, i_global, i_local, nindex, comm)

determine the global segments on this processor just once, so the info be used repeatedly below

this loop is optimized for the case that the indices in iglobal() are in the same order that they appear in the global segments, which seems usually (always?) to be the case.

note that the j loop exit condition is only executed when the index is not found in the current segment, which saves a factor of 2 since many consecutive indices are in the same segment.

!! this version vectorizes (outer loop) !! performance for in-order input is slightly slower than the above !! but performance on out-of-order input is probably much better !! at the moment we are going on the assumption that caller is !! likely providing in-order, so we won't use this version.

do i=1,nindex

i_local(i)= -1

SEARCH_LOOP: do j=1,nlseg

  if ( (mygs_lb(j) <= i_global(i)) .and. &
       (i_global(i) <= mygs_ub(j))) then
    i_local(i) = mygs_lstart(j) + (i_global(i) - mygs_lb(j))
  endif

end do SEARCH_LOOP

end do

Arguments

private function GlobalSegMapToIndex_(GSMap, i_g, comm)

Arguments

Return Value integer

private function GlobalMapToIndex_(GMap, i_g, comm)

Arguments

Return Value integer

private subroutine GlobalSegMapToLocalMatrix_(sMat, GSMap, RCFlag, comm)

rml new code from here down - do the mapping all in one function call which has been tuned for speed

Arguments