| /petsc/src/ts/tests/ |
| H A D | ex4.c | 26 PetscReal epsilon; /* the diffusion coefficient */ member
63 data.epsilon = 0.1; in main()
277 PetscReal dx, dy, a, epsilon, xc, xl, xr, yl, yr; in RHSJacobian() local
286 epsilon = data->epsilon; in RHSJacobian()
288 xc = -2.0 * epsilon * (1.0 / (dx * dx) + 1.0 / (dy * dy)); in RHSJacobian()
289 xl = 0.5 * a / dx + epsilon / (dx * dx); in RHSJacobian()
290 xr = -0.5 * a / dx + epsilon / (dx * dx); in RHSJacobian()
291 yl = 0.5 * a / dy + epsilon / (dy * dy); in RHSJacobian()
292 yr = -0.5 * a / dy + epsilon / (dy * dy); in RHSJacobian()
407 PetscReal a, epsilon; in RHSFunction() local
[all …]
|
| H A D | ex27.c | 28 PetscReal epsilon; /* gaussian regularization parameter */ member
39 options->epsilon = -1.0; in ProcessOptions()
45 …ilon", "Mollifier regularization parameter", "ex27.c", options->epsilon, &options->epsilon, NULL)); in ProcessOptions()
178 if (user->epsilon < 0.) user->epsilon = 0.64 * pow(user->h, 1.98); in CreateParticles()
180 …t_FMT " L: %g h: %g eps: %g\n", user->N, (double)user->L, (double)user->h, (double)user->epsilon)); in CreateParticles()
289 PetscReal vc_l[3], L = ctx->L, h = ctx->h, epsilon = ctx->epsilon, init = 0.5 * h - L; in ComputeGradS() local
304 for (q = 0; q < Np; ++q) sum += Gaussian(dim, &velocity[q * dim], epsilon, vc_l); in ComputeGradS()
306 … = 0; d < dim; ++d) integral[d] += (-1. / epsilon) * PetscAbsReal(vp[d] - vc_l[d]) * Gaussian(dim,… in ComputeGradS()
|
| /petsc/src/tao/leastsquares/impls/brgn/ |
| H A D | brgn.c | 143 PetscCall(VecShift(gn->y_work, gn->epsilon * gn->epsilon)); in GNObjectiveGradientEval()
147 *fcn += gn->lambda * (yESum - K * gn->epsilon); in GNObjectiveGradientEval()
194 PetscCall(VecShift(gn->y_work, gn->epsilon * gn->epsilon)); in GNComputeHessian()
199 PetscCall(VecScale(gn->diag, gn->epsilon * gn->epsilon)); in GNComputeHessian()
358 …||x||_1 = sum(sqrt(x.^2+epsilon^2)-epsilon) (default 1e-6)", "", gn->epsilon, &gn->epsilon, NULL)); in TaoSetFromOptions_BRGN()
385 PetscCall(PetscViewerASCIIPrintf(viewer, "L1 smooth epsilon: %g\n", (double)gn->epsilon)); in TaoView_BRGN()
587 PetscErrorCode TaoBRGNSetL1SmoothEpsilon(Tao tao, PetscReal epsilon) in TaoBRGNSetL1SmoothEpsilon() argument
591 PetscValidLogicalCollectiveReal(tao, epsilon, 2); in TaoBRGNSetL1SmoothEpsilon()
592 PetscTryMethod((PetscObject)tao, "TaoBRGNSetL1SmoothEpsilon_C", (Tao, PetscReal), (tao, epsilon)); in TaoBRGNSetL1SmoothEpsilon()
596 static PetscErrorCode TaoBRGNSetL1SmoothEpsilon_BRGN(Tao tao, PetscReal epsilon) in TaoBRGNSetL1SmoothEpsilon_BRGN() argument
[all …]
|
| H A D | brgn.h | 29 …PetscReal lambda, epsilon, fc_old; /* lambda is regularizer w… member
|
| /petsc/src/mat/tests/ |
| H A D | ex246.cxx | 49 PetscReal *coords, *gcoords, norm, epsilon, relative; in main() local
60 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mat_htool_epsilon", &epsilon, NULL)); in main()
91 …) <= epsilon, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "||A(!symmetric)-A(symmetric)|| = %g (> %g)", (dou… in main()
|
| H A D | ex241.c | 43 …PetscReal *coords, *gcoords, *scoords, *gscoords, *ctx[2], norm, epsilon = PetscSqrtReal(PE… in main() local
58 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mat_htool_epsilon", &epsilon, NULL)); in main()
110 …if (PetscAbsReal(epsilon) >= PETSC_SMALL) { /* when there is compression, it is more difficult to … in main()
125 …) <= epsilon, PETSC_COMM_WORLD, PETSC_ERR_PLIB, "||A(!symmetric)-A(symmetric)|| = %g (> %g)", (dou… in main()
|
| H A D | ex268.c | 27 PetscReal *coords, *gcoords, norm, epsilon; in main() local
40 PetscCall(PetscOptionsGetReal(NULL, NULL, "-mat_htool_epsilon", &epsilon, NULL)); in main()
|
| /petsc/src/binding/petsc4py/demo/legacy/kspsolve/ |
| H A D | petsc-cg.tex | 3 & cg(A,x,b,i_{max},\epsilon): \\
10 & \quad\quad\qquad \delta_{ } > \delta_{0} \epsilon^2 \text{ do} :\\
|
| /petsc/src/ts/tutorials/ |
| H A D | ex10.c | 412 const PetscScalar epsilon = x[i].T * PETSC_SQRT_MACHINE_EPSILON; in RDIJacobian_FD() local
418 n1.T = x[i].T + epsilon; in RDIJacobian_FD()
423 Em_TT = (dEm1.T - dEm.T) / epsilon; in RDIJacobian_FD()
871 PetscScalar epsilon; in RDTestDifferentiation() local
875 epsilon = 1e-8; in RDTestDifferentiation()
878 PetscScalar T0 = 1000., T1 = T0 * (1. + epsilon), Em0, Em1; in RDTestDifferentiation()
882 n.E = 1. + epsilon; in RDTestDifferentiation()
885 fdEm.E = (Em1 - Em0) / epsilon; in RDTestDifferentiation()
889 fdEm.T = (Em1 - Em0) / (T0 * epsilon); in RDTestDifferentiation()
901 n.E = 1. + epsilon; in RDTestDifferentiation()
[all …]
|
| H A D | ex34.c | 25 PetscReal epsilon; /* inverse of seismic ratio, B-A / A */ member
56 f[i].th = -hx * (u[i].v + 1) * (u[i].th + (1 + user->epsilon) * PetscLogScalar(u[i].v + 1)); in FormRHSFunction()
228 user.epsilon = 0.1; in main()
233 …l(PetscOptionsReal("-epsilon", "Inverse of seismic ratio", "", user.epsilon, &user.epsilon, NULL)); in main()
|
| /petsc/src/dm/dt/tests/ |
| H A D | ex3.c | 109 const PetscReal epsilon = 2.2204460492503131e-16; in main() local
111 const PetscReal epsilon = 2500. * PETSC_MACHINE_EPSILON; in main() local
134 …if (PetscAbsReal(integral - analytic[f]) > PetscMax(epsilon, PetscPowRealInt(10.0, -digits)) || Pe… in main()
|
| /petsc/src/snes/interface/ |
| H A D | snesj.c | 54 PetscReal amax, epsilon = PETSC_SQRT_MACHINE_EPSILON; in SNESComputeJacobianDefault() local
68 …Real(((PetscObject)snes)->options, ((PetscObject)snes)->prefix, "-snes_test_err", &epsilon, NULL)); in SNESComputeJacobianDefault()
116 dx *= epsilon; in SNESComputeJacobianDefault()
|
| /petsc/src/tao/unconstrained/impls/ntr/ |
| H A D | ntr.c | 139 if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) { in TaoSolve_NTR()
307 … if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) { in TaoSolve_NTR()
354 … if ((PetscAbsScalar(actred) <= tr->epsilon) && (PetscAbsScalar(prered) <= tr->epsilon)) { in TaoSolve_NTR()
491 …lerance used when computing actual and predicted reduction", "", tr->epsilon, &tr->epsilon, NULL)); in TaoSetFromOptions_NTR()
590 tr->epsilon = 1.0e-6; in TaoCreate_NTR()
|
| H A D | ntrimpl.h | 87 PetscReal epsilon; /* tolerance used when computing actred/prered */ member
|
| /petsc/src/tao/bound/impls/bqnls/ |
| H A D | bqnls.c | 48 …rance used when computing actual and predicted reduction", "", bnk->epsilon, &bnk->epsilon, NULL)); in TaoSetFromOptions_BQNLS()
|
| /petsc/src/ksp/pc/impls/spai/ |
| H A D | ispai.c | 45 double epsilon; /* tolerance */ member
90 …PetscCallExternal(bspai, ispai->B, &ispai->M, stdout, ispai->epsilon, ispai->nbsteps, ispai->max, … in PCSetUp_SPAI()
147 PetscCall(PetscViewerASCIIPrintf(viewer, " epsilon %g\n", ispai->epsilon)); in PCView_SPAI()
164 ispai->epsilon = (double)epsilon1; in PCSPAISetEpsilon_SPAI()
442 …PetscCall(PetscOptionsReal("-pc_spai_epsilon", "", "PCSPAISetEpsilon", ispai->epsilon, &epsilon1, … in PCSetFromOptions_SPAI()
502 ispai->epsilon = .4; in PCCreate_SPAI()
|
| /petsc/src/ts/tutorials/advection-diffusion-reaction/ |
| H A D | ex4.c | 33 PetscScalar epsilon, delta, alpha, beta, gamma, kappa, lambda, mu, cstar; member
55 appctx.epsilon = 1.0e-3; in main()
203 …f[i].rho = udot[i].rho - appctx->epsilon * rhoxx + kcxrhox - appctx->mu * PetscAbsScalar(rho) * (1… in IFunction()
|
| /petsc/src/ksp/ksp/tutorials/output/ |
| H A D | ex82_1.out | 33 epsilon: 0.01
82 epsilon: 0.01
127 epsilon: 0.01
|
| /petsc/src/tao/unconstrained/impls/bmrm/ |
| H A D | bmrm.c | 306 PetscReal jtwt = 0.0, max_jtwt, pre_epsilon, epsilon, jw, min_jw; in TaoSolve_BMRM() local
320 epsilon = 0.0; in TaoSolve_BMRM()
409 pre_epsilon = epsilon; in TaoSolve_BMRM()
410 epsilon = min_jw - jtwt; in TaoSolve_BMRM()
413 if (innerSolverTol > epsilon) innerSolverTol = epsilon; in TaoSolve_BMRM()
417 if (pre_epsilon < epsilon) innerSolverTol *= 0.2; in TaoSolve_BMRM()
423 PetscCall(TaoLogConvergenceHistory(tao, min_jw, epsilon, 0.0, tao->ksp_its)); in TaoSolve_BMRM()
424 PetscCall(TaoMonitor(tao, tao->niter, min_jw, epsilon, 0.0, tao->step)); in TaoSolve_BMRM()
|
| /petsc/src/tao/unconstrained/impls/ntl/ |
| H A D | ntl.c | 133 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
312 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
362 if ((PetscAbsScalar(actred) <= tl->epsilon) && (PetscAbsScalar(prered) <= tl->epsilon)) { in TaoSolve_NTL()
676 …lerance used when computing actual and predicted reduction", "", tl->epsilon, &tl->epsilon, NULL)); in TaoSetFromOptions_NTL()
809 tl->epsilon = 1.0e-6; in TaoCreate_NTL()
|
| H A D | ntlimpl.h | 112 PetscReal epsilon; /* tolerance used when computing ared/pred */ member
|
| /petsc/src/mat/impls/htool/ |
| H A D | htool.hpp | 42 PetscReal epsilon; member
|
| /petsc/src/tao/unconstrained/impls/nls/ |
| H A D | nls.c | 159 … if ((PetscAbsScalar(actred) <= nlsP->epsilon) && (PetscAbsScalar(prered) <= nlsP->epsilon)) { in TaoSolve_NLS()
576 … if ((PetscAbsScalar(actred) <= nlsP->epsilon) && (PetscAbsScalar(prered) <= nlsP->epsilon)) { in TaoSolve_NLS()
625 … if ((PetscAbsScalar(actred) <= nlsP->epsilon) && (PetscAbsScalar(prered) <= nlsP->epsilon)) { in TaoSolve_NLS()
781 …nce used when computing actual and predicted reduction", "", nlsP->epsilon, &nlsP->epsilon, NULL)); in TaoSetFromOptions_NLS()
949 nlsP->epsilon = 1.0e-6; in TaoCreate_NLS()
|
| H A D | nlsimpl.h | 156 PetscReal epsilon; /* tolerance used when computing ared/pred */ member
|
| /petsc/src/tao/bound/impls/bncg/ |
| H A D | bncg.h | 30 PetscReal epsilon; /* Machine precision unless changed by user */ member
|