| /petsc/src/snes/tutorials/output/ |
| H A D | ex20_2d_p1_gmg_fcycle_rate.out | 1 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.698164 R^2 0.5911… 2 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.69301 R^2 0.541848 3 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.73243 R^2 0.805662 4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.696497 R^2 0.7674… 5 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.700915 R^2 0.6791… 6 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.760929 R^2 0.5625… 7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.728701 R^2 0.8390… 8 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.742781 R^2 0.6120… 9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.694817 R^2 0.7814… 10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.700569 R^2 0.7003… [all …]
|
| H A D | ex20_2d_p1_gmg_vcycle_rate.out | 1 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.774356 R^2 0.9846… 2 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.613993 R^2 0.9107… 3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.689662 R^2 0.5206… 4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.688226 R^2 0.5129… 5 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.729655 R^2 0.8322… 6 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.712464 R^2 0.8676… 7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.773593 R^2 0.9848… 8 … Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.6147 R^2 0.91186 9 … Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.6897 R^2 0.520871 10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.688241 R^2 0.5129… [all …]
|
| H A D | ex20_2d_p1_gmg_vcycle_adapt_rate.out | 1 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.774356 R^2 0.9846… 2 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.719169 R^2 0.6944… 3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.765723 R^2 0.6269… 4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.767943 R^2 0.5237… 5 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.729265 R^2 0.84772 6 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.72737 R^2 0.943032 7 …Linear mg_levels_3_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.770778 R^2 0.9834… 8 …Linear mg_levels_2_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.718876 R^2 0.7020… 9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.76585 R^2 0.625507 10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.768003 R^2 0.5232… [all …]
|
| H A D | ex20_2d_p1_scalable_rate.out | 1 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.785382 R^2 0.9647… 2 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.648211 R^2 0.9771… 3 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.736109 R^2 0.9663… 4 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.656634 R^2 0.9768… 5 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.764618 R^2 0.9750… 6 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.644732 R^2 0.9831… 7 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.776005 R^2 0.9734… 8 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.636778 R^2 0.9840… 9 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.75345 R^2 0.943944 10 …Linear mg_levels_1_ solve converged due to CONVERGED_ITS iterations 5 res rate 0.647287 R^2 0.9823… [all …]
|
| /petsc/src/binding/petsc4py/src/petsc4py/include/petsc4py/ |
| H A D | petsc4py.i | 21 int res = SWIG_ConvertPtr(obj, ptr, ty, flags); in SWIG_convert_ptr() local 22 if (!SWIG_IsOK(res)) { in SWIG_convert_ptr() 24 res = SWIG_ConvertPtr(_this, ptr, ty, flags); in SWIG_convert_ptr() 27 return res; in SWIG_convert_ptr() 95 long v; int res = SWIG_AsVal(long)(obj, &v); in PetscInt64() local 97 long long v; int res = SWIG_AsVal(long long)(obj, &v); in PetscInt64() local 99 if (SWIG_IsOK(res) && val) if (val) *val = %numeric_cast(v,PetscInt64); in PetscInt64() 100 return res; in PetscInt64() 139 double v; int res = SWIG_AsVal(double)(obj, &v); in SWIG_AsVal_dec() local 140 if (SWIG_IsOK(res) && val) if (val) *val = %numeric_cast(v,long double); in SWIG_AsVal_dec() [all …]
|
| /petsc/src/sys/tests/ |
| H A D | ex49f.F90 | 12 subroutine CompareIntegers(a, b, ctx, res) argument 17 integer :: res 20 res = -1 22 res = 0 24 res = 1
|
| H A D | ex15.cxx | 9 double res; in TestComplexOperators() local 34 res = PetscRealPartComplex(z); in TestComplexOperators() 35 … PetscCheck(PetscAbs(ans - res) < 1e-5, PETSC_COMM_SELF, PETSC_ERR_PLIB, "Expected %g, but get inc… in TestComplexOperators() 36 else ans = res; in TestComplexOperators()
|
| /petsc/src/vec/vec/impls/hypre/ |
| H A D | vhyp.c | 85 …v, int rw, HYPRE_MemoryLocation hmem, PetscScalar **ptr, PetscErrorCode (**res)(Vec, PetscScalar *… in VecGetArrayForHYPRE() 95 *res = NULL; in VecGetArrayForHYPRE() 101 *res = (PetscErrorCode (*)(Vec, PetscScalar **))VecRestoreArrayRead; in VecGetArrayForHYPRE() 105 *res = (PetscErrorCode (*)(Vec, PetscScalar **))VecRestoreArrayReadAndMemType; in VecGetArrayForHYPRE() 111 *res = VecRestoreArrayWrite; in VecGetArrayForHYPRE() 115 *res = VecRestoreArrayWriteAndMemType; in VecGetArrayForHYPRE() 121 *res = VecRestoreArray; in VecGetArrayForHYPRE() 125 *res = VecRestoreArrayAndMemType; in VecGetArrayForHYPRE()
|
| /petsc/src/ksp/ksp/impls/gmres/ |
| H A D | gmres.c | 91 PetscReal res, hapbnd, tt; in KSPGMRESCycle() local 97 PetscCall(VecNormalize(VEC_VV(0), &res)); in KSPGMRESCycle() 98 KSPCheckNorm(ksp, res); in KSPGMRESCycle() 101 …if ((ksp->rnorm > 0.0) && (PetscAbsReal(res - ksp->rnorm) > gmres->breakdowntol * gmres->rnorm0)) { in KSPGMRESCycle() 103 (double)ksp->rnorm, (double)res, (double)gmres->rnorm0); in KSPGMRESCycle() 104 …estart, residual norm at start of cycle %g\n", (double)ksp->rnorm, (double)res, (double)gmres->rno… in KSPGMRESCycle() 108 *GRS(0) = gmres->rnorm0 = res; in KSPGMRESCycle() 111 ksp->rnorm = res; in KSPGMRESCycle() 114 PetscCall(KSPLogResidualHistory(ksp, res)); in KSPGMRESCycle() 116 PetscCall(KSPMonitor(ksp, ksp->its, res)); in KSPGMRESCycle() [all …]
|
| /petsc/src/ksp/ksp/impls/gmres/lgmres/ |
| H A D | lgmres.c | 91 PetscReal res_norm, res; in KSPLGMRESCycle() local 125 res = res_norm; in KSPLGMRESCycle() 131 if (!res) { in KSPLGMRESCycle() 142 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPLGMRESCycle() 155 PetscCall((*ksp->converged)(ksp, ksp->its, res, &ksp->reason, ksp->cnvP)); in KSPLGMRESCycle() 158 PetscCall(KSPLogResidualHistory(ksp, res)); in KSPLGMRESCycle() 160 PetscCall(KSPMonitor(ksp, ksp->its, res)); in KSPLGMRESCycle() 207 PetscCall(KSPLGMRESUpdateHessenberg(ksp, loc_it, hapend, &res)); in KSPLGMRESCycle() 215 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPLGMRESCycle() 219 PetscCall((*ksp->converged)(ksp, ksp->its, res, &ksp->reason, ksp->cnvP)); in KSPLGMRESCycle() [all …]
|
| /petsc/src/snes/tests/ |
| H A D | ex21f.F90 | 47 Vec res, x 65 PetscCallA(MatCreateVecs(user%A, x, res, ierr)) 68 PetscCallA(SNESSetPicard(snes, res, FormFunction, user%A, user%A, FormJacobian, user, ierr)) 72 PetscCallA(VecDestroy(res, ierr))
|
| /petsc/src/ksp/ksp/impls/gmres/pgmres/ |
| H A D | pgmres.c | 20 PetscReal res_norm, res, newnorm; in KSPPGMRESCycle() local 28 res = res_norm; in KSPPGMRESCycle() 33 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle() 39 if (!res) { in KSPPGMRESCycle() 70 PetscCall(KSPPGMRESUpdateHessenberg(ksp, it - 2, &hapend, &res)); in KSPPGMRESCycle() 73 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPPGMRESCycle() 84 …D, "Reached happy break down, but convergence was not indicated. Residual norm = %g", (double)res); in KSPPGMRESCycle() 236 …c PetscErrorCode KSPPGMRESUpdateHessenberg(KSP ksp, PetscInt it, PetscBool *hapend, PetscReal *res) in KSPPGMRESUpdateHessenberg() argument 294 *res = PetscAbsScalar(rs[it + 1]); in KSPPGMRESUpdateHessenberg() 301 *res = 0.0; in KSPPGMRESUpdateHessenberg()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex55.c | 233 Vec res; in main() local 238 PetscCall(VecDuplicate(xx, &res)); in main() 239 PetscCall(MatMult(Amat, xx, res)); in main() 240 PetscCall(VecAXPY(bb, -1.0, res)); in main() 241 PetscCall(VecDestroy(&res)); in main()
|
| H A D | ex54.c | 176 Vec res; in main() local 190 PetscCall(VecDuplicate(xx, &res)); in main() 191 PetscCall(MatMult(Amat, xx, res)); in main() 192 PetscCall(VecAXPY(bb, -1.0, res)); in main() 193 PetscCall(VecDestroy(&res)); in main()
|
| H A D | ex56.c | 376 Vec res; in main() local 405 PetscCall(VecDuplicate(xx, &res)); in main() 406 PetscCall(MatMult(Amat, xx, res)); in main() 407 PetscCall(VecAXPY(bb, -1.0, res)); in main() 408 PetscCall(VecDestroy(&res)); in main()
|
| /petsc/src/ksp/ksp/impls/gmres/agmres/ |
| H A D | agmres.c | 400 PetscReal res; in KSPAGMRESCycle() local 405 res = ksp->rnorm; /* Norm of the initial residual vector */ in KSPAGMRESCycle() 406 if (!res) { in KSPAGMRESCycle() 412 PetscCall((*ksp->converged)(ksp, ksp->its, res, &ksp->reason, ksp->cnvP)); in KSPAGMRESCycle() 423 res = ksp->rnorm; in KSPAGMRESCycle() 427 PetscCall((*ksp->converged)(ksp, ksp->its, res, &ksp->reason, ksp->cnvP)); in KSPAGMRESCycle() 428 PetscCall(KSPLogResidualHistory(ksp, res)); in KSPAGMRESCycle() 429 PetscCall(KSPMonitor(ksp, ksp->its, res)); in KSPAGMRESCycle() 440 PetscReal res_old, res; in KSPSolve_AGMRES() local 473 res = ksp->rnorm; in KSPSolve_AGMRES() [all …]
|
| /petsc/src/tao/leastsquares/tutorials/ |
| H A D | tomography.c | 48 Vec x, res; /* solution, function res(x) = A*x-b */ in main() local 69 PetscCall(VecCreateSeq(PETSC_COMM_SELF, user.M, &res)); in main() 83 PetscCall(TaoSetResidualRoutine(tao, res, EvaluateResidual, (void *)&user)); in main() 121 PetscCall(VecDestroy(&res)); in main()
|
| /petsc/src/sys/tests/linter/ |
| H A D | testValidFunctionDocStrings.c | 515 PetscErrorCode testFunctionParmsSameName(int tao, double res, PetscErrorCode (*func)(int tao, doubl… in testFunctionParmsSameName() argument 545 PetscErrorCode testFunctionParmsSameNameInOut(int *tao, double res, PetscErrorCode (*func)(int tao,… in testFunctionParmsSameNameInOut() argument
|
| /petsc/src/ts/utils/ |
| H A D | dmplexts.c | 312 PetscReal res; in DMTSCheckResidual() local 323 PetscCall(VecNorm(r, NORM_2, &res)); in DMTSCheckResidual() 325 …PetscCheck(res <= tol, comm, PETSC_ERR_ARG_WRONG, "L_2 Residual %g exceeds tolerance %g", (double)… in DMTSCheckResidual() 327 *residual = res; in DMTSCheckResidual() 329 PetscCall(PetscPrintf(comm, "L_2 Residual: %g\n", (double)res)); in DMTSCheckResidual()
|
| /petsc/include/ |
| H A D | petsctao_deprecations.h | 33 … static inline PetscErrorCode TaoSetSeparableObjectiveRoutine(Tao tao, Vec res, PetscErrorCode (*f… in TaoSetSeparableObjectiveRoutine() argument 35 return TaoSetResidualRoutine(tao, res, func, ctx); in TaoSetSeparableObjectiveRoutine()
|
| /petsc/src/ksp/ksp/impls/gmres/dgmres/ |
| H A D | dgmres.c | 107 PetscReal res_norm, res, hapbnd, tt; in KSPDGMRESCycle() local 117 res = res_norm; in KSPDGMRESCycle() 122 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPDGMRESCycle() 128 if (!res) { in KSPDGMRESCycle() 135 res_old = res; in KSPDGMRESCycle() 175 PetscCall(KSPDGMRESUpdateHessenberg(ksp, it, hapend, &res)); in KSPDGMRESCycle() 180 if (ksp->normtype != KSP_NORM_NONE) ksp->rnorm = res; in KSPDGMRESCycle() 189 …D, "Reached happy break down, but convergence was not indicated. Residual norm = %g", (double)res); in KSPDGMRESCycle() 215 test = max_k * PetscLogReal(ksp->rtol / res) / PetscLogReal(res / res_old); in KSPDGMRESCycle() 364 …ic PetscErrorCode KSPDGMRESUpdateHessenberg(KSP ksp, PetscInt it, PetscBool hapend, PetscReal *res) in KSPDGMRESUpdateHessenberg() argument [all …]
|
| /petsc/src/tao/interface/ |
| H A D | taosolver_fg.c | 341 PetscErrorCode TaoSetResidualRoutine(Tao tao, Vec res, PetscErrorCode (*func)(Tao tao, Vec x, Vec r… in TaoSetResidualRoutine() argument 345 PetscValidHeaderSpecific(res, VEC_CLASSID, 2); in TaoSetResidualRoutine() 346 PetscCall(PetscObjectReference((PetscObject)res)); in TaoSetResidualRoutine() 348 tao->ls_res = res; in TaoSetResidualRoutine()
|
| /petsc/src/binding/petsc4py/demo/python_types/ |
| H A D | tao.py | 38 tao.monitor(f=f, res=gradient.norm())
|
| /petsc/src/snes/utils/ |
| H A D | dmplexsnes.c | 187 Vec res; in SNESMonitorFields() local 196 PetscCall(SNESGetFunction(snes, &res, NULL, NULL)); in SNESMonitorFields() 202 PetscCall(VecGetArrayRead(res, &r)); in SNESMonitorFields() 212 PetscCall(VecRestoreArrayRead(res, &r)); in SNESMonitorFields() 850 PetscReal res; in DMSNESCheckResidual() local 861 PetscCall(VecNorm(r, NORM_2, &res)); in DMSNESCheckResidual() 863 …PetscCheck(res <= tol, comm, PETSC_ERR_ARG_WRONG, "L_2 Residual %g exceeds tolerance %g", (double)… in DMSNESCheckResidual() 865 *residual = res; in DMSNESCheckResidual() 867 PetscCall(PetscPrintf(comm, "L_2 Residual: %g\n", (double)res)); in DMSNESCheckResidual()
|
| /petsc/src/sys/utils/ftn-custom/ |
| H A D | zsortsof.c | 13 void (*fcmp)(const void *a, const void *b, void *c, int *res);
|