Home
last modified time | relevance | path

Searched refs:g2 (Results 1 – 25 of 34) sorted by relevance

12

/petsc/src/dm/impls/forest/tests/
H A Dex1.c11 Vec g, g2; in main() local
51 PetscCall(DMCreateGlobalVector(forest, &g2)); in main()
52 PetscCall(PetscObjectSetName((PetscObject)g2, "g")); in main()
54 PetscCall(VecLoad(g2, viewer)); in main()
58 PetscCall(VecAXPY(g2, -1.0, g)); in main()
59 PetscCall(VecNorm(g2, NORM_INFINITY, &diff)); in main()
63 PetscCall(VecDestroy(&g2)); in main()
H A Dex3.c12 Vec g = NULL, g2 = NULL; in main() local
124 PetscCall(DMCreateGlobalVector(forest, &g2)); in main()
125 PetscCall(PetscObjectSetName((PetscObject)g2, "g")); in main()
126 PetscCall(VecZeroEntries(g2)); in main()
130 PetscCall(VecLoad(g2, viewer)); in main()
134 PetscCall(VecAXPY(g2, -1.0, g)); in main()
135 PetscCall(VecNorm(g2, NORM_INFINITY, &nrm)); in main()
139 PetscCall(VecDestroy(&g2)); in main()
/petsc/src/tao/interface/
H A Dtaosolver_fg.c29 Vec g2, g3; in TaoTestGradient() local
65 PetscCall(VecDuplicate(x, &g2)); in TaoTestGradient()
69 PetscCall(TaoDefaultComputeGradient(tao, x, g2, NULL)); in TaoTestGradient()
71 PetscCall(VecNorm(g2, NORM_2, &fdnorm)); in TaoTestGradient()
73 PetscCall(VecNorm(g2, NORM_INFINITY, &fdmax)); in TaoTestGradient()
75 PetscCall(VecDot(g1, g2, &dot)); in TaoTestGradient()
77 PetscCall(VecAXPY(g3, -1.0, g2)); in TaoTestGradient()
88 PetscCall(VecView(g2, mviewer)); in TaoTestGradient()
92 PetscCall(VecDestroy(&g2)); in TaoTestGradient()
/petsc/src/dm/dt/interface/
H A Ddtweakform.c556 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n2, void (***g2)(PetscInt, PetscInt,… in PetscWeakFormGetJacobian()
563 …mGetFunction_Private(wf, wf->form[PETSC_WF_G2], label, val, find, part, n2, (void (***)(void))g2)); in PetscWeakFormGetJacobian()
568 …, const PetscReal[], PetscInt, const PetscScalar[], PetscScalar[]), void (*g2)(PetscInt, PetscInt,… in PetscWeakFormAddJacobian()
575 …WeakFormAddFunction_Private(wf, wf->form[PETSC_WF_G2], label, val, find, part, (PetscVoidFn *)g2)); in PetscWeakFormAddJacobian()
580 …eal[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt n2, void (**g2)(PetscInt, PetscInt,… in PetscWeakFormSetJacobian()
587 …ormSetFunction_Private(wf, wf->form[PETSC_WF_G2], label, val, find, part, n2, (PetscVoidFn **)g2)); in PetscWeakFormSetJacobian()
592 …Real[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt i2, void (*g2)(PetscInt, PetscInt,… in PetscWeakFormSetIndexJacobian()
599 …etIndexFunction_Private(wf, wf->form[PETSC_WF_G2], label, val, find, part, i2, (PetscVoidFn *)g2)); in PetscWeakFormSetIndexJacobian()
619 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n2, void (***g2)(PetscInt, PetscInt,… in PetscWeakFormGetJacobianPreconditioner()
626 …GetFunction_Private(wf, wf->form[PETSC_WF_GP2], label, val, find, part, n2, (void (***)(void))g2)); in PetscWeakFormGetJacobianPreconditioner()
[all …]
H A Ddtds.c519 … &prob->g2, NsMax * NsMax * NqMax * NcMax * NcMax * dimEmbed * dimEmbed, &prob->g3)); in PetscDSSetUp()
534 PetscCall(PetscFree6(prob->f0, prob->f1, prob->g0, prob->g1, prob->g2, prob->g3)); in PetscDSDestroyStructs_Static()
1508 … PetscInt g, PetscPointJacFn **g0, PetscPointJacFn **g1, PetscPointJacFn **g2, PetscPointJacFn **g… in PetscDSGetJacobian() argument
1520 *g2 = tmp2 ? tmp2[0] : NULL; in PetscDSGetJacobian()
1551 … f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2, PetscPointJacFn *g3) in PetscDSSetJacobian() argument
1557 if (g2) PetscValidFunction(g2, 6); in PetscDSSetJacobian()
1561 PetscCall(PetscWeakFormSetIndexJacobian(ds->wf, NULL, 0, f, g, 0, 0, g0, 0, g1, 0, g2, 0, g3)); in PetscDSSetJacobian()
1646 … PetscInt g, PetscPointJacFn **g0, PetscPointJacFn **g1, PetscPointJacFn **g2, PetscPointJacFn **g… in PetscDSGetJacobianPreconditioner() argument
1658 *g2 = tmp2 ? tmp2[0] : NULL; in PetscDSGetJacobianPreconditioner()
1693 … f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2, PetscPointJacFn *g3) in PetscDSSetJacobianPreconditioner() argument
[all …]
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h131 …p}, x1[3] = {a_x1[0], a_x1[1], a_x1[2]}, c02 = c0 * c0, g1 = GAMMA3(x1, c02), g2 = GAMMA3(x2, c02); in LandauTensor3DRelativistic() local
147 …tt = 2. * u1u2 * (1. - g1 * g2) + (u12 * u22 + u1u2 * u1u2) / c02; // these two terms are about t… in LandauTensor3DRelativistic()
151 …fact = -rsq / (g1 * g2 * PetscSqrtReal(wsq)); /* flip sign. papers use du/dt = C, PETSc uses form … in LandauTensor3DRelativistic()
/petsc/src/dm/dt/fe/impls/basic/
H A Dfebasic.c757 …PetscScalar *g0 = NULL, *g1 = NULL, *g2 = NULL, *g3 = NULL, *u, *u_t = NULL, *u_x, *a, *a_x,… in PetscFEIntegrateJacobian_Basic() local
796 …PetscCall(PetscDSGetWeakFormArrays(rds, NULL, NULL, n0 ? &g0 : NULL, n1 ? &g1 : NULL, n2 ? &g2 : N… in PetscFEIntegrateJacobian_Basic()
863 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateJacobian_Basic()
864 …ff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, constants, g2); in PetscFEIntegrateJacobian_Basic()
865 for (PetscInt c = 0; c < NcI * NcJ * dE; ++c) g2[c] *= w; in PetscFEIntegrateJacobian_Basic()
873 …sisReal, basisDerReal, cT[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, ctotDim, offset… in PetscFEIntegrateJacobian_Basic()
906 …PetscScalar *g0, *g1, *g2, *g3, *u, *u_t = NULL, *u_x, *a, *a_x, *basisReal, *basisDerReal,… in PetscFEIntegrateBdJacobian_Basic() local
943 PetscCall(PetscDSGetWeakFormArrays(ds, NULL, NULL, &g0, &g1, &g2, &g3)); in PetscFEIntegrateBdJacobian_Basic()
962 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
1023 PetscCall(PetscArrayzero(g2, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
[all …]
/petsc/src/tao/unconstrained/tutorials/
H A Drosenbrock4.h53 PetscScalar g2 = alpha * d * 2.0; in RosenbrockGradient() local
55 g[0] = -2.0 * x_1 * g2 - 2.0 * e; in RosenbrockGradient()
56 g[1] = g2; in RosenbrockGradient()
66 PetscScalar g2 = ad * 2.0; in RosenbrockObjectiveGradient() local
68 g[0] = -2.0 * x_1 * g2 - 2.0 * e; in RosenbrockObjectiveGradient()
69 g[1] = g2; in RosenbrockObjectiveGradient()
78 PetscScalar g2 = alpha * d * 2.0; in RosenbrockHessian() local
81 h[0] = -2.0 * (g2 + x_1 * h2) + 2.0; in RosenbrockHessian()
/petsc/src/tao/pde_constrained/impls/lcl/
H A Dlcl.h52 Vec g1, g2; member
H A Dlcl.c45 PetscCall(VecDestroy(&lclP->g2)); in TaoDestroy_LCL()
152 PetscCall(VecDuplicate(lclP->V, &lclP->g2)); in TaoSetup_LCL()
514 PetscCall(MatMultTranspose(tao->jacobian_design, lclP->lambda, lclP->g2)); in TaoSolve_LCL()
515 PetscCall(VecAXPY(lclP->g2, -1.0, lclP->GV)); in TaoSolve_LCL()
517 PetscCall(VecScale(lclP->g2, -1.0)); in TaoSolve_LCL()
520 PetscCall(MatLMVMUpdate(lclP->R, lclP->V, lclP->g2)); in TaoSolve_LCL()
/petsc/src/tao/tutorials/
H A Dex1.c76 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_ua() argument
79 for (d = 0; d < dim; ++d) g2[d] = u_x[dim * 2 + d]; in g2_ua()
115 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_la() argument
118 for (d = 0; d < dim; ++d) g2[d] = u_x[d]; in g2_la()
/petsc/include/petsc/private/
H A Dpetscdsimpl.h143 PetscScalar *g0, *g1, *g2, *g3; /* Point evaluations of weak form Jacobian integrands */ member
/petsc/src/ts/impls/arkimex/
H A Darkimex.c482 #define g2 PetscSqr(g) in TSARKIMEXRegisterAll() macro
483 #define g3 g *g2 in TSARKIMEXRegisterAll()
484 #define g4 PetscSqr(g2) in TSARKIMEXRegisterAll()
489 #define b3 (RC(1.0) - RC(6.0) * g + RC(6.0) * g2) / (RC(3.0) * c3 * (c3 - RC(2.0) * g)) in TSARKIMEXRegisterAll()
493 …g2)) + (RC(3.0) - RC(27.0) * g + RC(68.0) * g2 - RC(55.0) * g3 + RC(21.0) * g4 - RC(6.0) * g5) / (… in TSARKIMEXRegisterAll()
494 …g - RC(68.0) * g2 + RC(79.0) * g3 - RC(33.0) * g4 + RC(12.0) * g5) / (RC(2.0) * (RC(2.0) * g - c3)… in TSARKIMEXRegisterAll()
495 …#define bh4 -RC(3.0) * g2 * (-RC(1.0) + RC(4.0) * g - RC(2.0) * g2 + g3) / (RC(1.0) - RC(6.0) * g … in TSARKIMEXRegisterAll()
514 #undef g2 in TSARKIMEXRegisterAll()
/petsc/src/ts/tutorials/
H A Dex47.c110 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_ibp_phi() argument
113 for (d = 0; d < dim; ++d) g2[d] = a[d]; in g2_ibp_phi()
H A Dex46.c173 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_up() argument
176 for (d = 0; d < dim; ++d) g2[d * dim + d] = -1.0; /* \frac{\partial\psi^{u_d}}{\partial x_d} */ in g2_up()
/petsc/src/dm/dt/tests/
H A Dex10.c96 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_vp() argument
100 for (c = 0; c < dim; ++c) g2[c * dim + c] = -1.0; in g2_vp()
/petsc/src/snes/tutorials/
H A Dex71.c148 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_up() argument
151 for (d = 0; d < dim; ++d) g2[d * dim + d] = -1.0; /* \frac{\partial\psi^{u_d}}{\partial x_d} */ in g2_up()
H A Dex24.c191 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_qu() argument
194 for (d = 0; d < dim; ++d) g2[d * dim + d] = 1.0; in g2_qu()
H A Dex62.c77 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_up() argument
80 for (d = 0; d < dim; ++d) g2[d * dim + d] = -1.0; /* -< \nabla\cdot v, p > */ in g2_up()
H A Dex77.c225 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_up_3d() argument
232 for (d = 0; d < dim; ++d) g2[compI * dim + d] = cofu_x[compI * dim + d]; in g2_up_3d()
H A Dex76.c243 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_vp() argument
246 for (d = 0; d < dim; ++d) g2[d * dim + d] = -1.0; in g2_vp()
/petsc/src/dm/dt/fe/interface/
H A Dfe.c2611 const PetscScalar *G = g2 + (fc * (_NcJ) + gc) * _dE; \
2645 …*fegeom, const PetscScalar g0[], const PetscScalar g1[], const PetscScalar g2[], const PetscScalar… in PetscFEUpdateElementMat_Internal() argument
2731 if (PetscUnlikely(g2)) { in PetscFEUpdateElementMat_Internal()
2754 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ … in PetscFEUpdateElementMat_Internal()
2800 …*fegeom, const PetscScalar g0[], const PetscScalar g1[], const PetscScalar g2[], const PetscScalar… in PetscFEUpdateElementMat_Hybrid_Internal() argument
2851 … elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g2[(fc * NcJ + gc) * dE + df] * tmpBasisJ[gidx]; in PetscFEUpdateElementMat_Hybrid_Internal()
2899 …etscInt dimEmbed, const PetscInt g0[], const PetscInt g1[], const PetscInt g2[], const PetscInt g3…
2908 PetscBool has_g2 = g2 ? PETSC_TRUE : PETSC_FALSE;
2923 if (g2[i]) g2_i += NbI * NbJ;
2973 if (g2[(fc * NcJ + gc) * dE + df]) {
/petsc/src/ts/utils/dmplexlandau/kokkos/
H A Dlandau.kokkos.cxx598 g2_scr_t g2(team.team_scratch(jac_shared_level), dim, loc_Nf, Nq); in LandauKokkosJacobian() local
679 g2(d, fieldA, myQi) = 0; in LandauKokkosJacobian()
681 g2(d, fieldA, myQi) += invJ[d * dim + d2] * gg2(d2, fieldA, myQi); in LandauKokkosJacobian()
688 g2(d, fieldA, myQi) *= wj; in LandauKokkosJacobian()
702 t += DIq[f * dim + d] * g2(d, fieldA, qj) * BJq[g]; in LandauKokkosJacobian()
/petsc/src/snes/tests/
H A Dex15.c175 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g2[]) in g2_qphi() argument
177 for (PetscInt d = 0; d < dim; ++d) g2[d * dim + d] = 1.0; in g2_qphi()
/petsc/src/snes/interface/
H A Dsnes.c2661 Vec x, g1, g2, g3; in SNESTestFunction() local
2698 PetscCall(VecDuplicate(x, &g2)); in SNESTestFunction()
2701 PetscCall(SNESComputeFunction_FD(snes, x, g2)); in SNESTestFunction()
2703 PetscCall(VecNorm(g2, NORM_2, &fdnorm)); in SNESTestFunction()
2705 PetscCall(VecNorm(g2, NORM_INFINITY, &fdmax)); in SNESTestFunction()
2707 PetscCall(VecDot(g1, g2, &dot)); in SNESTestFunction()
2709 PetscCall(VecAXPY(g3, -1.0, g2)); in SNESTestFunction()
2720 PetscCall(VecView(g2, mviewer)); in SNESTestFunction()
2725 PetscCall(VecDestroy(&g2)); in SNESTestFunction()

12