subroutine alloc_buffers_srl(self, nwgc, &
ixS_srl_min1, ixS_srl_max1, &
ixS_srl_min2, ixS_srl_max2, &
ixS_srl_min3, ixS_srl_max3, &
ixR_srl_min1, ixR_srl_max1, &
ixR_srl_min2, ixR_srl_max2, &
ixR_srl_min3, ixR_srl_max3)
class(nbprocs_info_t) :: self
integer, intent(in) :: nwgc
integer, dimension(-1:2,-1:1), intent(in) :: &
ixS_srl_min1, ixS_srl_max1, &
ixS_srl_min2, ixS_srl_max2, &
ixS_srl_min3, ixS_srl_max3, &
ixR_srl_min1, ixR_srl_max1, &
ixR_srl_min2, ixR_srl_max2, &
ixR_srl_min3, ixR_srl_max3
! .. local ..
integer :: inb, igrid, isize_S, isize_R
integer :: n_i1, n_i2, n_i3, i1, i2, i3
integer :: ixSmin1, ixSmin2, ixSmin3
integer :: ixRmin1, ixRmin2, ixRmin3
integer :: ixSmax1, ixSmax2, ixSmax3
integer :: ixRmax1, ixRmax2, ixRmax3
integer :: iib1, iib2, iib3, ibuf_start
self%imaxigrids_srl = 0
do inb = 1, self%nbprocs_srl
self%imaxigrids_srl = max(self%srl_nb(inb)%info%nigrids, self%imaxigrids_srl)
isize_S = 0; isize_R = 0; ibuf_start = 1
do igrid = 1, self%srl_nb(inb)%info%nigrids
iib1 = idphyb(1,self%srl_nb(inb)%info%igrid(igrid))
iib2 = idphyb(2,self%srl_nb(inb)%info%igrid(igrid))
iib3 = idphyb(3,self%srl_nb(inb)%info%igrid(igrid))
call self%idecode( i1, i2, i3, self%srl_nb(inb)%info%iencode(igrid) )
n_i1=-i1; n_i2=-i2; n_i3=-i3
ixSmin1=ixS_srl_min1(iib1,i1);ixSmin2=ixS_srl_min2(iib2,i2)
ixSmin3=ixS_srl_min3(iib3,i3);ixSmax1=ixS_srl_max1(iib1,i1)
ixSmax2=ixS_srl_max2(iib2,i2);ixSmax3=ixS_srl_max3(iib3,i3)
ixRmin1=ixR_srl_min1(iib1,n_i1);ixRmin2=ixR_srl_min2(iib2,n_i2)
ixRmin3=ixR_srl_min3(iib3,n_i3);ixRmax1=ixR_srl_max1(iib1,n_i1)
ixRmax2=ixR_srl_max2(iib2,n_i2);ixRmax3=ixR_srl_max3(iib3,n_i3)
isize_S = isize_S + (ixSmax1-ixSmin1+1) * (ixSmax2-ixSmin2+1) * (ixSmax3-ixSmin3+1)
isize_R = isize_R + (ixRmax1-ixRmin1+1) * (ixRmax2-ixRmin2+1) * (ixRmax3-ixRmin3+1)
self%srl_nb(inb)%info%ibuf_start(igrid) = ibuf_start
self%srl_nb(inb)%info%isize(igrid) = (ixSmax1-ixSmin1+1) * (ixSmax2-ixSmin2+1) * (ixSmax3-ixSmin3+1) * nwgc
ibuf_start = ibuf_start + self%srl_nb(inb)%info%isize(igrid)
end do
self%srl_nb(inb)%rcv%size = isize_R * nwgc
self%srl_nb(inb)%send%size = isize_S * nwgc
self%srl_nb(inb)%info_send%size = 3 * self%srl_nb(inb)%info%nigrids
self%srl_nb(inb)%info_rcv%size = 3 * self%srl_nb(inb)%info%nigrids
allocate(self%srl_nb(inb)%rcv%buffer(isize_R * nwgc))
allocate(self%srl_nb(inb)%send%buffer(isize_S * nwgc))
allocate(self%srl_nb(inb)%info_rcv%buffer(3 * self%srl_nb(inb)%info%nigrids))
allocate(self%srl_nb(inb)%info_send%buffer(3 * self%srl_nb(inb)%info%nigrids))
end do
end subroutine alloc_buffers_srl