36       integer                           :: i,j,k,p  local
44          do k=I_cf_colm(level1)%p(i),I_cf_colm(level1)%p(i+1)-1
45             j = I_cf_rowp(level1)%p(k)
46             vf(i) = vf(i) + VC(j)*I_cf(level1)%p(k)
72       integer                           :: i,j,k,p  local
80          do k=I_fc_colm(level1)%p(i),I_fc_colm(level1)%p(i+1)-1
81             j = I_fc_rowp(level1)%p(k)
82             VC(i) = VC(i) + vf(j)*I_fc(level1)%p(k)
117       integer                             :: i,j,k,p  local
157       integer                             :: i,j,k,p,mi,mj  local
171          mi = amg_paramap(level)%p(i)
172          do k=amg_A_colm(level)%p(i),amg_A_colm(level)%p(i+1)-1
173             j = amg_A_rowp(level)%p(k)
174             mj = amg_paramap(level)%p(j)
176                 u(i) = u(i)+amg_A_lhs(level)%p(k,1)*v(j)
217       integer                             :: i,j,k,p,mi,mj,mk  variable
227          mi = amg_paramap(level)%p(i)
228          mk = amg_paraext(level)%p(i)
230          do k=amg_A_colm(level)%p(i)+diag,amg_A_colm(level)%p(i+1)-1
231             j = amg_A_rowp(level)%p(k)
232             mj = amg_paramap(level)%p(j)
234                 u(i) = u(i)+amg_A_lhs(level)%p(k,1)*v(j)
266          allocate(tarray(i)%p(amg_nshg(i)))
268             call random_number(tarray(i)%p)
274          call ramg_dump(tarray(i)%p,amg_nshg(i),fname)
278          call ramg_commIn(tarray(i)%p,i,ilwork,1,iper,iBC,BC)
280          call ramg_dump(tarray(i)%p,amg_nshg(i),fname)
284          call ramg_commOut(tarray(i)%p,i,ilwork,1,iper,iBC,BC)
286          call ramg_dump(tarray(i)%p,amg_nshg(i),fname)
298          deallocate(tarray(i)%p)
316       integer i,j,k,p  local
353       integer      :: i,p  local
372                 p = ilwork(itkbeg+3+j*2)
373                 amg_F(p) = 2
449       integer                              :: i,j,p  local
456          p = colm(i)
457          if (rowp(p).ne.i) then
461          diagrow = lhs(p)
467             p = rowp(j)
468             if (lhs(p).gt.diagrow) then
471                 write(*,*)'diag:',diagrow,p,lhs(p)
566       allocate(amg_A_lhs(level)%p(amg_nnz(level),n_sol),
569       allocate(amg_A_rowp(level)%p(amg_nnz(level)),
575       allocate(amg_A_colm(level)%p(amg_nshg(level)+1),
578       allocate(amg_A_rhs(level)%p(amg_nshg(level)),
581       allocate(amg_ppeDiag(level)%p(amg_nshg(level)),
591       amg_A_lhs(level)%p(:,:)  = 0
592       amg_A_rowp(level)%p(:) = 0
595       amg_A_colm(level)%p(:) = 0
596       amg_A_rhs(level)%p(:)  = 0
632       if (associated(amg_A_lhs(level)%p)) then
633       deallocate(amg_A_lhs(level)%p,
637       if (associated(amg_A_rowp(level)%p)) then
638       deallocate(amg_A_rowp(level)%p,
642       if (associated(amg_A_colm(level)%p)) then
643       deallocate(amg_A_colm(level)%p,
647       if (associated(amg_ppeDiag(level)%p)) then
648       deallocate(amg_ppeDiag(level)%p,
652       if (associated(amg_A_rhs(level)%p)) then
653       deallocate(amg_A_rhs(level)%p,
657       if (associated(amg_ilwork(level)%p)) then
658       deallocate(amg_ilwork(level)%p,
664       if (associated(amg_paramap(level)%p)) then
665       deallocate(amg_paramap(level)%p,
669       if (associated(amg_paraext(level)%p)) then
670       deallocate(amg_paraext(level)%p,
679       if (associated(I_fc_colm(level)%p)) then
680       deallocate(CF_map(level)%p,stat=mem_err)
682       deallocate(CF_revmap(level)%p,stat=mem_err)
684       deallocate(I_cf_colm(level)%p,stat=mem_err)
686       deallocate(I_cf_rowp(level)%p,stat=mem_err)
688       deallocate(I_cf(level)%p,stat=mem_err)
690       deallocate(I_fc_colm(level)%p,stat=mem_err)
692       deallocate(I_fc_rowp(level)%p,stat=mem_err)
694       deallocate(I_fc(level)%p,stat=mem_err)
903       integer i,j,p,k  local
914          do p=acolm(i),acolm(i+1)-1
915             j = arowp(p)
920             write(264,mformat)i,j,(alhs(k,p),k=1,redun)
945       integer i,j,p,k  local
965          do p=acolm(i),acolm(i+1)-1
966             j = arowp(p)
969             write(264,mformat)ii,jj,(alhs(k,p),k=1,redun)
986       integer i,j,p,k  local
995          do p=acolm(i),acolm(i+1)-1
996             j = arowp(p)
997             write(264,*)mAname,'[',i,',',j,']:=',alhs(p),':'