| <- HREF="index.html" PREV | Index | Next -> |
NHSE Review^TM: Comments · Archive · Search

      program two_triples
      common /pghpf_0/pghpf_0, /pghpf_0c/pghpf_0c
     +, /pghpf_lineno/pghpf_lineno, /pghpf_0l/pghpf_0l
      integer pghpf_0(4), pghpf_lineno, pghpf_0l(8)
      character*1 pghpf_0c
      common /pghpf_local_mode / pghpf_local_mode
      integer pghpf_local_mode
      integer i$i, i$i1, i$i2, i$i3, i$i4, i$i5, i$i6, i$$l, i$$u, i$$l1
     +, i$$u1, i$$l2, i$$u2, i$$l3, i$$u3, a1$tl$d(1), a1$d(1), b2$d(1)
     +, a1$d1(1), b2$d1(1), a1$d2(1), l$$b, u$$b, l$$b1, u$$b1, l$$b2, 
     +u$$b2, l$$b3, u$$b3
      real a1(1), a2(1), b1(1), b2(1)
      integer cp$$com, xfer$$com, cp$$com1, xfer$$com1
      integer a1$p, a2$p, b1$p, b2$p, a1$tl$dp, a1$dp, b2$dp, a1$dp1, 
     +b2$dp1, a1$dp2
      integer pghpf_olap_shift, pghpf_comm_start
      data pghpf_0l /852078511,-1726180828,927164202,-291367439,1,
     +1261038733,0,1627158493/
      call pghpf_init(0)
      pghpf_local_mode = 0
      call pghpf_template(a1$dp1,2,1024,0,0,0,1,30,1,40)
      call pghpf_ptr_offset(a1$dp1,a1$dp1,a1$d1,25)
      call pghpf_instance(a1$dp2,a1$d1(a1$dp1),27,4,0,4,4,4,4)
      call pghpf_ptr_offset(a1$dp2,a1$dp2,a1$d2,25)
      call pghpf_template(b2$dp,2,20480,a1$d1(a1$dp1),0,1,30,1,40)
      call pghpf_ptr_offset(b2$dp,b2$dp,b2$d,25)
      call pghpf_instance(b2$dp1,b2$d(b2$dp),27,4,0,4,4,4,4)
      call pghpf_ptr_offset(b2$dp1,b2$dp1,b2$d1,25)
      call pghpf_allobnds(a1$d2(a1$dp2),l$$b,u$$b,l$$b1,u$$b1)
      call pgf90_alloc((u$$b-l$$b+1)*(u$$b1-l$$b1+1),27,pghpf_0(3),
     +pghpf_0(3),a1$p,a1)
      call pghpf_allobnds(b2$d1(b2$dp1),l$$b2,u$$b2,l$$b3,u$$b3)
      call pgf90_alloc((u$$b2-l$$b2+1)*(u$$b3-l$$b3+1),27,pghpf_0(3),
     +pghpf_0(3),a2$p,a2)
      call pgf90_alloc((u$$b2-l$$b2+1)*(u$$b3-l$$b3+1),27,pghpf_0(3),
     +pghpf_0(3),b1$p,b1)
      call pgf90_alloc((u$$b2-l$$b2+1)*(u$$b3-l$$b3+1),27,pghpf_0(3),
     +pghpf_0(3),b2$p,b2)
      call init(a1(a1$p),b1(b1$p),a1$d2(a1$dp2),b2$d1(b2$dp1))


      cp$$com = pghpf_olap_shift(a1(a1$p),a1$d2(a1$dp2),1,1,1,1)
      xfer$$com = pghpf_comm_start(cp$$com,a1(a1$p),a1$d2(a1$dp2),a1(
     +a1$p),a1$d2(a1$dp2))
      call pghpf_comm_finish(xfer$$com)
      call pghpf_localize_bounds(b2$d1(b2$dp1),1,1,29,1,i$$l,i$$u)
      call pghpf_localize_bounds(b2$d1(b2$dp1),2,1,39,1,i$$l1,i$$u1)
!     forall (i$i=i$$l1:i$$u1:1, i$i1=i$$l:i$$u:1) a2((u$$b2-l$$b2+1)*(
!    +i$i-l$$b3)+i$i1-l$$b2+a2$p) = a1(i$i1+1+(u$$b-l$$b+1)*(i$i+1-l$$b1
!    +)-l$$b+a1$p)
      do i$i = i$$l1, i$$u1
         do i$i1 = i$$l, i$$u
            a2((u$$b2-l$$b2+1)*(i$i-l$$b3)+i$i1-l$$b2+a2$p) = a1(i$i1+1+
     +(u$$b-l$$b+1)*(i$i+1-l$$b1)-l$$b+a1$p)
         enddo
      enddo
      xfer$$com1 = pghpf_comm_start(cp$$com,b1(b1$p),b2$d1(b2$dp1),b1(
     +b1$p),b2$d1(b2$dp1))
      call pghpf_comm_finish(xfer$$com1)
      call pghpf_localize_bounds(b2$d1(b2$dp1),1,2,30,1,i$$l2,i$$u2)
      call pghpf_localize_bounds(b2$d1(b2$dp1),2,2,40,1,i$$l3,i$$u3)
!     forall (i$i=i$$l3:i$$u3:1, i$i1=i$$l2:i$$u2:1) b2((u$$b2-l$$b2+1)*
!    +(i$i-l$$b3)+i$i1-l$$b2+b2$p) = b1(i$i1-1+(u$$b2-l$$b2+1)*(i$i-1-
!    +l$$b3)-l$$b2+b1$p)
      do i$i = i$$l3, i$$u3
         do i$i1 = i$$l2, i$$u2
            b2((u$$b2-l$$b2+1)*(i$i-l$$b3)+i$i1-l$$b2+b2$p) = b1(i$i1-1+
     +(u$$b2-l$$b2+1)*(i$i-1-l$$b3)-l$$b2+b1$p)
         enddo
      enddo
      call pghpf_comm_free(1,cp$$com)

      call use(a2(a2$p),b2(b2$p),b2$d1(b2$dp1),b2$d1(b2$dp1))
99999 continue
      call pgf90_deallocate(a1(a1$p),pghpf_0(3))
      call pgf90_deallocate(a2(a2$p),pghpf_0(3))
      call pgf90_deallocate(b1(b1$p),pghpf_0(3))
      call pgf90_deallocate(b2(b2$p),pghpf_0(3))
      call pghpf_free(b2$d1(b2$dp1))
      call pghpf_free(b2$d(b2$dp))
      call pghpf_free(a1$d2(a1$dp2))
      call pghpf_free(a1$d1(a1$dp1))
      call pgf90_exit(0)
      end

| <- HREF="index.html" PREV | Index | Next -> |
NHSE Review^TM: Comments · Archive · Search