Home
last modified time | relevance | path

Searched refs:g3 (Results 1 – 25 of 44) sorted by relevance

12

/petsc/src/snes/tutorials/
H A Dex56.c75 void g3_uu_3d_private(PetscScalar g3[], const PetscReal mu, const PetscReal lambda) in g3_uu_3d_private() argument
78 g3[0] += lambda; in g3_uu_3d_private()
79 g3[0] += mu; in g3_uu_3d_private()
80 g3[0] += mu; in g3_uu_3d_private()
81 g3[4] += lambda; in g3_uu_3d_private()
82 g3[8] += lambda; in g3_uu_3d_private()
83 g3[10] += mu; in g3_uu_3d_private()
84 g3[12] += mu; in g3_uu_3d_private()
85 g3[20] += mu; in g3_uu_3d_private()
86 g3[24] += mu; in g3_uu_3d_private()
[all …]
H A Dex20.c54 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
57 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
H A Dex17.c314 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_vlap_uu() argument
320 for (d = 0; d < dim; ++d) g3[((c * Nc + c) * dim + d) * dim + d] = 1.0; in g3_vlap_uu()
330 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_elas_uu() argument
338 g3[((c * Nc + c) * dim + d) * dim + d] += mu; in g3_elas_uu()
339 g3[((c * Nc + d) * dim + d) * dim + c] += mu; in g3_elas_uu()
340 g3[((c * Nc + d) * dim + c) * dim + d] += lambda; in g3_elas_uu()
H A Dex12.c178 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
181 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
283 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_analytic_uu() argument
286 for (d = 0; d < dim; ++d) g3[d * dim + d] = x[0] + x[1]; in g3_analytic_uu()
289 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_field_uu() argument
292 for (d = 0; d < dim; ++d) g3[d * dim + d] = a[0]; in g3_field_uu()
328 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_analytic_nonlinear_uu() argument
334 g3[d * dim + d] = 0.5 * nu; in g3_analytic_nonlinear_uu()
335 for (e = 0; e < dim; ++e) g3[d * dim + e] += u_x[d] * u_x[e]; in g3_analytic_nonlinear_uu()
H A Dex26.c32 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
35 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
H A Dex34.c89 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
91 for (PetscInt d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
H A Dex23.c43 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
46 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
H A Dex76.c249 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_vu() argument
257 g3[((c * Nc + c) * dim + d) * dim + d] += nu; // gradU in g3_vu()
258 g3[((c * Nc + d) * dim + d) * dim + c] += nu; // gradU transpose in g3_vu()
275 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_wT() argument
280 for (d = 0; d < dim; ++d) g3[d * dim + d] = alpha; in g3_wT()
H A Dex36.c87 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
90 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
105 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_oscillatory_uu() argument
109 g3[0] = 1. / (2. + PetscCosReal(2. * PETSC_PI * x[0] / eps)); in g3_oscillatory_uu()
H A Dex62.c83 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
91 g3[((c * Nc + c) * dim + d) * dim + d] += mu; /* < \nabla v, \nabla u > */ in g3_uu()
92 g3[((c * Nc + d) * dim + d) * dim + c] += mu; /* < \nabla v, {\nabla u}^T > */ in g3_uu()
H A Dex71.c156 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
163 for (d = 0; d < dim; ++d) g3[((c * Nc + c) * dim + d) * dim + d] = nu; in g3_uu()
/petsc/src/tao/interface/
H A Dtaosolver_fg.c29 Vec g2, g3; in TaoTestGradient() local
66 PetscCall(VecDuplicate(x, &g3)); in TaoTestGradient()
76 PetscCall(VecCopy(g1, g3)); in TaoTestGradient()
77 PetscCall(VecAXPY(g3, -1.0, g2)); in TaoTestGradient()
78 PetscCall(VecNorm(g3, NORM_2, &diffnorm)); in TaoTestGradient()
79 PetscCall(VecNorm(g3, NORM_INFINITY, &diffmax)); in TaoTestGradient()
90 PetscCall(VecView(g3, mviewer)); in TaoTestGradient()
93 PetscCall(VecDestroy(&g3)); in TaoTestGradient()
/petsc/src/dm/dt/interface/
H A Ddtweakform.c556 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n3, void (***g3)(PetscInt, PetscInt,… in PetscWeakFormGetJacobian()
564 …mGetFunction_Private(wf, wf->form[PETSC_WF_G3], label, val, find, part, n3, (void (***)(void))g3)); in PetscWeakFormGetJacobian()
568 …, const PetscReal[], PetscInt, const PetscScalar[], PetscScalar[]), void (*g3)(PetscInt, PetscInt,… in PetscWeakFormAddJacobian()
576 …WeakFormAddFunction_Private(wf, wf->form[PETSC_WF_G3], label, val, find, part, (PetscVoidFn *)g3)); in PetscWeakFormAddJacobian()
580 …eal[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt n3, void (**g3)(PetscInt, PetscInt,… in PetscWeakFormSetJacobian()
588 …ormSetFunction_Private(wf, wf->form[PETSC_WF_G3], label, val, find, part, n3, (PetscVoidFn **)g3)); in PetscWeakFormSetJacobian()
592 …Real[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt i3, void (*g3)(PetscInt, PetscInt,… in PetscWeakFormSetIndexJacobian()
600 …etIndexFunction_Private(wf, wf->form[PETSC_WF_G3], label, val, find, part, i3, (PetscVoidFn *)g3)); in PetscWeakFormSetIndexJacobian()
619 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n3, void (***g3)(PetscInt, PetscInt,… in PetscWeakFormGetJacobianPreconditioner()
627 …GetFunction_Private(wf, wf->form[PETSC_WF_GP3], label, val, find, part, n3, (void (***)(void))g3)); 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 **g3) in PetscDSGetJacobian() argument
1521 *g3 = tmp3 ? tmp3[0] : NULL; in PetscDSGetJacobian()
1551 … f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2, PetscPointJacFn *g3) in PetscDSSetJacobian() argument
1558 if (g3) PetscValidFunction(g3, 7); 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 **g3) in PetscDSGetJacobianPreconditioner() argument
1659 *g3 = tmp3 ? tmp3[0] : NULL; in PetscDSGetJacobianPreconditioner()
1693 … f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2, PetscPointJacFn *g3) in PetscDSSetJacobianPreconditioner() argument
[all …]
/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 …akFormArrays(rds, NULL, NULL, n0 ? &g0 : NULL, n1 ? &g1 : NULL, n2 ? &g2 : NULL, n3 ? &g3 : NULL)); in PetscFEIntegrateJacobian_Basic()
868 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); in PetscFEIntegrateJacobian_Basic()
869 …ff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, constants, g3); in PetscFEIntegrateJacobian_Basic()
870 for (PetscInt c = 0; c < NcI * NcJ * dE * dE; ++c) g3[c] *= w; in PetscFEIntegrateJacobian_Basic()
873 …eal, basisDerReal, cT[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, ctotDim, offsetI, o… 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()
963 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); in PetscFEIntegrateBdJacobian_Basic()
1028 PetscCall(PetscArrayzero(g3, NcI * NcJ * dE * dE)); in PetscFEIntegrateBdJacobian_Basic()
[all …]
/petsc/src/tao/tutorials/
H A Dex1.c82 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_ul() argument
85 for (d = 0; d < dim; ++d) g3[d * dim + d] = u[1]; in g3_ul()
121 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_lu() argument
124 for (d = 0; d < dim; ++d) g3[d * dim + d] = u[1]; in g3_lu()
H A Dex2.c83 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_ul() argument
86 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_ul()
115 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_lu() argument
118 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_lu()
H A Dex3.c158 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in laplace_kernel() argument
161 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in laplace_kernel()
/petsc/src/snes/tests/
H A Dex13.c38 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
41 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
/petsc/include/petsc/private/
H A Dpetscdsimpl.h143 PetscScalar *g0, *g1, *g2, *g3; /* Point evaluations of weak form Jacobian integrands */ member
/petsc/src/ts/tutorials/
H A Dex76.c955 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_vu() argument
963 g3[((c * Nc + c) * dim + d) * dim + d] += nu; in g3_vu()
964 g3[((c * Nc + d) * dim + d) * dim + c] += nu; in g3_vu()
1023 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_conduct_vu() argument
1033 g3[((c * Nc + c) * dim + d) * dim + d] += mu / Re; // gradU in g3_conduct_vu()
1035 g3[((c * Nc + d) * dim + d) * dim + c] += mu / Re; // gradU transpose in g3_conduct_vu()
1037 g3[((c * Nc + d) * dim + c) * dim + d] -= 2.0 / 3.0 * mu / Re; in g3_conduct_vu()
1065 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_wT() argument
1070 for (d = 0; d < dim; ++d) g3[d * dim + d] = alpha; in g3_wT()
1108 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_conduct_wT() argument
[all …]
H A Dex48.c178 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_nmu() argument
181 for (PetscInt d = 0; d < dim; ++d) g3[d * dim + d] = mu; in g3_nmu()
185 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_n1() argument
188 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1; in g3_n1()
/petsc/src/dm/dt/fe/tests/
H A Dex1.c51 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
54 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
H A Dex2.c42 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g3[]) in g3_uu() argument
45 for (d = 0; d < dim; ++d) g3[d * dim + d] = 1.0; in g3_uu()
/petsc/src/dm/dt/fe/interface/
H A Dfe.c2629 const PetscScalar *G = g3 + (fc * (_NcJ) + gc) * (_dE) * (_dE); \
2645 …ar g0[], const PetscScalar g1[], const PetscScalar g2[], const PetscScalar g3[], PetscInt totDim, … in PetscFEUpdateElementMat_Internal() argument
2761 if (PetscUnlikely(g3)) { in PetscFEUpdateElementMat_Internal()
2785 …Int dg = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * d… in PetscFEUpdateElementMat_Internal()
2800 …ar g0[], const PetscScalar g1[], const PetscScalar g2[], const PetscScalar g3[], PetscInt eOffset,… in PetscFEUpdateElementMat_Hybrid_Internal() argument
2852 …for (dg = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * … in PetscFEUpdateElementMat_Hybrid_Internal()
2899 …st PetscInt g0[], const PetscInt g1[], const PetscInt g2[], const PetscInt g3[], PetscInt totDim, …
2909 PetscBool has_g3 = g3 ? PETSC_TRUE : PETSC_FALSE;
2926 if (g3[i]) g3_i += NbI * NbJ;
2983 if (g3[((fc * NcJ + gc) * dE + df) * dE + dg]) {

12