| /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
198 d5 = (xr - xrb); in FormConstraints()
214 d5 /= hy; 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
323 d5 = (xrb - xr) / hy; in FormJacobian()
332 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormJacobian()
336 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobian()
340 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobian()
343 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… 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
184 d5 = (xr - xrb); in FormFunctionGradient()
200 d5 *= rhy; in FormFunctionGradient()
209 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); 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
369 d5 = (xrb - xr) * rhy; in QuadraticH()
378 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in QuadraticH()
382 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in QuadraticH()
386 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in QuadraticH()
389 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… 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
364 d5 = (xr - xrb); in FormFunctionGradient()
380 d5 *= rhy; in FormFunctionGradient()
389 f5 = PetscSqrtReal(1.0 + d2 * d2 + d5 * d5); 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
578 d5 = (xrb - xr) / hy; in QuadraticH()
587 f5 = PetscSqrtReal(1.0 + d2 * d2 + d5 * d5); in QuadraticH()
591 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in QuadraticH()
595 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in QuadraticH()
598 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… 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
210 d5 = (xr - xrb); in FormGradient()
226 d5 /= hy; 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
348 d5 = (xrb - xr) / hy; in FormJacobian()
357 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormJacobian()
361 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobian()
365 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobian()
368 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… 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
259 d5 = (xr - xrb); in FormFunctionLocal()
275 d5 /= hy; in FormFunctionLocal()
284 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); 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
359 d5 = (xrb - xr) / hy; in FormJacobianLocal()
368 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormJacobianLocal()
372 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormJacobianLocal()
376 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobianLocal()
379 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… in FormJacobianLocal()
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 73 PetscReal d4, d5, d6, d7, d8
163 d5 = xr - xrb
179 d5 = d5*rhy
188 f5 = sqrt(1.0 + d2*d2 + d5*d5)
307 PetscReal d4, d5, d6, d7, d8
402 d5 = (xrb - xr)*rhy
411 f5 = sqrt(1.0 + d2*d2 + d5*d5)
416 hr = (-hydhx*(1.0 + d5*d5) + d2*d5)/(f5*f5*f5) + (-hydhx*(1.0 + d3*d3) + d2*d3)/(f4*f4*f4)
422 hbr = -d2*d5/(f5*f5*f5) - d4*d6/(f6*f6*f6)
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
275 d5 = (xr - xrb); in FormFunctionGradient()
291 d5 *= rhy; in FormFunctionGradient()
300 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); 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
498 d5 = (xrb - xr) * rhy; in FormHessian()
507 f5 = PetscSqrtScalar(1.0 + d2 * d2 + d5 * d5); in FormHessian()
511 …hr = (-hydhx * (1.0 + d5 * d5) + d2 * d5) / (f5 * f5 * f5) + (-hydhx * (1.0 + d3 * d3) + d2 * d3) … in FormHessian()
515 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormHessian()
518 …1 * f1 * f1) + hxdhy * (1.0 + d8 * d8) / (f3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5)… 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
108 d5 = d[5]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
145 m5 = uik[5] = -(d5 * u0 + d11 * u1 + d17 * u2 + d23 * u3 + d29 * u4 + d35 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
152 m11 = uik[11] = -(d5 * u6 + d11 * u7 + d17 * u8 + d23 * u9 + d29 * u10 + d35 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
159 m17 = uik[17] = -(d5 * u12 + d11 * u13 + d17 * u14 + d23 * u15 + d29 * u16 + d35 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
166 m23 = uik[23] = -(d5 * u18 + d11 * u19 + d17 * u20 + d23 * u21 + d29 * u22 + d35 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
173 m29 = uik[29] = -(d5 * u24 + d11 * u25 + d17 * u26 + d23 * u27 + d29 * u28 + d35 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
180 m35 = uik[35] = -(d5 * u30 + d11 * u31 + d17 * u32 + d23 * u33 + d29 * u34 + d35 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|