| /phasta/phSolver/AMG/ |
| H A D | ramg_tools.f | 137 subroutine ramg_calcAv_g(level,u,v,colm,rowp,lhsK,lhsP, argument
152 integer,intent(in) :: level local
153 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: v
154 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u
162 IF (level.eq.1) THEN
167 call ramg_commOut(v,level,ilwork,1,iper,iBC,BC)
170 do i=1,amg_nshg(level)
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)
[all …]
|
| H A D | ramg_mls.f | 8 subroutine ramg_mls_eigen(evmax,level,flagfb, argument
18 integer,intent(in) :: level local
55 real(kind=8),dimension(amg_nshg(level)) :: twork1,twork2
64 asize = amg_nshg(level)
68 call ramg_generate_gmap(ilwork,asize,nsize,gmap,grevmap,level)
136 call ramg_calcAv_g(level,twork2,twork1,colm,rowp,lhsK,lhsP,
139 call ramg_mls_calcPAv(level,twork2,twork1,
156 ramg_levelx = level-1
188 subroutine ramg_mls_min_eigen(evmax,level,flagfb, argument
198 integer,intent(in) :: level local
[all …]
|
| H A D | ramg_cheby.f | 9 subroutine ramg_cheby_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument
26 real(kind=8),intent(in),dimension(amg_nshg(level)) :: v,r
27 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: u
31 real(kind=8),dimension(amg_nshg(level)) :: pAux
32 real(kind=8),dimension(amg_nshg(level)) :: epdk,epx
34 real(kind=8),dimension(amg_nshg(level)) :: cy1,cz1
41 beta = 1.1*mlsCf(level,1) ! 1.1 * Ev_max
42 alpha = mlsCf(level,1)*ramg_chebyratio
53 call ramg_calcAv_g(level,pAux,v,colm,rowp,lhsK,lhsP,
58 call ramg_calcAv_g(level,pAux,cy1,colm,rowp,lhsK,lhsP,
[all …]
|
| H A D | ramg_smoother.f | 395 subroutine ramg_smoother(level,xnew,xold,xrhs, argument
403 integer level local
404 real(kind=8),dimension(amg_nshg(level)),intent(in) ::xrhs,xold
405 real(kind=8),dimension(amg_nshg(level)),intent(inout) :: xnew
429 call ramg_mls_apply(xnew,xold,xrhs,level,
433 call ramg_gauss(amg_A_colm(level)%p,amg_A_rowp(level)%p,
434 & amg_A_lhs(level)%p,amg_nshg(level),
435 & amg_nnz(level),xrhs,xnew,xold,
436 & fwdbck,level,
439 call ramg_cheby_apply(xnew,xold,xrhs,level,
[all …]
|
| H A D | ramg_paratools.f | 893 subroutine ramg_commOut(global,level,ilwork,redun,iper,iBC,BC) argument
899 integer,intent(in) :: level local
901 real(kind=8),intent(inout),dimension(amg_nshg(level),redun)
916 if (level.eq.1) then
921 numtask = amg_ilwork(level)%p(1)
931 iother(i)=amg_ilwork(level)%p(itkbeg+1)
932 numseg(i)=amg_ilwork(level)%p(itkbeg+2)
937 tmparray(i)%p(j,:)=global(amg_ilwork(level)%p(itkbeg),:)
941 if ((myrank.eq.1).and.(level.eq.2)) then ! debug
975 numseg(i)=amg_ilwork(level)%p(itkbeg+2)
[all …]
|
| H A D | ramg_driver.f | 106 & sol,rhs,level, argument
113 integer,intent(in) :: level local
114 real(kind=8),intent(in),dimension(amg_nshg(level)) :: rhs
115 real(kind=8),intent(inout),dimension(amg_nshg(level)) :: sol
129 real(kind=8),dimension(amg_nshg(level)) :: cgP,cgQ,cgZ,cgR
135 call ramg_L2_norm_p(level,cgR,1,norm_0)
141 call ramg_dot_p(level,cgZ,cgR,1,rz)
149 call ramg_calcAv_g(level,cgQ,cgP,colm,rowp,lhsK,lhsP,
151 call ramg_dot_p(level,cgP,cgQ,1,pq)
155 call ramg_L2_norm_p(level,cgR,1,tmp)
[all …]
|
| H A D | amgread.txt | 85 information is computed at very beginning for the finest level. Then
86 after each coarsening, the information is carried on to higher level in
109 Also in ramg_control, we do lu decomposition for the coarsest level, if
110 direct solve for coarsest level is selected. Thus, each iteration will
112 especially when you have a bit larger coarest level.
128 Pay attention to diagonal scaling in V-cycle, at the finest level, the
130 second level and above.
133 improve efficiency: For second level and above, scaling can be
136 coarse level vector. We don't need to do this every iteration. This can
230 init_ilwork initialize all the ilwork structure for higher level AMG
[all …]
|
| H A D | ramg_ggb.f | 194 & level) argument
205 integer,intent(in) :: level local
249 m = iabs(amg_paramap(level)%p(i))
|
| /phasta/phSolver/incompressible/e3source/ |
| H A D | e3source.fimplicit | 493 elseif (isclr.eq.2) then !we are redistancing level-sets
558 endif ! close of scalar 2 of level set
560 else ! NOT turbulence and NOT level set so this is a simple
|
| H A D | e3source.fRHSonly | 479 elseif (isclr.eq.2) then !we are redistancing level-sets
544 endif ! close of scalar 2 of level set
546 else ! NOT turbulence and NOT level set so this is a simple
|
| /phasta/phSolver/ |
| H A D | CHANGES | 367 …" throughout the code with iprec and ilhs flags being set at the itrdrv.f level with the exception…
|