Lines Matching refs:p
48 integer :: i,j,k,m,n,p,q,ki,kj,ni,nj,qq
69 deallocate(ramg_flowDiag%p)
75 allocate(amg_paramap(1)%p(nshg))
76 allocate(amg_paraext(1)%p(nshg))
97 allocate(ramg_flowDiag%p(nshg,4))
109 ramg_flowDiag%p(:,i) = mflowDiag(:,i)
118 call ramg_initBCflag(amg_paramap(1)%p,ilwork,BC,iBC,iper)
127 amg_paraext(1)%p = amg_paramap(1)%p
129 if (amg_paramap(1)%p(i).ne.(myrank+1)) then
132 do m=lhsGPcolm%p(i),lhsGPcolm%p(i+1)-1
133 k = lhsGProwp%p(m)
134 if (amg_paraext(1)%p(k).eq.(myrank+1)) then
136 amg_paraext(1)%p(k)=amg_paramap(1)%p(i)
157 amg_A_colm(1)%p(i) = mnnz
159 do m = lhsGPcolm%p(i),lhsGPcolm%p(i+1)-1
160 k = lhsGProwp%p(m)
161 do j = lhsGPcolm%p(k),lhsGPcolm%p(k+1)-1
162 p = lhsGProwp%p(j)
165 if (temprow(p).ne.i) then
175 temprow(p) = i
182 amg_A_colm(1)%p(nshg+1) = mnnz
204 amg_A_rowp(1)%p(mnnz) = i
208 do m = lhsGPcolm%p(i),lhsGPcolm%p(i+1)-1
209 j = lhsGProwp%p(m)
210 do n = lhsGPcolm%p(j),lhsGPcolm%p(j+1)-1
211 k = lhsGProwp%p(n)
227 amg_A_rowp(1)%p(mnnz) = k
233 do m=amg_A_colm(1)%p(i)+1,amg_A_colm(1)%p(i+1)-1
234 ki = amg_A_rowp(1)%p(m)
235 do n=m+1,amg_A_colm(1)%p(i+1)-1
236 kj = amg_A_rowp(1)%p(n)
238 amg_A_rowp(1)%p(m) = kj
239 amg_A_rowp(1)%p(n) = ki
263 do m = amg_A_colm(1)%p(i),amg_A_colm(1)%p(i+1)-1
264 j = amg_A_rowp(1)%p(m)
269 ki = lhsGPcolm%p(i)
270 ni = lhsGPcolm%p(i+1)
271 kj = lhsGPcolm%p(j)
272 nj = lhsGPcolm%p(j+1)
278 & (lhsGProwp%p(ki).lt.lhsGProwp%p(kj)) )
283 & (lhsGProwp%p(kj).lt.lhsGProwp%p(ki)) )
287 p = lhsGProwp%p(ki)
288 q = lhsGProwp%p(kj)
289 if (p.eq.q) then
292 cki = amg_paramap(1)%p(i)
293 ckj = amg_paramap(1)%p(j)
294 ckk = amg_paramap(1)%p(k)
299 & + lhsGP%p(1,ki)*lhsGP%p(1,kj)*mflowDiag(k,1)**2
300 & + lhsGP%p(2,ki)*lhsGP%p(2,kj)*mflowDiag(k,2)**2
301 & + lhsGP%p(3,ki)*lhsGP%p(3,kj)*mflowDiag(k,3)**2
308 amg_A_lhs(1)%p(m,1)=rtemp
312 ki = amg_A_colm(1)%p(i)+1
313 mloop: do m=lhsGPcolm%p(i),lhsGPcolm%p(i+1)-1
314 j = lhsGProwp%p(m)
316 cki = amg_paramap(1)%p(i)
317 ckj = amg_paramap(1)%p(j)
320 amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1) =
321 & + amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1) + lhsGP%p(4,m)
325 do while(amg_A_rowp(1)%p(ki).lt.j)
328 cki = amg_paramap(1)%p(i)
329 ckj = amg_paramap(1)%p(j)
330 ckk = amg_paramap(1)%p(amg_A_rowp(1)%p(ki))
334 amg_A_lhs(1)%p(ki,1) = amg_A_lhs(1)%p(ki,1) + lhsGP%p(4,m)
342 deallocate(lhsGP%p)
343 deallocate(lhsGProwp%p)
344 deallocate(lhsGPcolm%p)
348 call ramg_global_lhs(amg_A_colm(1)%p,amg_A_rowp(1)%p,
349 & amg_A_lhs(1)%p,amg_nnz(1),
353 amg_A_colm(1)%p = lhsGPcolm%p
354 amg_nnz(1) = lhsGPcolm%p(amg_nshg(1)+1)-1
355 deallocate(amg_A_rowp(1)%p)
356 deallocate(amg_A_lhs(1)%p)
357 allocate(amg_A_rowp(1)%p(amg_nnz(1)))
358 allocate(amg_A_lhs(1)%p(amg_nnz(1),1))
359 amg_A_rowp(1)%p = lhsGProwp%p
360 amg_A_lhs(1)%p = lhsGP%p
361 deallocate(lhsGP%p)
362 deallocate(lhsGProwp%p)
363 deallocate(lhsGPcolm%p)
369 call ramg_dump_matlab_map(amg_A_colm(1)%p,amg_A_rowp(1)%p,
370 & amg_A_lhs(1)%p,
385 rtemp=amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1)
388 do j=amg_A_colm(1)%p(i)+1,amg_A_colm(1)%p(i+1)-1
389 rtemp = rtemp + amg_A_lhs(1)%p(j,1)
390 rt = amg_A_lhs(1)%p(j,1)
398 rtp = rtp/ (amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1))
399 rtn = rtn/amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1)
408 amg_ppeDiag(1)%p(i) =
409 & 1.0/sqrt(amg_A_lhs(1)%p(amg_A_colm(1)%p(i),1))
413 do m=amg_A_colm(1)%p(i),amg_A_colm(1)%p(i+1)-1
414 j = amg_A_rowp(1)%p(m)
415 amg_A_lhs(1)%p(m,1) = amg_A_lhs(1)%p(m,1) *
416 & amg_ppeDiag(1)%p(i)*amg_ppeDiag(1)%p(j)
418 amg_A_rhs(1)%p(i) = amg_A_rhs(1)%p(i)*amg_ppeDiag(1)%p(i)
424 amg_ppeDiag(1)%p(i) = amg_ppeDiag(1)%p(i)/mflowDiag(i,4)
429 call ramg_dump_matlab_map(amg_A_colm(1)%p,amg_A_rowp(1)%p,
430 & amg_A_lhs(1)%p,