Home
last modified time | relevance | path

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

12

/petsc/src/dm/ftn-mod/
H A Dpetscdmcomposite.h9025 subroutine DMCompositeGetEntries1(dm, d1, ierr)
29 DM d1
32 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/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
180 d1 = (xc - xl); in FormFunctionGradient()
189 df1dxc = d1 * hydhx; in FormFunctionGradient()
190 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient()
196 d1 *= rhx; in FormFunctionGradient()
205 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormFunctionGradient()
206 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient()
238 d1 = (x[(j + 1) * mx - 1] - user->right[j + 1]) * rhx; in FormFunctionGradient()
240 ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient()
244 d1 = (x[(my - 1) * mx + i] - user->top[i + 1]) * rhy; in FormFunctionGradient()
[all …]
H A Dminsurf2.c170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local
216 d1 = (xc - xl); in FormFunction()
221 d1 *= rhx; in FormFunction()
226 f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
250 d1 = (x[j][mx - 1] - user->right[j - ys + 1]) * rhx; in FormFunction()
252 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
257 d1 = (x[my - 1][i] - user->top[i - xs + 1]) * rhy; in FormFunction()
259 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction()
263 d1 = (user->left[0] - user->left[1]) / hy; 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
148 d1 = (xc - xl); in FormObjectiveLocal()
153 d1 /= hx; in FormObjectiveLocal()
158 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal()
182 d1 = (x[j][mx - 1] - user->right[j + 1]) / hx; in FormObjectiveLocal()
184 ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal()
189 d1 = (x[my - 1][i] - user->top[i + 1]) / hy; in FormObjectiveLocal()
191 ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal()
195 d1 = (user->left[0] - user->left[1]) / hy; 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
206 d1 = (xc - xl); in FormGradient()
215 df1dxc = d1 * hydhx; in FormGradient()
216 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormGradient()
222 d1 /= hx; in FormGradient()
231 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormGradient()
232 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); 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
344 d1 = (xc - xl) / hx; in FormJacobian()
353 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormJacobian()
[all …]
H A Dex11.c69 for (PetscInt d1 = 0; d1 < dim; ++d1) { in f_1() local
71 for (PetscInt d2 = 0; d2 < dim; ++d2) v += K[d1 * dim + d2] * u_x[d2]; in f_1()
72 f1[d1] = v; in f_1()
81 for (PetscInt d1 = 0; d1 < dim; ++d1) { in J_11() local
82 for (PetscInt d2 = 0; d2 < dim; ++d2) J11[d1 * dim + d2] = K[d1 * dim + d2]; in J_11()
H A Dex77.c150 PetscInt compI, compJ, d1, d2; in g3_uu_3d() local
163 for (d1 = 0; d1 < dim; ++d1) { in g3_uu_3d()
165 const PetscReal g = (d1 == d2) ? 1.0 : 0.0; in g3_uu_3d()
167 …m + d1) * dim + d2] = g * G * mu + pp * cofu_x[compI * dim + d1] * cofu_x[compJ * dim + d2] - pm *… in g3_uu_3d()
/petsc/src/sys/fileio/
H A Dghome.c28 const char *d1; in PetscGetHomeDirectory() local
31 d1 = getenv("HOME"); in PetscGetHomeDirectory()
32 if (d1) PetscCall(PetscStrncpy(dir, d1, maxlen)); in PetscGetHomeDirectory()
/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
194 d1 = (xc - xl); in FormConstraints()
203 df1dxc = d1 * hydhx; in FormConstraints()
204 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormConstraints()
210 d1 /= hx; in FormConstraints()
219 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormConstraints()
220 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); 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
319 d1 = (xc - xl) / hx; in FormJacobian()
328 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormJacobian()
[all …]
/petsc/src/sys/ftn-mod/
H A Dpetscsys.h9075 subroutine PetscShmgetAllocateArrayScalar(start, len, d1, ierr)
78 PetscScalar, pointer :: d1(:)
82 subroutine PetscShmgetDeallocateArrayScalar(d1, ierr)
84 PetscScalar, pointer :: d1(:)
88 subroutine PetscShmgetAllocateArrayInt(start, len, d1, ierr)
91 PetscInt, pointer :: d1(:)
95 subroutine PetscShmgetDeallocateArrayInt(d1, ierr)
97 PetscInt, pointer :: d1(:)
/petsc/src/tao/bound/tutorials/
H A Dplate2f.F9072 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3
159 d1 = xc - xl
168 df1dxc = d1*hydhx
169 df2dxc = d1*hydhx + d4*hxdhy
175 d1 = d1*rhx
184 f1 = sqrt(1.0 + d1*d1 + d7*d7)
185 f2 = sqrt(1.0 + d1*d1 + d4*d4)
223 d1 = (x_v(1 + (j + 1 - gys)*gxm - 1) - right_v(1 + j - ys + 1))*rhx
225 ft = ft + sqrt(1.0 + d1*d1 + d4*d4)
231 d1 = (x_v(1 + (gym - 1)*gxm + i - gxs) - top_v(1 + i - xs + 1))*rhy
[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
271 d1 = (xc - xl); in FormFunctionGradient()
280 df1dxc = d1 * hydhx; in FormFunctionGradient()
281 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient()
287 d1 *= rhx; in FormFunctionGradient()
296 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormFunctionGradient()
297 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient()
334 d1 = (x[(j + 1 - gys) * gxm - 1] - right[j - ys + 1]) * rhx; in FormFunctionGradient()
336 ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient()
341 d1 = (x[(gym - 1) * gxm + i - gxs] - top[i - xs + 1]) * rhy; in FormFunctionGradient()
[all …]
/petsc/src/ksp/ksp/impls/bcgs/fbcgs/
H A Dfbcgs.c18 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_FBCGS() local
84 PetscCall(VecDot(V, RP, &d1)); in KSPSolve_FBCGS()
85 if (d1 == 0.0) { in KSPSolve_FBCGS()
91 alpha = rho / d1; /* alpha <- rho / (v,rp) */ in KSPSolve_FBCGS()
97 PetscCall(VecDotNorm2(S, T, &d1, &d2)); in KSPSolve_FBCGS()
100 PetscCall(VecDot(S, S, &d1)); in KSPSolve_FBCGS()
101 if (d1 != 0.0) { in KSPSolve_FBCGS()
117 omega = d1 / d2; /* omega <- (t's) / (t't) */ in KSPSolve_FBCGS()
/petsc/src/ksp/ksp/impls/bcgs/
H A Dbcgs.c21 PetscScalar rho, rhoold, alpha, beta, omega, omegaold, d1; in KSPSolve_BCGS() local
78 PetscCall(VecDot(V, RP, &d1)); in KSPSolve_BCGS()
79 KSPCheckDot(ksp, d1); in KSPSolve_BCGS()
80 if (d1 == 0.0) { in KSPSolve_BCGS()
86 alpha = rho / d1; /* a <- rho / (v,rp) */ in KSPSolve_BCGS()
89 PetscCall(VecDotNorm2(S, T, &d1, &d2)); in KSPSolve_BCGS()
93 PetscCall(VecDot(S, S, &d1)); in KSPSolve_BCGS()
94 if (d1 != 0.0) { in KSPSolve_BCGS()
110 omega = d1 / d2; /* w <- (t's) / (t't) */ in KSPSolve_BCGS()
/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
108 PetscCall(VecDotBegin(Q, Y, &d1)); /* d1 <- (q,y) */ in KSPSolve_PIPEBCGS()
115 PetscCall(VecDotEnd(Q, Y, &d1)); in KSPSolve_PIPEBCGS()
120 PetscCall(VecDot(Q, Q, &d1)); in KSPSolve_PIPEBCGS()
121 if (d1 != 0.0) { in KSPSolve_PIPEBCGS()
135 omega = d1 / d2; /* omega <- (y'q) / (y'y) */ in KSPSolve_PIPEBCGS()
146 …PetscCall(VecDotBegin(S, RP, &d1)); … in KSPSolve_PIPEBCGS()
156 PetscCall(VecDotEnd(S, RP, &d1)); in KSPSolve_PIPEBCGS()
160 …PetscCheck(d2 + beta * d1 - beta * omega * d3 != 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR… in KSPSolve_PIPEBCGS()
163 …alpha = rho / (d2 + beta * d1 - beta * omega * d3); /* alpha <- rho / (d2 + beta * d1 - beta * ome… in KSPSolve_PIPEBCGS()
/petsc/src/dm/dt/space/impls/poly/
H A Dspacepoly.c209 for (PetscInt d1 = 0; d1 < dim; d1++) { in PetscSpaceEvaluate_Polynomial() local
212 derivs[d1]++; in PetscSpaceEvaluate_Polynomial()
215 derivs[d1]--; 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.c215 for (PetscInt d1 = 0; d1 < dim; d1++) { in PetscSpaceEvaluate_Ptrimmed() local
218 derivs[d1]++; in PetscSpaceEvaluate_Ptrimmed()
221 derivs[d1]--; in PetscSpaceEvaluate_Ptrimmed()
225 …for (PetscInt p = 0; p < npoints; p++) H[p * p_strl + b * b_strl + v * v_strl + d1 * d1_strl + d2 … in PetscSpaceEvaluate_Ptrimmed()
/petsc/src/tao/linesearch/impls/gpcglinesearch/
H A Dgpcglinesearch.c32 PetscReal d1, finit, actred, prered, rho, gdx; in TaoLineSearchApply_GPCG() local
77 PetscCall(VecStepBoundInfo(x, s, ls->lower, ls->upper, &rho, &actred, &d1)); in TaoLineSearchApply_GPCG()
78 ls->step = PetscMin(ls->step, d1); in TaoLineSearchApply_GPCG()
/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
104 d1 = d[1]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
141 m1 = uik[1] = -(d1 * u0 + d7 * u1 + d13 * u2 + d19 * u3 + d25 * u4 + d31 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
148 m7 = uik[7] = -(d1 * u6 + d7 * u7 + d13 * u8 + d19 * u9 + d25 * u10 + d31 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
155 m13 = uik[13] = -(d1 * u12 + d7 * u13 + d13 * u14 + d19 * u15 + d25 * u16 + d31 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
162 m19 = uik[19] = -(d1 * u18 + d7 * u19 + d13 * u20 + d19 * u21 + d25 * u22 + d31 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
169 m25 = uik[25] = -(d1 * u24 + d7 * u25 + d13 * u26 + d19 * u27 + d25 * u28 + d31 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
176 m31 = uik[31] = -(d1 * u30 + d7 * u31 + d13 * u32 + d19 * u33 + d25 * u34 + d31 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
/petsc/src/tao/matrix/
H A Dadamat.c352 static PetscErrorCode MatCreateADA(Mat mat, Vec d1, Vec d2, Mat *J) in MatCreateADA() argument
361 ctx->D1 = d1; in MatCreateADA()
363 if (d1) { in MatCreateADA()
364 PetscCall(VecDuplicate(d1, &ctx->W)); in MatCreateADA()
365 PetscCall(PetscObjectReference((PetscObject)d1)); in MatCreateADA()
/petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/
H A Ddfovec.m44 c13 = 1.3d1; c14 = 1.4d1; c29 = 2.9d1; c45 = 4.5d1;
/petsc/src/tao/quadratic/impls/bqpip/
H A Dbqpip.c241 PetscReal d1, d2, ksptol, sigmamu; in TaoSolve_BQPIP() local
305 PetscCall(QPIPComputeNormFromCentralPath(qp, &d1)); in TaoSolve_BQPIP()
434 PetscCall(VecDot(tao->solution, qp->Work, &d1)); 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
2459 d1 = PetscLog10Real(PetscAbsReal(sum - osum)); 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
2572 d1 = 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()
2617 PetscInt dim1, Nc1, Np1, order1, qa, d1; in PetscDTTensorQuadratureCreate() local
2773 for (d1 = 0, d = 0; d1 < dim1; ++d1, ++d) x[qc * dim + d] = x1[qa * dim1 + d1]; in PetscDTTensorQuadratureCreate()
/petsc/lib/petsc/bin/maint/
H A Dtoclapack.sh2316 static doublereal b1, b2, c1, c2, d1, d2, one, zero, rbase;
2363 d1 = a;
2368 if (c1 == a && c2 == a && d1 == a && d2 == a) {
2375 d1 = zero;
2378 d1 += b1;
4307 static halfreal b1, b2, c1, c2, d1, d2, one, zero, rbase;
4354 d1 = a;
4359 if (c1 == a && c2 == a && d1 == a && d2 == a) {
4366 d1 = zero;
4369 d1 += b1;
/petsc/src/tao/pde_constrained/tutorials/
H A Delliptic.c185 PetscReal d1 = 0, d2 = 0; in FormFunction() local
192 PetscCall(VecDot(user->dwork, user->dwork, &d1)); in FormFunction()
196 *f = 0.5 * (d1 + user->alpha * d2); in FormFunction()
224 PetscReal d1, d2; in FormFunctionGradient() local
231 PetscCall(VecDot(user->dwork, user->dwork, &d1)); in FormFunctionGradient()
239 *f = 0.5 * (d1 + user->alpha * d2); in FormFunctionGradient()

12