| /phasta/phSolver/AMG/ |
| H A D | ramg_driver.f | 17 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper argument 32 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 52 & colm,rowp,lhsK,lhsP, 64 & colm,rowp,lhsK,lhsP, 68 call ramg_calcAv_g(clevel,myvF,ramg_sol,colm,rowp,lhsK,lhsP, 75 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper) 83 & colm,rowp,lhsK,lhsP, 107 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper,iterNum) argument 121 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 149 call ramg_calcAv_g(level,cgQ,cgP,colm,rowp,lhsK,lhsP, [all …]
|
| H A D | ramg_mls.f | 9 & colm,rowp,lhsK,lhsP, argument 24 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 136 call ramg_calcAv_g(level,twork2,twork1,colm,rowp,lhsK,lhsP, 140 & colm,rowp,lhsK,lhsP, 189 & colm,rowp,lhsK,lhsP, argument 204 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 316 call ramg_calcAv_g(level,twork2,twork1,colm,rowp,lhsK,lhsP, 320 & colm,rowp,lhsK,lhsP, 371 subroutine ramg_mls_setup(colm,rowp,lhsK,lhsP, argument 381 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local [all …]
|
| H A D | ramg_cheby.f | 9 subroutine ramg_cheby_apply(u,v,r,level,colm,rowp,lhsK,lhsP, argument 23 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 53 call ramg_calcAv_g(level,pAux,v,colm,rowp,lhsK,lhsP, 58 call ramg_calcAv_g(level,pAux,cy1,colm,rowp,lhsK,lhsP, 82 call ramg_calcAv_g(level,pAux,v,colm,rowp,lhsK,lhsP, 90 call ramg_calcAv_g(level,pAux,epx,colm,rowp,lhsK,lhsP, 106 subroutine ramg_cheby_setup(colm,rowp,lhsK,lhsP, argument 116 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 129 call ramg_mls_eigen(mlsCf(lvl,1),lvl,1,colm,rowp,lhsK,lhsP,
|
| H A D | ramg_smoother.f | 5 subroutine ramg_calcPPEv(colm,rowp,lhsK,lhsP, argument 15 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 28 call fLesSparseApG(colm,rowp,lhsP,v,utmp,nshg,nnz_tot) 42 call fLesSparseApNGtC(colm,rowp,lhsP,utmp4,u,nshg,nnz_tot) 325 & colm,rowp,lhsK,lhsP, argument 338 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 362 & colm,rowp,lhsK,lhsP, 365 & colm,rowp,lhsK,lhsP, 372 & colm,rowp,lhsK,lhsP, 375 & colm,rowp,lhsK,lhsP, [all …]
|
| H A D | ramg_ggb.f | 1 subroutine ramg_ggb_setup(colm,rowp,lhsK,lhsP, argument 14 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 112 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper) 152 & colm,rowp,lhsK,lhsP, 270 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper argument 287 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 300 call ramg_calcAv_g(clevel,vF,ramg_rhs,colm,rowp,lhsK,lhsP, 304 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper) 316 & colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper argument 330 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local [all …]
|
| H A D | ramg_control.f | 7 subroutine ramg_control(colm,rowp,lhsK,lhsP, argument 15 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 35 call ramg_extract_ppe(colm,rowp,lhsK,lhsP, 42 call ramg_cheby_setup(colm,rowp,lhsK,lhsP, 45 call ramg_mls_setup(colm,rowp,lhsK,lhsP, 57 call ramg_ggb_setup(colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper)
|
| H A D | ramg_extract.f | 13 & colm,rowp,lhsK,lhsP, 31 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP 87 & rowp,colm,lhsK,lhsP) 120 call ramg_global_lhs(colm,rowp,lhsP,nnz_tot,
|
| H A D | ramg_paratools.f | 487 & colm,rowp,lhsK,lhsP, argument 498 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 516 call fLesSparseApG(colm,rowp,lhsP,t1,t2,nshg,nnz_tot) 527 call fLesSparseApNGtC(colm,rowp,lhsP,t2b,t1,nshg,nnz_tot) 543 & colm,rowp,lhsK,lhsP, argument 555 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 570 call fLesSparseApG(colm,rowp,lhsP,t1,t2,nshg,nnz_tot) 581 call fLesSparseApNGtC(colm,rowp,lhsP,t2b,t1,nshg,nnz_tot) 594 & colm,rowp,lhsK,lhsP, argument 606 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local [all …]
|
| H A D | ramg_tools.f | 137 subroutine ramg_calcAv_g(level,u,v,colm,rowp,lhsK,lhsP, argument 151 real(kind=8),intent(in),dimension(4,nnz_tot) :: lhsP local 164 call ramg_PPEAp(u,v,colm,rowp,lhsK,lhsP,ilwork,BC,iBC,iper)
|
| H A D | amgread.txt | 152 Now to parallel. lhsGP* is a sparse-communicated duplicate of lhsP 210 global_lhs is a function that reads in lhsP (in sparse format so there 215 alter the structure of lhsP too. This is done task by task. There are
|
| /phasta/phSolver/common/ |
| H A D | fillsparse.f | 2 & xGoC, lhsP, argument 10 real*8 lhsK(9,nnz_tot), lhsP(4,nnz_tot) local 54 lhsP(1,k) = lhsP(1,k) + xGoC(e,1,aa,b) 55 lhsP(2,k) = lhsP(2,k) + xGoC(e,2,aa,b) 56 lhsP(3,k) = lhsP(3,k) + xGoC(e,3,aa,b) 57 lhsP(4,k) = lhsP(4,k) + xGoC(e,4,aa,b)
|
| H A D | fillsparse.c | 21 void fillsparsecpetscs(gcorp_t* ieng, double* EGmass, Mat* lhsP) in fillsparsecpetscs() argument 40 MatSetValues(*lhsP, petsc_nshl , locat, petsc_nshl, locat, mb, ADD_VALUES); in fillsparsecpetscs() 45 void fillsparsecpetscc(gcorp_t* ieng, double* EGmass, Mat* lhsP) in fillsparsecpetscc() argument 66 MatSetValuesBlocked(*lhsP, petsc_nshl , locat, petsc_nshl, locat, mb, ADD_VALUES); in fillsparsecpetscc()
|
| /phasta/phSolver/incompressible/ |
| H A D | solvedata.f | 5 real*8, allocatable :: lhsP(:,:), lhsK(:,:), lhsS(:,:) variable 17 allocate (lhsP(4,nnz_tot)) 42 deallocate (lhsP)
|
| H A D | elmStats.f | 82 & lhsK, lhsP ) argument 93 & lhsP(4,nnz_tot), res(nshg,ndof) local 112 & lhsP, rjunk )
|
| H A D | solfar.f | 150 & lhsP, rerr, GradV ) 181 Val4(13:15,i) = lhsP(1:3,i) 182 Val4(16,i) = lhsP(4,i) 191 Val4(4:12:4,l) = -lhsP(1:3,j) 225 call ramg_control(colm,rowp,lhsK,lhsP, 252 & lhsK, lhsP, rdtmp, 278 & lhsK, lhsP, rdtmp,
|
| H A D | lestools.c | 30 int *colm, double *lhsK, double *lhsP) ; 91 double *lhsP, double *lesP, double *lesQ, 106 double *lhsP, double *lesP, double *lesQ, 110 double *lhsP, double *lesP, double *lesQ, 114 double *lhsP, double *lesP, double *lesQ, 118 double *lhsP, double *lesP, double *lesQ, 120 void RAMG_INTERFACE ( int *colm, int *rowp, double *lhsK,double *lhsP, 142 usrHd->lhsK, usrHd->lhsP) ; in lesPrepDiag() 914 fLesSparseApG( usrHd->colm, usrHd->rowp, usrHd->lhsP, in lesApG() 1000 usrHd->lhsP, usrHd->lesP, usrHd->lesQ, in lesApKG() [all …]
|
| H A D | usr.h | 41 double* lhsP ; /* sparse GoC matrix (4,nnzeros) */ member 131 double* lhsP,
|
| H A D | lesSparse.f | 16 & lhsK, lhsP) 28 real*8 lhsK(9,nnz_tot), lhsP(4,nnz_tot) 46 flowdiag(n,4) = lhsP(4,k) 61 flowdiag(n,4) = flowdiag(n,4) + abs(lhsP(4,k))
|
| H A D | usr.c | 62 double* lhsP, in usrNew() argument 95 usrHd->lhsP = lhsP ; in usrNew()
|
| H A D | elmgmr.f | 6 & lhsP, rerr, GradV) argument 48 real*8 lhsK(9,nnz_tot), lhsP(4,nnz_tot) local 185 & xGoC, lhsP, 311 & xGoC, lhsP,
|
| H A D | soldir.f | 56 & lhsP, rerr )
|
| H A D | stats.f | 191 & rowp, colm, lhsK, lhsP ) argument 204 & lhsP(4,nnz_tot), x(numnp,nsd) local 236 & lhsK, lhsP )
|
| H A D | itrdrv.f | 501 & rowp, colm, lhsK, lhsP )
|
| /phasta/phSolver/compressible/ |
| H A D | solgmrpetsc.c | 33 static Mat lhsP; variable 203 0, idiagnz, 0, iodiagnz, &lhsP); in SolGMRp() 205 ierr = MatSetOption(lhsP, MAT_NEW_NONZERO_ALLOCATION_ERR, PETSC_FALSE); in SolGMRp() 209 ierr = MatSetOption(lhsP, MAT_SUBSET_OFF_PROC_ENTRIES, PETSC_TRUE); in SolGMRp() 211 ierr = MatSetUp(lhsP); in SolGMRp() 214 ierr = MatGetOwnershipRange(lhsP, &myMatStart, &myMatEnd); in SolGMRp() 218 if(genpar.lhs ==1) ierr = MatZeroEntries(lhsP); in SolGMRp() 238 ilwork, rerr , &lhsP); in SolGMRp() 293 ierr = MatAssemblyBegin(lhsP, MAT_FINAL_ASSEMBLY); in SolGMRp() 294 ierr = MatAssemblyEnd(lhsP, MAT_FINAL_ASSEMBLY); in SolGMRp() [all …]
|
| H A D | elmgmrpetsc.f | 11 & rerr, lhsP) argument 195 call fillsparsecpetscc( mieng(iblk)%p, EGmass, lhsP)
|