| /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 199 d6 = (xrb - xb); in FormConstraints() 215 d6 /= hx; in FormConstraints() 224 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); 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 324 d6 = (xrb - xb) / hx; in FormJacobian() 333 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormJacobian() 338 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobian() 340 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobian() 343 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… 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 185 d6 = (xrb - xb); in FormFunctionGradient() 201 d6 *= rhx; in FormFunctionGradient() 210 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); 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 370 d6 = (xrb - xb) * rhx; in QuadraticH() 379 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in QuadraticH() 384 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in QuadraticH() 386 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in QuadraticH() 389 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… 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 365 d6 = (xrb - xb); in FormFunctionGradient() 381 d6 *= rhx; in FormFunctionGradient() 390 f6 = PetscSqrtReal(1.0 + d4 * d4 + d6 * d6); 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 579 d6 = (xrb - xb) / hx; in QuadraticH() 588 f6 = PetscSqrtReal(1.0 + d4 * d4 + d6 * d6); in QuadraticH() 593 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in QuadraticH() 595 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in QuadraticH() 598 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… 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 211 d6 = (xrb - xb); in FormGradient() 227 d6 /= hx; in FormGradient() 236 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); 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 349 d6 = (xrb - xb) / hx; in FormJacobian() 358 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormJacobian() 363 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobian() 365 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobian() 368 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… 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 260 d6 = (xrb - xb); in FormFunctionLocal() 276 d6 /= hx; in FormFunctionLocal() 285 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); 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 360 d6 = (xrb - xb) / hx; in FormJacobianLocal() 369 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormJacobianLocal() 374 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormJacobianLocal() 376 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormJacobianLocal() 379 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… in FormJacobianLocal()
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 73 PetscReal d4, d5, d6, d7, d8 164 d6 = xrb - xb 180 d6 = d6*rhx 189 f6 = sqrt(1.0 + d4*d4 + d6*d6) 307 PetscReal d4, d5, d6, d7, d8 403 d6 = (xrb - xb)*rhx 412 f6 = sqrt(1.0 + d4*d4 + d6*d6) 420 hb = (-hxdhy*(1.0 + d6*d6) + d4*d6)/(f6*f6*f6) + (-hxdhy*(1.0 + d1*d1) + d1*d4)/(f2*f2*f2) 422 hbr = -d2*d5/(f5*f5*f5) - d4*d6/(f6*f6*f6) 426 …& hxdhy*(1.0 + d6*d6)/(f6*f6*f6) + (hxdhy*(1.0 + d1*d1) + hydhx*(1.0 + d4*d4) - 2*d1*d4)…
|
| 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 276 d6 = (xrb - xb); in FormFunctionGradient() 292 d6 *= rhx; in FormFunctionGradient() 301 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); 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 499 d6 = (xrb - xb) * rhx; in FormHessian() 508 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormHessian() 513 …hb = (-hxdhy * (1.0 + d6 * d6) + d4 * d6) / (f6 * f6 * f6) + (-hxdhy * (1.0 + d1 * d1) + d1 * d4) … in FormHessian() 515 hbr = -d2 * d5 / (f5 * f5 * f5) - d4 * d6 / (f6 * f6 * f6); in FormHessian() 518 …3 * f3 * f3) + hydhx * (1.0 + d5 * d5) / (f5 * f5 * f5) + hxdhy * (1.0 + d6 * d6) / (f6 * f6 * f6)… 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 109 d6 = d[6]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 140 m0 = uik[0] = -(d0 * u0 + d6 * u1 + d12 * u2 + d18 * u3 + d24 * u4 + d30 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 147 m6 = uik[6] = -(d0 * u6 + d6 * u7 + d12 * u8 + d18 * u9 + d24 * u10 + d30 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 154 m12 = uik[12] = -(d0 * u12 + d6 * u13 + d12 * u14 + d18 * u15 + d24 * u16 + d30 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 161 m18 = uik[18] = -(d0 * u18 + d6 * u19 + d12 * u20 + d18 * u21 + d24 * u22 + d30 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 168 m24 = uik[24] = -(d0 * u24 + d6 * u25 + d12 * u26 + d18 * u27 + d24 * u28 + d30 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 175 m30 = uik[30] = -(d0 * u30 + d6 * u31 + d12 * u32 + d18 * u33 + d24 * u34 + d30 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|