/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/interface/ |
H A D | itfunc.c | 330 if (ksp->dmActive && !ksp->setupstage) { in KSPSetUp() 375 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPSetUp() 382 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPSetUp() 633 …if (ksp->viewPre) {ierr = ObjectView((PetscObject) ksp, ksp->viewerPre, ksp->formatPre);CHKERRQ(ie… in KSPSolve_Private() 635 …if (ksp->presolve) {ierr = (*ksp->presolve)(ksp,ksp->vec_rhs,ksp->vec_sol,ksp->prectx);CHKERRQ(ier… in KSPSolve_Private() 730 ksp->totalits += ksp->its; in KSPSolve_Private() 757 ierr = KSPGuessUpdate(ksp->guess,ksp->vec_rhs,ksp->vec_sol);CHKERRQ(ierr); in KSPSolve_Private() 760 ierr = (*ksp->postsolve)(ksp,ksp->vec_rhs,ksp->vec_sol,ksp->postctx);CHKERRQ(ierr); in KSPSolve_Private() 766 …if (ksp->viewSV) {ierr = KSPViewSingularvalues_Internal(ksp, ksp->viewerSV, ksp->formatSV);C… in KSPSolve_Private() 1286 ksp->pc_side = ksp->pc_side_set = side; in KSPSetPCSide() [all …]
|
H A D | itcl.c | 62 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPSetOptionsPrefix() 92 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPAppendOptionsPrefix() 160 ksp->guess->ksp = ksp; in KSPSetGuess() 194 ksp->guess->ksp = ksp; in KSPGetGuess() 333 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPSetFromOptions() 358 …= PetscOptionsGetViewer(comm,((PetscObject) ksp)->options,prefix,"-ksp_view",&ksp->viewer, &ksp->f… in KSPSetFromOptions() 360 …tscOptionsGetViewer(comm,((PetscObject) ksp)->options,prefix,"-ksp_view_mat",&ksp->viewerMat,&ksp-… in KSPSetFromOptions() 362 …tscOptionsGetViewer(comm,((PetscObject) ksp)->options,prefix,"-ksp_view_rhs",&ksp->viewerRhs,&ksp-… in KSPSetFromOptions() 405 ierr = KSPGetGuess(ksp,&ksp->guess);CHKERRQ(ierr); in KSPSetFromOptions() 544 … = PetscOptionsGetViewer(comm,((PetscObject) ksp)->options,prefix,"-ksp_view",&ksp->viewer,&ksp->f… in KSPSetFromOptions() [all …]
|
H A D | itcreate.c | 128 ierr = (*ksp->ops->view)(ksp,viewer);CHKERRQ(ierr); in KSPView() 172 if (ksp->ops->view) {ierr = (*ksp->ops->view)(ksp,viewer);CHKERRQ(ierr);} in KSPView() 212 ierr = (*ksp->ops->view)(ksp,viewer);CHKERRQ(ierr); in KSPView() 288 ksp->normtype = ksp->normtype_set = normtype; in KSPSetNormType() 391 ksp->pc_side = ksp->pc_side_set; in KSPNormSupportTableReset_Private() 392 ksp->normtype = ksp->normtype_set; in KSPNormSupportTableReset_Private() 404 …if ((ksp->normtype == KSP_NORM_DEFAULT || ksp->normtype == i) && (ksp->pc_side == PC_SIDE_DEFAULT … in KSPSetUpNorms_Private() 528 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPSetOperators() 560 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPGetOperators() 588 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPGetOperatorsSet() [all …]
|
H A D | iterativ.c | 32 *rnorm = ksp->rnorm; in KSPGetResidualNorm() 60 *its = ksp->its; in KSPGetIterationNumber() 744 ierr = KSP_PCApply(ksp,ksp->vec_rhs,z);CHKERRQ(ierr); in KSPConvergedDefault() 766 ksp->ttol = PetscMax(ksp->rtol*ksp->rnorm0,ksp->abstol); in KSPConvergedDefault() 797 } else if (rnorm >= ksp->divtol*ksp->rnorm0) { in KSPConvergedDefault() 852 if (ksp->pc) { in KSPBuildSolutionDefault() 906 if (!ksp->pc) {ierr = KSPGetPC(ksp,&ksp->pc);CHKERRQ(ierr);} in KSPBuildResidualDefault() 952 if (ksp->vec_sol) vecr = ksp->vec_sol; in KSPCreateVecs() 982 if (ksp->vec_rhs) vecl = ksp->vec_rhs; in KSPCreateVecs() 1203 if (!ksp->dm) { in KSPGetDM() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cg/ |
H A D | cg.c | 123 ksp->its = 0; in KSPSolve_CG() 154 ksp->rnorm = dp; in KSPSolve_CG() 156 …ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for convergen… in KSPSolve_CG() 169 ksp->its = i+1; in KSPSolve_CG() 228 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CG() 241 if (i >= ksp->max_it) ksp->reason = KSP_DIVERGED_ITS; in KSPSolve_CG() 284 ksp->its = 0; in KSPSolve_CG_SingleReduction() 317 ksp->rnorm = dp; in KSPSolve_CG_SingleReduction() 319 …ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for convergen… in KSPSolve_CG_SingleReduction() 402 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CG_SingleReduction() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/rich/ |
H A D | rich.c | 38 ksp->its = 0; in KSPSolve_Richardson() 52 r = ksp->work[0]; in KSPSolve_Richardson() 53 z = ksp->work[1]; in KSPSolve_Richardson() 65 …ierr = PCApplyRichardson(ksp->pc,b,x,r,ksp->rtol,ksp->abstol,ksp->divtol,maxit,ksp->guess_z… in KSPSolve_Richardson() 79 ksp->its = 0; in KSPSolve_Richardson() 94 ierr = (*ksp->converged)(ksp,i,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Richardson() 103 ksp->its++; in KSPSolve_Richardson() 117 ierr = (*ksp->converged)(ksp,i,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Richardson() 144 if (ksp->its >= ksp->max_it) { in KSPSolve_Richardson() 146 ierr = (*ksp->converged)(ksp,i,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Richardson() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/bcgs/ |
H A D | bcgs.c | 34 X = ksp->vec_sol; in KSPSolve_BCGS() 35 B = ksp->vec_rhs; in KSPSolve_BCGS() 47 if (ksp->pc_side == PC_RIGHT && !ksp->guess_zero) { in KSPSolve_BCGS() 61 ksp->its = 0; in KSPSolve_BCGS() 62 ksp->rnorm = dp; in KSPSolve_BCGS() 66 ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_BCGS() 123 if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i+2) { in KSPSolve_BCGS() 132 ksp->its++; in KSPSolve_BCGS() 137 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_BCGS() 146 if (i >= ksp->max_it) ksp->reason = KSP_DIVERGED_ITS; in KSPSolve_BCGS() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cr/ |
H A D | cr.c | 26 X = ksp->vec_sol; in KSPSolve_CR() 27 B = ksp->vec_rhs; in KSPSolve_CR() 28 R = ksp->work[0]; in KSPSolve_CR() 29 RT = ksp->work[1]; in KSPSolve_CR() 30 P = ksp->work[2]; in KSPSolve_CR() 31 AP = ksp->work[3]; in KSPSolve_CR() 72 ksp->its = 0; in KSPSolve_CR() 78 ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CR() 125 ksp->its++; in KSPSolve_CR() 131 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CR() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cg/cgne/ |
H A D | cgne.c | 39 if (ksp->calc_sings) { in KSPSetUp_CGNE() 93 ksp->its = 0; in KSPSolve_CGNE() 95 if (!ksp->guess_zero) { in KSPSolve_CGNE() 120 ksp->rnorm = dp; in KSPSolve_CGNE() 121 …ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for converg… in KSPSolve_CGNE() 126 ksp->its = i+1; in KSPSolve_CGNE() 155 KSPCheckDot(ksp,dpi); in KSPSolve_CGNE() 173 ksp->rnorm = dp; in KSPSolve_CGNE() 177 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CGNE() 189 if (i >= ksp->max_it) ksp->reason = KSP_DIVERGED_ITS; in KSPSolve_CGNE() [all …]
|
/dports/misc/ncnn/ncnn-20211208/tests/ |
H A D | test_pooling1d.cpp | 63 || test_pooling1d(9, 1, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling1d_0() 64 || test_pooling1d(9, 2, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling1d_0() 65 || test_pooling1d(9, 3, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling1d_0() 66 || test_pooling1d(9, 4, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 3, 0, 0, 0) in test_pooling1d_0() 67 || test_pooling1d(9, 7, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling1d_0() 68 || test_pooling1d(9, 8, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling1d_0() 69 || test_pooling1d(9, 15, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling1d_0() 98 || test_pooling1d(9, 1, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling1d_1() 99 || test_pooling1d(9, 2, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling1d_1() 100 || test_pooling1d(9, 3, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 1, 0, 0) in test_pooling1d_1() [all …]
|
H A D | test_pooling.cpp | 63 || test_pooling(9, 7, 1, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling_0() 64 || test_pooling(9, 7, 2, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling_0() 65 || test_pooling(9, 7, 3, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling_0() 66 || test_pooling(9, 7, 4, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 3, 0, 0, 0) in test_pooling_0() 67 || test_pooling(9, 7, 7, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling_0() 68 || test_pooling(9, 7, 8, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling_0() 69 || test_pooling(9, 7, 15, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling_0() 98 || test_pooling(9, 7, 1, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling_1() 99 || test_pooling(9, 7, 2, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling_1() 100 || test_pooling(9, 7, 3, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 1, 0, 0) in test_pooling_1() [all …]
|
H A D | test_pooling3d.cpp | 63 || test_pooling3d(9, 8, 7, 1, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling3d_0() 64 || test_pooling3d(9, 8, 7, 2, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling3d_0() 65 || test_pooling3d(9, 8, 7, 3, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling3d_0() 66 || test_pooling3d(9, 8, 7, 4, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 3, 0, 0, 0) in test_pooling3d_0() 67 || test_pooling3d(9, 8, 7, 7, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling3d_0() 68 || test_pooling3d(9, 8, 7, 8, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling3d_0() 69 || test_pooling3d(9, 8, 7, 15, 0, ksp[i][0], ksp[i][1], ksp[i][2], 0, 2, 0, 0, 0) in test_pooling3d_0() 98 || test_pooling3d(9, 8, 7, 1, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 0, 0, 0) in test_pooling3d_1() 99 || test_pooling3d(9, 8, 7, 2, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 1, 0, 0, 0) in test_pooling3d_1() 100 || test_pooling3d(9, 8, 7, 3, 1, ksp[i][0], ksp[i][1], ksp[i][2], 0, 0, 1, 0, 0) in test_pooling3d_1() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/lcd/ |
H A D | lcd.c | 47 X = ksp->vec_sol; in KSPSolve_LCD() 48 B = ksp->vec_rhs; in KSPSolve_LCD() 56 ksp->its = 0; in KSPSolve_LCD() 71 ksp->rnorm = rnorm; in KSPSolve_LCD() 74 ierr = (*ksp->converged)(ksp,0,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_LCD() 79 while (!ksp->reason && ksp->its < ksp->max_it) { in KSPSolve_LCD() 84 while (!ksp->reason && it < max_k && ksp->its < ksp->max_it) { in KSPSolve_LCD() 85 ksp->its++; in KSPSolve_LCD() 99 ierr = KSPMonitor(ksp,ksp->its,rnorm);CHKERRQ(ierr); in KSPSolve_LCD() 100 ierr = (*ksp->converged)(ksp,ksp->its,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_LCD() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/bcgs/fbcgs/ |
H A D | fbcgs.c | 31 X = ksp->vec_sol; in KSPSolve_FBCGS() 32 B = ksp->vec_rhs; in KSPSolve_FBCGS() 33 R = ksp->work[0]; in KSPSolve_FBCGS() 34 RP = ksp->work[1]; in KSPSolve_FBCGS() 70 ksp->its = 0; in KSPSolve_FBCGS() 71 ksp->rnorm = dp; in KSPSolve_FBCGS() 75 ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_FBCGS() 126 if (ksp->normtype != KSP_NORM_NONE && ksp->chknorm < i+2) { in KSPSolve_FBCGS() 134 ksp->its++; in KSPSolve_FBCGS() 139 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_FBCGS() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/gmres/ |
H A D | gmres.c | 144 ierr = KSPMonitor(ksp,ksp->its,res);CHKERRQ(ierr); in KSPGMRESCycle() 151 ierr = (*ksp->converged)(ksp,ksp->its,res,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPGMRESCycle() 152 while (!ksp->reason && it < max_k && ksp->its < ksp->max_it) { in KSPGMRESCycle() 155 ierr = KSPMonitor(ksp,ksp->its,res);CHKERRQ(ierr); in KSPGMRESCycle() 186 ksp->its++; in KSPGMRESCycle() 190 ierr = (*ksp->converged)(ksp,ksp->its,res,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPGMRESCycle() 207 if (it && (ksp->reason || ksp->its >= ksp->max_it)) { in KSPGMRESCycle() 209 ierr = KSPMonitor(ksp,ksp->its,res);CHKERRQ(ierr); in KSPGMRESCycle() 221 ierr = KSPGMRESBuildSoln(GRS(0),ksp->vec_sol,ksp->vec_sol,ksp,it-1);CHKERRQ(ierr); in KSPGMRESCycle() 237 ksp->its = 0; in KSPSolve_GMRES() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/tfqmr/ |
H A D | tfqmr.c | 23 X = ksp->vec_sol; in KSPSolve_TFQMR() 24 B = ksp->vec_rhs; in KSPSolve_TFQMR() 43 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dp; in KSPSolve_TFQMR() 45 ksp->its = 0; in KSPSolve_TFQMR() 47 ierr = KSPMonitor(ksp,0,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_TFQMR() 48 ierr = (*ksp->converged)(ksp,0,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_TFQMR() 69 ksp->its++; in KSPSolve_TFQMR() 97 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = dpest; in KSPSolve_TFQMR() 101 ierr = KSPMonitor(ksp,i+1,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_TFQMR() 102 ierr = (*ksp->converged)(ksp,i+1,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_TFQMR() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/gcr/ |
H A D | gcr.c | 44 ierr = (*ctx->modifypc)(ksp,ksp->its,ksp->rnorm,ctx->modifypc_ctx);CHKERRQ(ierr); in KSPSolve_GCR_cycle() 62 if (ksp->its > ksp->chknorm && ksp->normtype != KSP_NORM_NONE) { in KSPSolve_GCR_cycle() 67 ksp->its++; in KSPSolve_GCR_cycle() 71 ierr = KSPMonitor(ksp,ksp->its,norm_r);CHKERRQ(ierr); in KSPSolve_GCR_cycle() 73 if (ksp->its-1 > ksp->chknorm) { in KSPSolve_GCR_cycle() 74 ierr = (*ksp->converged)(ksp,ksp->its,norm_r,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_GCR_cycle() 78 if (ksp->its >= ksp->max_it) { in KSPSolve_GCR_cycle() 112 ierr = KSPMonitor(ksp,ksp->its,ksp->rnorm0);CHKERRQ(ierr); in KSPSolve_GCR() 113 ierr = (*ksp->converged)(ksp,ksp->its,ksp->rnorm0,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_GCR() 119 } while (ksp->its < ksp->max_it);CHKERRQ(ierr); in KSPSolve_GCR() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cgs/ |
H A D | cgs.c | 34 X = ksp->vec_sol; in KSPSolve_CGS() 35 B = ksp->vec_rhs; in KSPSolve_CGS() 36 R = ksp->work[0]; in KSPSolve_CGS() 37 RP = ksp->work[1]; in KSPSolve_CGS() 38 V = ksp->work[2]; in KSPSolve_CGS() 56 ksp->its = 0; in KSPSolve_CGS() 57 ksp->rnorm = dp; in KSPSolve_CGS() 61 ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CGS() 112 ksp->its++; in KSPSolve_CGS() 117 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_CGS() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/bicg/ |
H A D | bicg.c | 30 X = ksp->vec_sol; in KSPSolve_BiCG() 31 B = ksp->vec_rhs; in KSPSolve_BiCG() 32 Rl = ksp->work[0]; in KSPSolve_BiCG() 33 Zl = ksp->work[1]; in KSPSolve_BiCG() 34 Pl = ksp->work[2]; in KSPSolve_BiCG() 35 Rr = ksp->work[3]; in KSPSolve_BiCG() 62 ksp->its = 0; in KSPSolve_BiCG() 63 ksp->rnorm = dp; in KSPSolve_BiCG() 66 ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_BiCG() 118 ierr = (*ksp->converged)(ksp,i+1,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_BiCG() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cheby/ |
H A D | cheby.c | 411 ksp->its = 0; in KSPSolve_Chebyshev() 443 if (ksp->normtype) { in KSPSolve_Chebyshev() 460 ierr = (*ksp->converged)(ksp,0,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Chebyshev() 462 if (ksp->reason || ksp->max_it==0) { in KSPSolve_Chebyshev() 471 ksp->its = 1; in KSPSolve_Chebyshev() 476 ksp->its++; in KSPSolve_Chebyshev() 502 ierr = (*ksp->converged)(ksp,i,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Chebyshev() 523 if (!ksp->reason) { in KSPSolve_Chebyshev() 547 if (ksp->its >= ksp->max_it) { in KSPSolve_Chebyshev() 549 ierr = (*ksp->converged)(ksp,i,rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_Chebyshev() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/gmres/fgmres/ |
H A D | fgmres.c | 132 ierr = (*ksp->converged)(ksp,ksp->its,res_norm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPFGMRESCycle() 144 while (!ksp->reason && loc_it < max_k && ksp->its < ksp->max_it) { in KSPFGMRESCycle() 147 ierr = KSPMonitor(ksp,ksp->its,res_norm);CHKERRQ(ierr); in KSPFGMRESCycle() 213 ksp->its++; in KSPFGMRESCycle() 217 ierr = (*ksp->converged)(ksp,ksp->its,res_norm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPFGMRESCycle() 235 if (loc_it && (ksp->reason || ksp->its >= ksp->max_it)) { in KSPFGMRESCycle() 236 ierr = KSPMonitor(ksp,ksp->its,res_norm);CHKERRQ(ierr); in KSPFGMRESCycle() 251 ierr = KSPFGMRESBuildSoln(RS(0),ksp->vec_sol,ksp->vec_sol,ksp,loc_it-1);CHKERRQ(ierr); in KSPFGMRESCycle() 279 ksp->its = 0; in KSPSolve_FGMRES() 294 if (ksp->its >= ksp->max_it) break; in KSPSolve_FGMRES() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/preonly/ |
H A D | preonly.c | 4 static PetscErrorCode KSPSetUp_PREONLY(KSP ksp) in KSPSetUp_PREONLY() argument 21 ksp->its = 0; in KSPSolve_PREONLY() 22 ierr = KSP_PCApply(ksp,ksp->vec_rhs,ksp->vec_sol);CHKERRQ(ierr); in KSPSolve_PREONLY() 29 ksp->its = 1; in KSPSolve_PREONLY() 30 ksp->reason = KSP_CONVERGED_ITS; in KSPSolve_PREONLY() 46 ksp->its = 0; in KSPMatSolve_PREONLY() 54 ksp->its = 1; in KSPMatSolve_PREONLY() 55 ksp->reason = KSP_CONVERGED_ITS; in KSPMatSolve_PREONLY() 98 ksp->data = NULL; in KSPCreate_PREONLY() 105 ksp->ops->setfromoptions = NULL; in KSPCreate_PREONLY() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/minres/ |
H A D | minres.c | 51 ksp->its = 0; in KSPSolve_MINRES() 77 ksp->rnorm = 0.0; in KSPSolve_MINRES() 78 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = np; in KSPSolve_MINRES() 79 ierr = KSPLogResidualHistory(ksp,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_MINRES() 80 ierr = KSPMonitor(ksp,0,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_MINRES() 81 …ierr = (*ksp->converged)(ksp,0,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for conve… in KSPSolve_MINRES() 160 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = np; in KSPSolve_MINRES() 161 ierr = KSPLogResidualHistory(ksp,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_MINRES() 162 ierr = KSPMonitor(ksp,i+1,ksp->rnorm);CHKERRQ(ierr); in KSPSolve_MINRES() 163 …ierr = (*ksp->converged)(ksp,i+1,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for con… in KSPSolve_MINRES() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/cg/groppcg/ |
H A D | groppcg.c | 39 x = ksp->vec_sol; in KSPSolve_GROPPCG() 40 b = ksp->vec_rhs; in KSPSolve_GROPPCG() 41 r = ksp->work[0]; in KSPSolve_GROPPCG() 42 p = ksp->work[1]; in KSPSolve_GROPPCG() 43 s = ksp->work[2]; in KSPSolve_GROPPCG() 44 S = ksp->work[3]; in KSPSolve_GROPPCG() 45 z = ksp->work[4]; in KSPSolve_GROPPCG() 50 ksp->its = 0; in KSPSolve_GROPPCG() 86 …ierr = (*ksp->converged)(ksp,0,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); /* test for converg… in KSPSolve_GROPPCG() 132 ierr = (*ksp->converged)(ksp,i,dp,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPSolve_GROPPCG() [all …]
|
/dports/science/PETSc/petsc-3.14.1/src/ksp/ksp/impls/gmres/pgmres/ |
H A D | pgmres.c | 65 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle() 69 ierr = KSPLogResidualHistory(ksp,ksp->rnorm);CHKERRQ(ierr); in KSPPGMRESCycle() 70 ierr = KSPMonitor(ksp,ksp->its,ksp->rnorm);CHKERRQ(ierr); in KSPPGMRESCycle() 77 ierr = (*ksp->converged)(ksp,ksp->its,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPPGMRESCycle() 106 ksp->its++; in KSPPGMRESCycle() 107 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle() 110 ierr = (*ksp->converged)(ksp,ksp->its,ksp->rnorm,&ksp->reason,ksp->cnvP);CHKERRQ(ierr); in KSPPGMRESCycle() 113 ierr = KSPMonitor(ksp,ksp->its,ksp->rnorm);CHKERRQ(ierr); in KSPPGMRESCycle() 192 ierr = KSPPGMRESBuildSoln(RS(0),ksp->vec_sol,ksp->vec_sol,ksp,it-2);CHKERRQ(ierr); in KSPPGMRESCycle() 217 ksp->its = 0; in KSPSolve_PGMRES() [all …]
|