Home
last modified time | relevance | path

Searched refs:d2 (Results 1 – 25 of 31) sorted by relevance

12

/petsc/src/tao/unconstrained/tutorials/
H A Dminsurf1.c131 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormFunctionGradient() local
181 d2 = (xc - xr); in FormFunctionGradient()
192 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient()
193 df5dxc = d2 * hydhx; in FormFunctionGradient()
197 d2 *= rhx; in FormFunctionGradient()
208 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient()
209 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormFunctionGradient()
227 d2 = (user->left[j + 1] - x[j * mx]) * rhx; in FormFunctionGradient()
228 ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient()
232 d2 = (user->bottom[i + 1] - user->bottom[i + 2]) * rhx; in FormFunctionGradient()
[all …]
H A Dminsurf2.c170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local
217 d2 = (xc - xr); in FormFunction()
222 d2 *= rhx; in FormFunction()
227 f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
237 d2 = (user->left[j - ys + 1] - x[j][0]) * rhx; in FormFunction()
238 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
243 d2 = (user->bottom[i + 1 - xs] - user->bottom[i - xs + 2]) * rhx; in FormFunction()
245 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction()
264 d2 = (user->bottom[0] - user->bottom[1]) * rhx; in FormFunction()
265 ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2); in FormFunction()
[all …]
/petsc/src/snes/tutorials/
H A Dex4.c119 PetscScalar f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormObjectiveLocal() local
149 d2 = (xc - xr); in FormObjectiveLocal()
154 d2 /= hx; in FormObjectiveLocal()
159 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormObjectiveLocal()
169 d2 = (user->left[j + 1] - x[j][0]) / hx; in FormObjectiveLocal()
170 ft += PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormObjectiveLocal()
175 d2 = (user->bottom[i + 1] - user->bottom[i + 2]) / hx; in FormObjectiveLocal()
177 ft += PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormObjectiveLocal()
196 d2 = (user->bottom[0] - user->bottom[1]) / hx; in FormObjectiveLocal()
197 ft += PetscSqrtReal(1.0 + d1 * d1 + d2 * d2); in FormObjectiveLocal()
[all …]
H A Dex58.c149 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormGradient() local
207 d2 = (xc - xr); in FormGradient()
218 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormGradient()
219 df5dxc = d2 * hydhx; in FormGradient()
223 d2 /= hx; in FormGradient()
234 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormGradient()
235 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormGradient()
277 …PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormJacobian() local
345 d2 = (xc - xr) / hx; in FormJacobian()
356 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormJacobian()
[all …]
H A Dex59.c130 PetscScalar *ff, *FF, d, d2; in FormFunction() local
139 d2 = d * d; in FormFunction()
144 …for (i = 1; i < n - 1; i++) ff[i] = d2 * (xx[i - 1] - 2. * xx[i] + xx[i + 1]) + xx[i] * xx[i] - FF… in FormFunction()
156 PetscScalar A[3], d, d2; in FormJacobian() local
163 d2 = d * d; in FormJacobian()
185 A[0] = d2; in FormJacobian()
186 A[1] = -2. * d2 + 2. * xx[i]; in FormJacobian()
187 A[2] = d2; in FormJacobian()
H A Dex11.c71 for (PetscInt d2 = 0; d2 < dim; ++d2) v += K[d1 * dim + d2] * u_x[d2]; in f_1() local
82 for (PetscInt d2 = 0; d2 < dim; ++d2) J11[d1 * dim + d2] = K[d1 * dim + d2]; in J_11() local
H A Dex77.c150 PetscInt compI, compJ, d1, d2; in g3_uu_3d() local
164 for (d2 = 0; d2 < dim; ++d2) { in g3_uu_3d()
165 const PetscReal g = (d1 == d2) ? 1.0 : 0.0; in g3_uu_3d()
167 …J) * dim + d1) * dim + d2] = g * G * mu + pp * cofu_x[compI * dim + d1] * cofu_x[compJ * dim + d2]… in g3_uu_3d()
/petsc/src/dm/ftn-mod/
H A Dpetscdmcomposite.h9032 subroutine DMCompositeGetEntries2(dm, d1, d2, ierr)
36 DM d1, d2
39 subroutine DMCompositeGetEntries3(dm, d1, d2, d3, ierr)
43 DM d1, d2, d3
46 subroutine DMCompositeGetEntries4(dm, d1, d2, d3, d4, ierr)
50 DM d1, d2, d4, d3
/petsc/src/tao/complementarity/tutorials/
H A Dminsurf1.c142 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormConstraints() local
195 d2 = (xc - xr); in FormConstraints()
206 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormConstraints()
207 df5dxc = d2 * hydhx; in FormConstraints()
211 d2 /= hx; in FormConstraints()
222 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormConstraints()
223 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormConstraints()
264 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormJacobian() local
320 d2 = (xc - xr) / hx; in FormJacobian()
331 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormJacobian()
[all …]
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F9072 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3
160 d2 = xc - xr
171 df4dxc = d2*hydhx + d3*hxdhy
172 df5dxc = d2*hydhx
176 d2 = d2*rhx
187 f4 = sqrt(1.0 + d3*d3 + d2*d2)
188 f5 = sqrt(1.0 + d2*d2 + d5*d5)
208 d2 = (left_v(1 + j - ys + 1) - x_v(1 + (j - gys)*gxm))*rhx
209 ft = ft + sqrt(1.0 + d3*d3 + d2*d2)
215 d2 = (bottom_v(1 + i + 1 - xs) - bottom_v(1 + i - xs + 2))*rhx
[all …]
H A Dplate2.c207 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionGradient() local
272 d2 = (xc - xr); in FormFunctionGradient()
283 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient()
284 df5dxc = d2 * hydhx; in FormFunctionGradient()
288 d2 *= rhx; in FormFunctionGradient()
299 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient()
300 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormFunctionGradient()
320 d2 = (left[j - ys + 1] - x[(j - gys) * gxm]) * rhx; in FormFunctionGradient()
321 ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient()
326 d2 = (bottom[i + 1 - xs] - bottom[i - xs + 2]) * rhx; in FormFunctionGradient()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/pipebcgs/
H A Dpipebcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega = 0.0, d1, d2, d3; in KSPSolve_PIPEBCGS() local
81 PetscCall(VecDotBegin(W, RP, &d2)); /* d2 <- (w,rp) */ in KSPSolve_PIPEBCGS()
85 PetscCall(VecDotEnd(W, RP, &d2)); in KSPSolve_PIPEBCGS()
87 alpha = rho / d2; in KSPSolve_PIPEBCGS()
109 PetscCall(VecDotBegin(Y, Y, &d2)); /* d2 <- (y,y) */ in KSPSolve_PIPEBCGS()
116 PetscCall(VecDotEnd(Y, Y, &d2)); in KSPSolve_PIPEBCGS()
118 if (d2 == 0.0) { in KSPSolve_PIPEBCGS()
135 omega = d1 / d2; /* omega <- (y'q) / (y'y) */ in KSPSolve_PIPEBCGS()
147 …PetscCall(VecDotBegin(W, RP, &d2)); … in KSPSolve_PIPEBCGS()
157 PetscCall(VecDotEnd(W, RP, &d2)); in KSPSolve_PIPEBCGS()
[all …]
/petsc/src/ts/utils/dmplexlandau/kokkos/
H A Dlandau.kokkos.cxx614 PetscInt fieldB, d2, d3, f_off_r, grid_r, ipidx_g, nip_loc_r, loc_Nf_r; in LandauKokkosJacobian() local
645 for (d2 = 0; d2 < 2; d2++) { in LandauKokkosJacobian()
648 ggg.gg2[d2] += Uk[d2][d3] * temp1[d3]; in LandauKokkosJacobian()
650 ggg.gg3[d2][d3] += Ud[d2][d3] * temp2; in LandauKokkosJacobian()
654 for (d2 = 0; d2 < 3; ++d2) { in LandauKokkosJacobian()
657 ggg.gg2[d2] += U[d2][d3]*temp1[d3]; in LandauKokkosJacobian()
659 ggg.gg3[d2][d3] += U[d2][d3]*temp2; in LandauKokkosJacobian()
667 PetscInt d2, d3; in LandauKokkosJacobian() local
668 for (d2 = 0; d2 < dim; d2++) { in LandauKokkosJacobian()
669 gg2(d2, fieldA, myQi) = gg_temp.gg2[d2] * d_alpha[fieldA + f_off]; in LandauKokkosJacobian()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c20 PetscReal dp = 0.0, d2; in KSPSolve_FBCGS() local
97 PetscCall(VecDotNorm2(S, T, &d1, &d2)); in KSPSolve_FBCGS()
98 if (d2 == 0.0) { in KSPSolve_FBCGS()
117 omega = d1 / d2; /* omega <- (t's) / (t't) */ in KSPSolve_FBCGS()
/petsc/src/tao/matrix/
H A Dadamat.c352 static PetscErrorCode MatCreateADA(Mat mat, Vec d1, Vec d2, Mat *J) in MatCreateADA() argument
362 ctx->D2 = d2; in MatCreateADA()
369 if (d2) { in MatCreateADA()
370 PetscCall(VecDuplicate(d2, &ctx->W2)); in MatCreateADA()
371 PetscCall(VecDuplicate(d2, &ctx->ADADiag)); in MatCreateADA()
372 PetscCall(PetscObjectReference((PetscObject)d2)); in MatCreateADA()
381 PetscCall(VecGetLocalSize(d2, &nloc)); in MatCreateADA()
382 PetscCall(VecGetSize(d2, &n)); in MatCreateADA()
/petsc/src/ksp/ksp/impls/bcgs/
H A Dbcgs.c23 PetscReal dp = 0.0, d2; in KSPSolve_BCGS() local
89 PetscCall(VecDotNorm2(S, T, &d1, &d2)); in KSPSolve_BCGS()
90 if (d2 == 0.0) { in KSPSolve_BCGS()
110 omega = d1 / d2; /* w <- (t's) / (t't) */ in KSPSolve_BCGS()
/petsc/src/dm/dt/tests/
H A Dex1.c19 PetscReal b, d, d2; in CheckPoints() local
22 d2 = D2[i * ndegrees + j]; in CheckPoints()
25 if (PetscAbsReal(d2) < PETSC_SMALL) d2 = 0; in CheckPoints()
26 …=%12.4g D=%12.4g D2=%12.4g\n", degrees[j], (double)points[i], (double)b, (double)d, (double)d2)); in CheckPoints()
/petsc/src/dm/dt/space/impls/poly/
H A Dspacepoly.c210 for (PetscInt d2 = 0; d2 < dim; d2++) { in PetscSpaceEvaluate_Polynomial() local
213 derivs[d2]++; in PetscSpaceEvaluate_Polynomial()
216 derivs[d2]--; in PetscSpaceEvaluate_Polynomial()
218 …for (PetscInt p = 0; p < npoints; p++) H[p * p_strl + b * b_strl + d1 * d1_strl + d2 * d2_strl] = … in PetscSpaceEvaluate_Polynomial()
/petsc/src/dm/dt/space/impls/ptrimmed/
H A Dspaceptrimmed.c216 for (PetscInt d2 = 0; d2 < dim; d2++) { in PetscSpaceEvaluate_Ptrimmed() local
219 derivs[d2]++; in PetscSpaceEvaluate_Ptrimmed()
222 derivs[d2]--; in PetscSpaceEvaluate_Ptrimmed()
225 … p < npoints; p++) H[p * p_strl + b * b_strl + v * v_strl + d1 * d1_strl + d2 * d2_strl] = eval[b … in PetscSpaceEvaluate_Ptrimmed()
/petsc/src/ts/utils/dmplexlandau/
H A Dplexland.c337 PetscInt d, d2, dp, d3, IPf_idx; in LandauFormJacobian_Internal() local
345 for (d2 = 0; d2 < LANDAU_DIM; d2++) gg3_temp[d][d2] = 0; in LandauFormJacobian_Internal()
383 for (d2 = 0; d2 < 2; d2++) { in LandauFormJacobian_Internal()
386 gg2_temp[d2] += Uk[d2][d3] * temp1[d3]; in LandauFormJacobian_Internal()
388 gg3_temp[d2][d3] += Ud[d2][d3] * temp2; in LandauFormJacobian_Internal()
392 for (d2 = 0; d2 < 3; ++d2) { in LandauFormJacobian_Internal()
395 gg2_temp[d2] += U[d2][d3] * temp1[d3]; in LandauFormJacobian_Internal()
397 gg3_temp[d2][d3] += U[d2][d3] * temp2; in LandauFormJacobian_Internal()
408 for (d2 = 0; d2 < LANDAU_DIM; d2++) { in LandauFormJacobian_Internal()
409 gg2[fieldA][d2] = gg2_temp[d2] * nu_alpha[fieldA + f_off]; in LandauFormJacobian_Internal()
[all …]
/petsc/src/mat/impls/sbaij/seq/
H A Dsbaijfact10.c16 MatScalar d0, d1, d2, d3, d4, d5, d6, d7, d8, d9, d10, d11, d12; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
105 d2 = d[2]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
142 m2 = uik[2] = -(d2 * u0 + d8 * u1 + d14 * u2 + d20 * u3 + d26 * u4 + d32 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
149 m8 = uik[8] = -(d2 * u6 + d8 * u7 + d14 * u8 + d20 * u9 + d26 * u10 + d32 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
156 m14 = uik[14] = -(d2 * u12 + d8 * u13 + d14 * u14 + d20 * u15 + d26 * u16 + d32 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
163 m20 = uik[20] = -(d2 * u18 + d8 * u19 + d14 * u20 + d20 * u21 + d26 * u22 + d32 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
170 m26 = uik[26] = -(d2 * u24 + d8 * u25 + d14 * u26 + d20 * u27 + d26 * u28 + d32 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
177 m32 = uik[32] = -(d2 * u30 + d8 * u31 + d14 * u32 + d20 * u33 + d26 * u34 + d32 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
/petsc/src/tao/quadratic/impls/bqpip/
H A Dbqpip.c241 PetscReal d1, d2, ksptol, sigmamu; in TaoSolve_BQPIP() local
435 PetscCall(VecDot(tao->solution, qp->C, &d2)); in TaoSolve_BQPIP()
440 qp->pobj = d1 / 2.0 + d2 + qp->d; in TaoSolve_BQPIP()
/petsc/src/dm/dt/interface/
H A Ddt.c2431 PetscReal lterm, rterm, maxTerm = 0.0, d1, d2, d3, d4; in PetscDTTanhSinhIntegrate() local
2460 d2 = PetscLog10Real(PetscAbsReal(sum - psum)); in PetscDTTanhSinhIntegrate()
2464 d = PetscAbsInt(PetscMin(0, PetscMax(PetscMax(PetscMax(PetscSqr(d1) / d2, 2 * d1), d3), d4))); in PetscDTTanhSinhIntegrate()
2509 PetscReal d1, d2, d3, d4; in PetscDTTanhSinhIntegrateMPFR() local
2576 d2 = mpfr_get_d(tmp, MPFR_RNDN); in PetscDTTanhSinhIntegrateMPFR()
2581 d = PetscAbsInt(PetscMin(0, PetscMax(PetscMax(PetscMax(PetscSqr(d1) / d2, 2 * d1), d3), d4))); in PetscDTTanhSinhIntegrateMPFR()
2618 PetscInt dim2, Nc2, Np2, order2, qb, d2; in PetscDTTensorQuadratureCreate() local
2774 for (d2 = 0; d2 < dim2; ++d2, ++d) x[qc * dim + d] = x2[qb * dim2 + d2]; in PetscDTTensorQuadratureCreate()
/petsc/lib/petsc/bin/maint/
H A Dtoclapack.sh2316 static doublereal b1, b2, c1, c2, d1, d2, one, zero, rbase;
2364 d2 = a;
2368 if (c1 == a && c2 == a && d1 == a && d2 == a) {
2385 d2 = zero;
2388 d2 += b2;
4307 static halfreal b1, b2, c1, c2, d1, d2, one, zero, rbase;
4355 d2 = a;
4359 if (c1 == a && c2 == a && d1 == a && d2 == a) {
4376 d2 = zero;
4379 d2 += b2;
/petsc/src/tao/pde_constrained/tutorials/
H A Delliptic.c185 PetscReal d1 = 0, d2 = 0; in FormFunction() local
195 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunction()
196 *f = 0.5 * (d1 + user->alpha * d2); in FormFunction()
224 PetscReal d1, d2; in FormFunctionGradient() local
236 PetscCall(VecDot(user->lwork, user->lwork, &d2)); in FormFunctionGradient()
239 *f = 0.5 * (d1 + user->alpha * d2); in FormFunctionGradient()

12