Home
last modified time | relevance | path

Searched refs:g1 (Results 1 – 25 of 32) sorted by relevance

12

/petsc/src/dm/dt/interface/
H A Ddtweakform.c556 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n1, void (***g1)(PetscInt, PetscInt,… in PetscWeakFormGetJacobian()
562 …mGetFunction_Private(wf, wf->form[PETSC_WF_G1], label, val, find, part, n1, (void (***)(void))g1)); in PetscWeakFormGetJacobian()
568 …, const PetscReal[], PetscInt, const PetscScalar[], PetscScalar[]), void (*g1)(PetscInt, PetscInt,… in PetscWeakFormAddJacobian()
574 …WeakFormAddFunction_Private(wf, wf->form[PETSC_WF_G1], label, val, find, part, (PetscVoidFn *)g1)); in PetscWeakFormAddJacobian()
580 …eal[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt n1, void (**g1)(PetscInt, PetscInt,… in PetscWeakFormSetJacobian()
586 …ormSetFunction_Private(wf, wf->form[PETSC_WF_G1], label, val, find, part, n1, (PetscVoidFn **)g1)); in PetscWeakFormSetJacobian()
592 …Real[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt i1, void (*g1)(PetscInt, PetscInt,… in PetscWeakFormSetIndexJacobian()
598 …etIndexFunction_Private(wf, wf->form[PETSC_WF_G1], label, val, find, part, i1, (PetscVoidFn *)g1)); in PetscWeakFormSetIndexJacobian()
619 …l[], PetscInt, const PetscScalar[], PetscScalar[]), PetscInt *n1, void (***g1)(PetscInt, PetscInt,… in PetscWeakFormGetJacobianPreconditioner()
625 …GetFunction_Private(wf, wf->form[PETSC_WF_GP1], label, val, find, part, n1, (void (***)(void))g1)); in PetscWeakFormGetJacobianPreconditioner()
[all …]
H A Ddtds.c518 …NcMax, &prob->g0, NsMax * NsMax * NqMax * NcMax * NcMax * dimEmbed, &prob->g1, NsMax * NsMax * NqM… in PetscDSSetUp()
534 PetscCall(PetscFree6(prob->f0, prob->f1, prob->g0, prob->g1, prob->g2, prob->g3)); in PetscDSDestroyStructs_Static()
1508 …etscDS ds, PetscInt f, PetscInt g, PetscPointJacFn **g0, PetscPointJacFn **g1, PetscPointJacFn **g… in PetscDSGetJacobian() argument
1519 *g1 = tmp1 ? tmp1[0] : NULL; in PetscDSGetJacobian()
1551 …(PetscDS ds, PetscInt f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2… in PetscDSSetJacobian() argument
1556 if (g1) PetscValidFunction(g1, 5); in PetscDSSetJacobian()
1561 PetscCall(PetscWeakFormSetIndexJacobian(ds->wf, NULL, 0, f, g, 0, 0, g0, 0, g1, 0, g2, 0, g3)); in PetscDSSetJacobian()
1646 …etscDS ds, PetscInt f, PetscInt g, PetscPointJacFn **g0, PetscPointJacFn **g1, PetscPointJacFn **g… in PetscDSGetJacobianPreconditioner() argument
1657 *g1 = tmp1 ? tmp1[0] : NULL; in PetscDSGetJacobianPreconditioner()
1693 …(PetscDS ds, PetscInt f, PetscInt g, PetscPointJacFn *g0, PetscPointJacFn *g1, PetscPointJacFn *g2… in PetscDSSetJacobianPreconditioner() argument
[all …]
/petsc/src/tao/interface/
H A Dtaosolver_fg.c27 PetscErrorCode TaoTestGradient(Tao tao, Vec x, Vec g1) in TaoTestGradient() argument
72 PetscCall(VecNorm(g1, NORM_2, &hcnorm)); in TaoTestGradient()
74 PetscCall(VecNorm(g1, NORM_INFINITY, &hcmax)); in TaoTestGradient()
75 PetscCall(VecDot(g1, g2, &dot)); in TaoTestGradient()
76 PetscCall(VecCopy(g1, g3)); in TaoTestGradient()
86 PetscCall(VecView(g1, mviewer)); in TaoTestGradient()
/petsc/src/benchmarks/streams/
H A Dmakefile60 ${MPIEXEC} ${MPI_BINDING} -n 1 -c$${i} -a$${i} -g1 ./CUDAVersion | tee -a scaling.log; \
64 ${MPIEXEC} ${MPI_BINDING} -n $${n} -c$${c} -a$${c} -g1 ./CUDAVersion | tee -a scaling.log; \
/petsc/src/ts/utils/dmplexlandau/
H A Dland_tensors.h131 …const PetscReal x2[3] = {xp, yp, zp}, x1[3] = {a_x1[0], a_x1[1], a_x1[2]}, c02 = c0 * c0, g1 = GAM… 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()
858 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateJacobian_Basic()
859 …ff_x, u, u_t, u_x, aOff, aOff_x, a, NULL, a_x, t, u_tshift, fegeom.v, numConstants, constants, g1); in PetscFEIntegrateJacobian_Basic()
860 for (PetscInt c = 0; c < NcI * NcJ * dE; ++c) g1[c] *= w; in PetscFEIntegrateJacobian_Basic()
873 …, basisReal, basisDerReal, cT[fieldJ], testReal, testDerReal, &fegeom, g0, g1, g2, g3, ctotDim, of… 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()
961 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
1018 PetscCall(PetscArrayzero(g1, NcI * NcJ * dE)); in PetscFEIntegrateBdJacobian_Basic()
[all …]
/petsc/src/ts/tutorials/
H A Dex48.c126 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_phi_right() argument
136 g1[i] += s_K[i][j] * pphiDer[j]; in g1_phi_right()
140 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_njz_left() argument
148 g1[j] += -jzDer[i] * s_K[i][j]; in g1_njz_left()
151 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_npsi_right() argument
156 for (j = 0; j < dim; ++j) g1[i] += -s_K[i][j] * psiDer[j]; in g1_npsi_right()
160 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_omega_left() argument
165 for (j = 0; j < dim; ++j) g1[j] += pOmegaDer[i] * s_K[i][j]; in g1_omega_left()
169 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_psi_left() argument
174 for (j = 0; j < dim; ++j) g1[j] += pPsiDer[i] * s_K[i][j]; in g1_psi_left()
H A Dex46.c148 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_uu() argument
157 if (fc == gc) g1[(fc * NcJ + gc) * dim + dg] += u[dg]; in g1_uu()
165 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_pu() argument
168 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; /* \frac{\partial\phi^{u_d}}{\partial x_d} */ in g1_pu()
H A Dex76.c875 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_qu() argument
878 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_qu()
893 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_vu() argument
902 if (c == d) g1[(c * NcJ + d) * dim + e] += u[e]; in g1_vu()
917 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_conduct_qu() argument
923 for (d = 0; d < dim; ++d) g1[d * dim + d] = p_th / u[uOff[TEMP]]; in g1_conduct_qu()
940 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_conduct_qT() argument
946 for (d = 0; d < dim; ++d) g1[d] = -p_th / PetscSqr(u[uOff[TEMP]]) * u[uOff[VEL] + d]; in g1_conduct_qT()
1006 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_conduct_vu() argument
1017 if (c == d) g1[(c * NcJ + d) * dim + e] += p_th / u[uOff[TEMP]] * u[uOff[VEL] + e]; in g1_conduct_vu()
[all …]
H A Dex77.c211 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_qu() argument
214 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_qu()
229 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_vu() argument
238 if (c == d) g1[(c * NcJ + d) * dim + e] += u[e]; in g1_vu()
276 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_wT() argument
279 for (d = 0; d < dim; ++d) g1[d] = u[uOff[0] + d]; in g1_wT()
H A Dex47.c103 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_prim_phi() argument
106 for (d = 0; d < dim; ++d) g1[d] = a[d]; in g1_prim_phi()
H A Dex18.c268 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_adv_pp() argument
271 for (d = 0; d < dim; ++d) g1[d] = u[d]; in g1_adv_pp()
280 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_adv_pu() argument
283 for (d = 0; d < dim; ++d) g1[d * dim + d] = u[dim]; in g1_adv_pu()
/petsc/src/tao/pde_constrained/impls/lcl/
H A Dlcl.c44 PetscCall(VecDestroy(&lclP->g1)); in TaoDestroy_LCL()
151 PetscCall(VecDuplicate(lclP->V, &lclP->g1)); in TaoSetup_LCL()
427 PetscCall(MatMultTranspose(tao->jacobian_design, lclP->lambda, lclP->g1)); in TaoSolve_LCL()
428 PetscCall(VecAXPY(lclP->g1, -1.0, lclP->GAugL_V)); in TaoSolve_LCL()
432 PetscCall(MatSolve(lclP->R, lclP->g1, lclP->s)); in TaoSolve_LCL()
433 PetscCall(VecDot(lclP->s, lclP->g1, &descent)); in TaoSolve_LCL()
436 PetscCall(VecCopy(lclP->g1, lclP->s)); in TaoSolve_LCL()
439 PetscCall(VecCopy(lclP->g1, lclP->s)); in TaoSolve_LCL()
441 PetscCall(VecScale(lclP->g1, -1.0)); in TaoSolve_LCL()
H A Dlcl.h52 Vec g1, g2; member
/petsc/src/snes/tutorials/
H A Dex76.c212 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_qu() argument
215 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_qu()
228 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_vu() argument
237 if (c == d) g1[(c * NcJ + d) * dim + e] = u[e]; in g1_vu()
269 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_wT() argument
272 for (d = 0; d < dim; ++d) g1[d] = u[uOff[0] + d]; in g1_wT()
H A Dex77.c188 …l x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_bd_uu_3d() argument
203 …for (d = 0; d < dim; ++d) g1[(compI * Ncomp + compJ) * dim + d] = p * (m[compI] * cofu_x[compJ * d… in g1_bd_uu_3d()
215 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_pu_3d() argument
222 for (d = 0; d < dim; ++d) g1[compI * dim + d] = cofu_x[compI * dim + d]; in g1_pu_3d()
H A Dex71.c140 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_pu() argument
143 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; /* \frac{\partial\phi^{u_d}}{\partial x_d} */ in g1_pu()
H A Dex62.c71 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_pu() argument
74 for (d = 0; d < dim; ++d) g1[d * dim + d] = -1.0; /* < q, -\nabla\cdot u > */ in g1_pu()
H A Dex24.c198 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_uq() argument
201 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_uq()
/petsc/src/dm/impls/plex/tests/
H A Dex39.c153 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_qu() argument
156 for (d = 0; d < dim; ++d) g1[d * dim + d] = -1.0; in g1_qu()
175 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_wu() argument
178 for (d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_wu()
/petsc/include/petsc/private/
H A Dpetscdsimpl.h143 PetscScalar *g0, *g1, *g2, *g3; /* Point evaluations of weak form Jacobian integrands */ member
/petsc/src/dm/dt/tests/
H A Dex10.c89 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_qu() argument
93 for (c = 0; c < dim; ++c) g1[c * dim + c] = 1.0; in g1_qu()
/petsc/src/dm/dt/fe/interface/
H A Dfe.c2593 const PetscScalar *G = g1 + (fc * (_NcJ) + gc) * _dE; \
2645 …asisDerJ[], PetscFEGeom *fegeom, const PetscScalar g0[], const PetscScalar g1[], const PetscScalar… in PetscFEUpdateElementMat_Internal() argument
2701 if (PetscUnlikely(g1)) { in PetscFEUpdateElementMat_Internal()
2724 …for (PetscInt df = 0; df < dE; ++df) elemMat[fOff] += bI * g1[(fc * NcJ + gc) * dE + df] * tmpBasi… in PetscFEUpdateElementMat_Internal()
2800 …asisDerJ[], PetscFEGeom *fegeom, const PetscScalar g0[], const PetscScalar g1[], const PetscScalar… in PetscFEUpdateElementMat_Hybrid_Internal() argument
2850 … elemMat[fOff] += tmpBasisI[fidx] * g1[(fc * NcJ + gc) * dE + df] * tmpBasisDerJ[gidx * dE + df]; in PetscFEUpdateElementMat_Hybrid_Internal()
2899 …PetscTabulation TJ, PetscInt dimEmbed, const PetscInt g0[], const PetscInt g1[], const PetscInt g2…
2907 PetscBool has_g1 = g1 ? PETSC_TRUE : PETSC_FALSE;
2920 if (g1[i]) g1_i += NbI * NbJ;
2964 if (g1[(fc * NcJ + gc) * dE + df]) {
/petsc/src/snes/tests/
H A Dex15.c180 …Shift, const PetscReal x[], PetscInt numConstants, const PetscScalar constants[], PetscScalar g1[]) in g1_phiq() argument
182 for (PetscInt d = 0; d < dim; ++d) g1[d * dim + d] = 1.0; in g1_phiq()
/petsc/src/snes/interface/
H A Dsnes.c2661 Vec x, g1, g2, g3; in SNESTestFunction() local
2697 PetscCall(VecDuplicate(x, &g1)); in SNESTestFunction()
2700 …PetscCall(SNESComputeFunction(snes, x, g1)); /* does not handle use of SNESSetFunctionDomainError(… in SNESTestFunction()
2704 PetscCall(VecNorm(g1, NORM_2, &hcnorm)); in SNESTestFunction()
2706 PetscCall(VecNorm(g1, NORM_INFINITY, &hcmax)); in SNESTestFunction()
2707 PetscCall(VecDot(g1, g2, &dot)); in SNESTestFunction()
2708 PetscCall(VecCopy(g1, g3)); in SNESTestFunction()
2718 PetscCall(VecView(g1, mviewer)); in SNESTestFunction()
2724 PetscCall(VecDestroy(&g1)); in SNESTestFunction()

12