| /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
200 d7 = (xlt - xl); in FormConstraints()
216 d7 /= hy; in FormConstraints()
219 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); 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
325 d7 = (xlt - xl) / hy; in FormJacobian()
328 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormJacobian()
335 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
186 d7 = (xlt - xl); in FormFunctionGradient()
202 d7 *= rhy; in FormFunctionGradient()
205 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); 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
371 d7 = (xlt - xl) * rhy; in QuadraticH()
374 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in QuadraticH()
381 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
366 d7 = (xlt - xl); in FormFunctionGradient()
382 d7 *= rhy; in FormFunctionGradient()
385 f1 = PetscSqrtReal(1.0 + d1 * d1 + d7 * d7); 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
580 d7 = (xlt - xl) / hy; in QuadraticH()
583 f1 = PetscSqrtReal(1.0 + d1 * d1 + d7 * d7); in QuadraticH()
590 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
212 d7 = (xlt - xl); in FormGradient()
228 d7 /= hy; in FormGradient()
231 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); 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
350 d7 = (xlt - xl) / hy; in FormJacobian()
353 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormJacobian()
360 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
261 d7 = (xlt - xl); in FormFunctionLocal()
277 d7 /= hy; in FormFunctionLocal()
280 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); 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
361 d7 = (xlt - xl) / hy; in FormJacobianLocal()
364 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormJacobianLocal()
371 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
165 d7 = xlt - xl
181 d7 = d7*rhy
184 f1 = sqrt(1.0 + d1*d1 + d7*d7)
307 PetscReal d4, d5, d6, d7, d8
404 d7 = (xlt - xl)*rhy
407 f1 = sqrt(1.0 + d1*d1 + d7*d7)
414 hl = (-hydhx*(1.0 + d7*d7) + d1*d7)/(f1*f1*f1) + (-hydhx*(1.0 + d4*d4) + d1*d4)/(f2*f2*f2)
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
277 d7 = (xlt - xl); in FormFunctionGradient()
293 d7 *= rhy; in FormFunctionGradient()
296 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); 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
500 d7 = (xlt - xl) * rhy; in FormHessian()
503 f1 = PetscSqrtScalar(1.0 + d1 * d1 + d7 * d7); in FormHessian()
510 …hl = (-hydhx * (1.0 + d7 * d7) + d1 * d7) / (f1 * f1 * f1) + (-hydhx * (1.0 + d4 * d4) + d1 * d4) … 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
110 d7 = d[7]; 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()
|