| /petsc/src/tao/complementarity/tutorials/ |
| H A D | minsurf1.c | 142 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
201 d8 = (xt - xlt); in FormConstraints()
217 d8 /= hx; in FormConstraints()
221 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); 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
326 d8 = (xlt - xt) / hx; in FormJacobian()
330 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormJacobian()
337 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobian()
341 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in FormJacobian()
343 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in FormJacobian()
|
| /petsc/src/tao/unconstrained/tutorials/ |
| H A D | minsurf1.c | 131 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in FormFunctionGradient() local
187 d8 = (xt - xlt); in FormFunctionGradient()
203 d8 *= rhx; in FormFunctionGradient()
207 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormFunctionGradient()
310 …PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt… in QuadraticH() local
372 d8 = (xlt - xt) * rhx; in QuadraticH()
376 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in QuadraticH()
383 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in QuadraticH()
387 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in QuadraticH()
389 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in QuadraticH()
|
| H A D | minsurf2.c | 304 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
367 d8 = (xt - xlt); in FormFunctionGradient()
383 d8 *= rhx; in FormFunctionGradient()
387 f3 = PetscSqrtReal(1.0 + d3 * d3 + d8 * d8); in FormFunctionGradient()
508 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in QuadraticH() local
581 d8 = (xlt - xt) / hx; in QuadraticH()
585 f3 = PetscSqrtReal(1.0 + d3 * d3 + d8 * d8); in QuadraticH()
592 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in QuadraticH()
596 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in QuadraticH()
598 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in QuadraticH()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex58.c | 149 …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
213 d8 = (xt - xlt); in FormGradient()
229 d8 /= hx; in FormGradient()
233 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); 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
351 d8 = (xlt - xt) / hx; in FormJacobian()
355 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormJacobian()
362 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobian()
366 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in FormJacobian()
368 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in FormJacobian()
|
| H A D | ex4.c | 217 PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormFunctionLocal() local
262 d8 = (xt - xlt); in FormFunctionLocal()
278 d8 /= hx; in FormFunctionLocal()
282 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormFunctionLocal()
309 PetscScalar f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormJacobianLocal() local
362 d8 = (xlt - xt) / hx; in FormJacobianLocal()
366 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormJacobianLocal()
373 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormJacobianLocal()
377 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in FormJacobianLocal()
379 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in FormJacobianLocal()
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 73 PetscReal d4, d5, d6, d7, d8
166 d8 = xt - xlt
182 d8 = d8*rhx
186 f3 = sqrt(1.0 + d3*d3 + d8*d8)
307 PetscReal d4, d5, d6, d7, d8
405 d8 = (xlt - xt)*rhx
409 f3 = sqrt(1.0 + d3*d3 + d8*d8)
418 ht = (-hxdhy*(1.0 + d8*d8) + d3*d8)/(f3*f3*f3) + (-hxdhy*(1.0 + d2*d2) + d2*d3)/(f4*f4*f4)
423 htl = -d1*d7/(f1*f1*f1) - d3*d8/(f3*f3*f3)
425 …hc = hydhx*(1.0 + d7*d7)/(f1*f1*f1) + hxdhy*(1.0 + d8*d8)/(f3*f3*f3) + hydhx*(1.0 + d5*d5)/(f5*f5*…
|
| H A D | plate2.c | 207 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
278 d8 = (xt - xlt); in FormFunctionGradient()
294 d8 *= rhx; in FormFunctionGradient()
298 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormFunctionGradient()
423 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3, d4, d5, d6, d7, d8, xc, xl, xr, xt, xb, xlt, xrb; in FormHessian() local
501 d8 = (xlt - xt) * rhx; in FormHessian()
505 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormHessian()
512 …ht = (-hxdhy * (1.0 + d8 * d8) + d3 * d8) / (f3 * f3 * f3) + (-hxdhy * (1.0 + d2 * d2) + d2 * d3) … in FormHessian()
516 htl = -d1 * d7 / (f1 * f1 * f1) - d3 * d8 / (f3 * f3 * f3); in FormHessian()
518 …hc = hydhx * (1.0 + d7 * d7) / (f1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx *… in FormHessian()
|
| /petsc/src/mat/impls/sbaij/seq/ |
| H A D | sbaijfact10.c | 16 MatScalar d0, d1, d2, d3, d4, d5, d6, d7, d8, d9, d10, d11, d12; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() local
111 d8 = d[8]; 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()
|