alloc_buffers_f Subroutine

public subroutine alloc_buffers_f(self, nwgc, ixS_p_min1, ixS_p_max1, ixS_p_min2, ixS_p_max2, ixS_p_min3, ixS_p_max3, ixR_r_min1, ixR_r_max1, ixR_r_min2, ixR_r_max2, ixR_r_min3, ixR_r_max3)

Type Bound

nbprocs_info_t

Arguments

Type IntentOptional Attributes Name
class(nbprocs_info_t) :: self
integer, intent(in) :: nwgc
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_min1
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_max1
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_min2
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_max2
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_min3
integer, intent(in), dimension(-1:1,0:3) :: ixS_p_max3
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_min1
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_max1
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_min2
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_max2
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_min3
integer, intent(in), dimension(-1:1,0:3) :: ixR_r_max3

Called by

proc~~alloc_buffers_f~~CalledByGraph proc~alloc_buffers_f nbprocs_info_t%alloc_buffers_f proc~build_connectivity build_connectivity proc~build_connectivity->proc~alloc_buffers_f proc~amr_coarsen_refine amr_coarsen_refine proc~amr_coarsen_refine->proc~build_connectivity proc~initlevelone initlevelone proc~initlevelone->proc~build_connectivity proc~read_forest read_forest proc~read_forest->proc~build_connectivity proc~main main proc~main->proc~initlevelone proc~read_snapshot read_snapshot proc~main->proc~read_snapshot proc~resettree_convert resettree_convert proc~main->proc~resettree_convert proc~settree settree proc~main->proc~settree proc~timeintegration timeintegration proc~main->proc~timeintegration proc~read_snapshot->proc~read_forest proc~read_snapshot_old read_snapshot_old proc~read_snapshot->proc~read_snapshot_old proc~read_snapshot_old->proc~read_forest proc~resettree resettree proc~resettree->proc~amr_coarsen_refine proc~resettree_convert->proc~amr_coarsen_refine proc~settree->proc~amr_coarsen_refine proc~timeintegration->proc~resettree program~agile agile program~agile->proc~main

Source Code

   subroutine alloc_buffers_f(self, nwgc, &
         ixS_p_min1, ixS_p_max1, &
         ixS_p_min2, ixS_p_max2, &
         ixS_p_min3, ixS_p_max3, &
         ixR_r_min1, ixR_r_max1, &
         ixR_r_min2, ixR_r_max2, &
         ixR_r_min3, ixR_r_max3)
     ! neighbor is fine
     ! required buffers:
     ! send: psb(ixS_p_) with indices type_send_p
     ! receive: psb(ixR_r_) with indices type_recv_r

     class(nbprocs_info_t)                     :: self
     integer, intent(in)                       :: nwgc
     integer, dimension(-1:1,0:3), intent(in)  :: &
         ixS_p_min1, ixS_p_max1, &
         ixS_p_min2, ixS_p_max2, &
         ixS_p_min3, ixS_p_max3, &
         ixR_r_min1, ixR_r_max1, &
         ixR_r_min2, ixR_r_max2, &
         ixR_r_min3, ixR_r_max3
     ! .. local ..
     integer         :: inb, isize_S, isize_R
     integer         :: igrid, inc1, inc2, inc3
     integer         :: iib1, iib2, iib3, ibuf_start
     integer         :: ixSmin1, ixSmin2, ixSmin3
     integer         :: ixRmin1, ixRmin2, ixRmin3
     integer         :: ixSmax1, ixSmax2, ixSmax3
     integer         :: ixRmax1, ixRmax2, ixRmax3

     self%imaxigrids_f = 0
     do inb = 1, self%nbprocs_f
        self%imaxigrids_f = max(self%fine_nb(inb)%info%nigrids, self%imaxigrids_f)

        isize_S = 0; isize_R = 0; ibuf_start = 1
        do igrid = 1, self%fine_nb(inb)%info%nigrids
           iib1 = idphyb(1,self%fine_nb(inb)%info%igrid(igrid))
           iib2 = idphyb(2,self%fine_nb(inb)%info%igrid(igrid))
           iib3 = idphyb(3,self%fine_nb(inb)%info%igrid(igrid))
           inc1 = self%fine_nb(inb)%info%inc1(igrid)
           inc2 = self%fine_nb(inb)%info%inc2(igrid)
           inc3 = self%fine_nb(inb)%info%inc3(igrid)

           ixSmin1=ixS_p_min1(iib1,inc1);ixSmin2=ixS_p_min2(iib2,inc2)
           ixSmin3=ixS_p_min3(iib3,inc3);ixSmax1=ixS_p_max1(iib1,inc1)
           ixSmax2=ixS_p_max2(iib2,inc2);ixSmax3=ixS_p_max3(iib3,inc3)

           ixRmin1=ixR_r_min1(iib1,inc1);ixRmin2=ixR_r_min2(iib2,inc2)
           ixRmin3=ixR_r_min3(iib3,inc3);ixRmax1=ixR_r_max1(iib1,inc1)
           ixRmax2=ixR_r_max2(iib2,inc2);ixRmax3=ixR_r_max3(iib3,inc3)

           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%fine_nb(inb)%info%ibuf_start(igrid) = ibuf_start
           self%fine_nb(inb)%info%isize(igrid) = (ixSmax1-ixSmin1+1) * (ixSmax2-ixSmin2+1) * (ixSmax3-ixSmin3+1) * nwgc

           ibuf_start = ibuf_start + self%fine_nb(inb)%info%isize(igrid)

        end do

        self%fine_nb(inb)%rcv%size       = isize_R*nwgc
        self%fine_nb(inb)%send%size      = isize_S*nwgc
        self%fine_nb(inb)%info_send%size = 5 * self%fine_nb(inb)%info%nigrids
        self%fine_nb(inb)%info_rcv%size  = 5 * self%fine_nb(inb)%info%nigrids

        allocate(self%fine_nb(inb)%rcv%buffer(isize_R * nwgc))
        allocate(self%fine_nb(inb)%send%buffer(isize_S * nwgc))
        allocate(self%fine_nb(inb)%info_rcv%buffer(5 * self%fine_nb(inb)%info%nigrids))
        allocate(self%fine_nb(inb)%info_send%buffer(5 * self%fine_nb(inb)%info%nigrids))

     end do

   end subroutine alloc_buffers_f