Home
last modified time | relevance | path

Searched refs:dg (Results 1 – 5 of 5) sorted by relevance

/petsc/src/ts/impls/implicit/discgrad/
H A Dtsdiscgrad.c31 TS_DiscGrad *dg = (TS_DiscGrad *)ts->data; in TSDiscGradGetX0AndXdot() local
40 else *Xdot = dg->Xdot; in TSDiscGradGetX0AndXdot()
108 TS_DiscGrad *dg = (TS_DiscGrad *)ts->data; in TSSetUp_DiscGrad() local
112 if (!dg->X) PetscCall(VecDuplicate(ts->vec_sol, &dg->X)); in TSSetUp_DiscGrad()
113 if (!dg->X0) PetscCall(VecDuplicate(ts->vec_sol, &dg->X0)); in TSSetUp_DiscGrad()
114 if (!dg->Xdot) PetscCall(VecDuplicate(ts->vec_sol, &dg->Xdot)); in TSSetUp_DiscGrad()
124 TS_DiscGrad *dg = (TS_DiscGrad *)ts->data; in TSSetFromOptions_DiscGrad() local
129 …gradient solver", "TSDiscGradSetDGType", DGTypes, (PetscEnum)dg->discgrad, (PetscEnum *)&dg->discg… in TSSetFromOptions_DiscGrad()
147 TS_DiscGrad *dg = (TS_DiscGrad *)ts->data; in TSDiscGradGetType_DiscGrad() local
150 *dgtype = dg->discgrad; in TSDiscGradGetType_DiscGrad()
[all …]
/petsc/src/tao/linesearch/impls/morethuente/
H A Dmorethuente.c38 PetscReal dgx, dgy, dg, dg2, fx, fy, stx, sty, dgtest; in TaoLineSearchApply_MT() local
144 PetscCall(TaoLineSearchComputeObjectiveAndGTS(ls, mt->work, f, &dg)); in TaoLineSearchApply_MT()
150 PetscCall(VecDot(g, x, &dg)); in TaoLineSearchApply_MT()
152 dg = (dg2 - dg) / ls->step; in TaoLineSearchApply_MT()
154 PetscCall(VecDot(g, s, &dg)); in TaoLineSearchApply_MT()
170 if (PetscIsInfOrNanReal(*f) || PetscIsInfOrNanReal(dg)) { in TaoLineSearchApply_MT()
175 dg = PETSC_INFINITY; in TaoLineSearchApply_MT()
182 …if ((*f - ftest1 <= PETSC_SMALL * PetscAbsReal(finit)) && (PetscAbsReal(dg) + ls->gtol * dginit <=… in TaoLineSearchApply_MT()
201 if (ls->step == ls->stepmax && *f <= ftest1 && dg <= dgtest) { in TaoLineSearchApply_MT()
206 if (ls->step == ls->stepmin && *f >= ftest1 && dg >= dgtest) { in TaoLineSearchApply_MT()
[all …]
/petsc/src/dm/dt/fe/interface/
H A Dfe.c2636 … for (PetscInt dg = 0; dg < (_dE); ++dg) s += tBDI[df] * G[df * (_dE) + dg] * tBDJ[dg]; \
2676 … for (PetscInt dg = 0; dg < cdim; ++dg) tmpBasisDerJ[gidx * dE + dg] = basisDerJ[gidx * cdim + dg]; in PetscFEUpdateElementMat_Internal() local
2785 …cInt dg = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * … in PetscFEUpdateElementMat_Internal()
2815 PetscInt f, fc, g, gc, df, dg; in PetscFEUpdateElementMat_Hybrid_Internal() local
2832 for (dg = 0; dg < dE; ++dg) tmpBasisDerJ[gidx * dE + dg] = basisDerJ[gidx * dE + dg]; in PetscFEUpdateElementMat_Hybrid_Internal()
2852 …for (dg = 0; dg < dE; ++dg) elemMat[fOff] += tmpBasisDerI[fidx * dE + df] * g3[((fc * NcJ + gc) * … in PetscFEUpdateElementMat_Hybrid_Internal()
2982 for (PetscInt dg = 0; dg < dE; ++dg) {
2983 if (g3[((fc * NcJ + gc) * dE + df) * dE + dg]) {
2985 g3_idxs[4 * g3_i + 1] = ((fc * NcJ + gc) * dE + df) * dE + dg;
2986 g3_idxs[4 * g3_i + 2] = gidx * dE + dg;
/petsc/src/ts/tutorials/
H A Dex46.c152 PetscInt fc, gc, dg; in g1_uu() local
155 for (dg = 0; dg < dim; ++dg) { in g1_uu()
157 if (fc == gc) g1[(fc * NcJ + gc) * dim + dg] += u[dg]; in g1_uu()
H A Dex48.c192 static PetscScalar poissonBracket(PetscInt dim, const PetscScalar df[], const PetscScalar dg[]) in poissonBracket() argument
194 PetscScalar ret = df[0] * dg[1] - df[1] * dg[0]; in poissonBracket()
196 ret += df[1] * dg[2] - df[2] * dg[1]; in poissonBracket()
197 ret += df[2] * dg[0] - df[0] * dg[2]; in poissonBracket()