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

      program sum_prog
      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, a1$tl$tl$d(1), a1$a$d(1), a1$a$d1(1), l$$b, u$$b, l$$b1, 
     +u$$b1, a1$tl$tl$d1(1), a1$d(1), a2$d(1), a1$d1(1), a2$d1(1), a1$d2
     +(1), l$$b2, u$$b2, l$$b3, u$$b3, l$$b4, u$$b4, l$$b5, u$$b5
      real a1(1), a2(1), a1$a(1), a1$a$r
      integer a1$p, a2$p, a1$a$p, a1$tl$tl$dp, a1$a$dp, a1$a$dp1, 
     +a1$tl$tl$dp1, a1$dp, a2$dp, a1$dp1, a2$dp1, a1$dp2
      data pghpf_0l /852079647,-1647718444,927163098,-265212932,1,
     +1261038733,0,1680421997/
      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,0,0,0,0)
      call pghpf_ptr_offset(a1$dp2,a1$dp2,a1$d2,25)
      call pghpf_template(a2$dp,2,20480,a1$d1(a1$dp1),0,1,30,1,40)
      call pghpf_ptr_offset(a2$dp,a2$dp,a2$d,25)
      call pghpf_instance(a2$dp1,a2$d(a2$dp),27,4,0,0,0,0,0)
      call pghpf_ptr_offset(a2$dp1,a2$dp1,a2$d1,25)
      call pghpf_allobnds(a1$d2(a1$dp2),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),a1$p,a1)
      call pghpf_allobnds(a2$d1(a2$dp1),l$$b4,u$$b4,l$$b5,u$$b5)
      call pgf90_alloc((u$$b4-l$$b4+1)*(u$$b5-l$$b5+1),27,pghpf_0(3),
     +pghpf_0(3),a2$p,a2)
      call pghpf_template(a1$tl$tl$dp,2,1024,0,0,0,1,30,1,40)
      call pghpf_ptr_offset(a1$tl$tl$dp,a1$tl$tl$dp,a1$tl$tl$d,25)

      call pghpf_template(a1$a$dp,2,20480,a1$tl$tl$d(a1$tl$tl$dp),0,1,30
     +,1,40)
      call pghpf_ptr_offset(a1$a$dp,a1$a$dp,a1$a$d,25)
      call pghpf_instance(a1$a$dp1,a1$a$d(a1$a$dp),27,4,0,0,0,0,0)
      call pghpf_ptr_offset(a1$a$dp1,a1$a$dp1,a1$a$d1,25)
      call pghpf_allobnds(a1$a$d1(a1$a$dp1),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$a$p,a1$a)

      call pghpf_localize_bounds(a1$a$d1(a1$a$dp1),1,1,30,1,i$$l,i$$u)
      call pghpf_localize_bounds(a1$a$d1(a1$a$dp1),2,1,40,1,i$$l1,i$$u1)
!     forall (i$i=i$$l1:i$$u1:1, i$i1=i$$l:i$$u:1) a1$a((u$$b-l$$b+1)*(
!    +i$i-l$$b1)+i$i1-l$$b+a1$a$p) = a1((u$$b2-l$$b2+1)*(i$i-l$$b3)+i$i1
!    +-l$$b2+a1$p) + a2((u$$b4-l$$b4+1)*(i$i-l$$b5)+i$i1-l$$b4+a2$p)
      do i$i = i$$l1, i$$u1
         do i$i1 = i$$l, i$$u
            a1$a((u$$b-l$$b+1)*(i$i-l$$b1)+i$i1-l$$b+a1$a$p) = a1((u$$b2
     +-l$$b2+1)*(i$i-l$$b3)+i$i1-l$$b2+a1$p) + a2((u$$b4-l$$b4+1)*(i$i-
     +l$$b5)+i$i1-l$$b4+a2$p)
         enddo
      enddo
      call pghpf_sums(a1$a$r,a1$a(a1$a$p),.true.,27,a1$a$d1(a1$a$dp1),19
     +)
      s = a1$a$r

      call pghpf_freen(3,a1$a$d1(a1$a$dp1),a1$a$d(a1$a$dp),a1$tl$tl$d(
     +a1$tl$tl$dp))
      call pgf90_deallocate(a1$a(a1$a$p),pghpf_0(3))
      call use_it(s,27)
99999 continue
      call pgf90_deallocate(a1(a1$p),pghpf_0(3))
      call pgf90_deallocate(a2(a2$p),pghpf_0(3))
      call pghpf_free(a2$d1(a2$dp1))
      call pghpf_free(a2$d(a2$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