| /petsc/src/snes/tutorials/ |
| H A D | ex4.c | 119 PetscScalar f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormObjectiveLocal() local 151 d4 = (xc - xb); in FormObjectiveLocal() 156 d4 /= hy; in FormObjectiveLocal() 158 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal() 183 d4 = (user->right[j] - user->right[j + 1]) / hy; in FormObjectiveLocal() 184 ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal() 190 d4 = (user->top[i + 1] - user->top[i]) / hx; in FormObjectiveLocal() 191 ft += PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormObjectiveLocal() 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 258 d4 = (xc - xb); in FormFunctionLocal() [all …]
|
| 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 209 d4 = (xc - xb); in FormGradient() 216 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormGradient() 220 df6dxc = d4 * hxdhy; in FormGradient() 225 d4 /= hy; in FormGradient() 232 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); 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 347 d4 = (xc - xb) / hy; in FormJacobian() 354 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormJacobian() [all …]
|
| /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 183 d4 = (xc - xb); in FormFunctionGradient() 190 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient() 194 df6dxc = d4 * hxdhy; in FormFunctionGradient() 199 d4 *= rhy; in FormFunctionGradient() 206 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient() 210 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormFunctionGradient() 239 d4 = (user->right[j] - user->right[j + 1]) * rhy; in FormFunctionGradient() 240 ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient() 245 d4 = (user->top[i + 1] - user->top[i]) * rhx; in FormFunctionGradient() [all …]
|
| H A D | minsurf2.c | 170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local 219 d4 = (xc - xb); in FormFunction() 224 d4 *= rhy; in FormFunction() 226 f2 = PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction() 251 d4 = (user->right[j - ys] - user->right[j - ys + 1]) * rhy; in FormFunction() 252 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction() 258 d4 = (user->top[i - xs + 1] - user->top[i - xs]) * rhx; in FormFunction() 259 ft = ft + PetscSqrtReal(1.0 + d1 * d1 + d4 * d4); in FormFunction() 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 363 d4 = (xc - xb); in FormFunctionGradient() [all …]
|
| /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 197 d4 = (xc - xb); in FormConstraints() 204 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormConstraints() 208 df6dxc = d4 * hxdhy; in FormConstraints() 213 d4 /= hy; in FormConstraints() 220 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); 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 322 d4 = (xc - xb) / hy; in FormJacobian() 329 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormJacobian() [all …]
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 73 PetscReal d4, d5, d6, d7, d8 162 d4 = xc - xb 169 df2dxc = d1*hydhx + d4*hxdhy 173 df6dxc = d4*hxdhy 178 d4 = d4*rhy 185 f2 = sqrt(1.0 + d1*d1 + d4*d4) 189 f6 = sqrt(1.0 + d4*d4 + d6*d6) 224 d4 = (right_v(1 + j - ys) - right_v(1 + j - ys + 1))*rhy 225 ft = ft + sqrt(1.0 + d1*d1 + d4*d4) 232 d4 = (top_v(1 + i - xs + 1) - top_v(1 + i - xs))*rhx [all …]
|
| 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 274 d4 = (xc - xb); in FormFunctionGradient() 281 df2dxc = (d1 * hydhx + d4 * hxdhy); in FormFunctionGradient() 285 df6dxc = d4 * hxdhy; in FormFunctionGradient() 290 d4 *= rhy; in FormFunctionGradient() 297 f2 = PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient() 301 f6 = PetscSqrtScalar(1.0 + d4 * d4 + d6 * d6); in FormFunctionGradient() 335 d4 = (right[j - ys] - right[j - ys + 1]) * rhy; in FormFunctionGradient() 336 ft = ft + PetscSqrtScalar(1.0 + d1 * d1 + d4 * d4); in FormFunctionGradient() 342 d4 = (top[i - xs + 1] - top[i - xs]) * rhx; in FormFunctionGradient() [all …]
|
| /petsc/src/dm/ftn-mod/ |
| H A D | petscdmcomposite.h90 | 46 subroutine DMCompositeGetEntries4(dm, d1, d2, d3, d4, ierr) 50 DM d1, d2, d4, d3
|
| /petsc/src/tao/leastsquares/tutorials/matlab/more_wild_probs/ |
| H A D | dfovec.m | 51 y3 = [3.478d4,2.861d4,2.365d4,1.963d4,1.637d4,1.372d4,1.154d4,9.744d3,...
|
| H A D | g_dfovec_1d.m | 42 …y3= [3.478d4, 2.861d4, 2.365d4, 1.963d4, 1.637d4, 1.372d4, 1.154d4, 9.744d3, 8.261d3, 7.03d3, 6.00…
|
| /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 107 d4 = d[4]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 144 m4 = uik[4] = -(d4 * u0 + d10 * u1 + d16 * u2 + d22 * u3 + d28 * u4 + d34 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 151 m10 = uik[10] = -(d4 * u6 + d10 * u7 + d16 * u8 + d22 * u9 + d28 * u10 + d34 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 158 m16 = uik[16] = -(d4 * u12 + d10 * u13 + d16 * u14 + d22 * u15 + d28 * u16 + d34 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 165 m22 = uik[22] = -(d4 * u18 + d10 * u19 + d16 * u20 + d22 * u21 + d28 * u22 + d34 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 172 m28 = uik[28] = -(d4 * u24 + d10 * u25 + d16 * u26 + d22 * u27 + d28 * u28 + d34 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 179 m34 = uik[34] = -(d4 * u30 + d10 * u31 + d16 * u32 + d22 * u33 + d28 * u34 + d34 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 2431 PetscReal lterm, rterm, maxTerm = 0.0, d1, d2, d3, d4; in PetscDTTanhSinhIntegrate() local 2462 if (PetscMax(PetscAbsReal(lterm), PetscAbsReal(rterm)) == 0.0) d4 = 0.0; in PetscDTTanhSinhIntegrate() 2463 else d4 = PetscLog10Real(PetscMax(PetscAbsReal(lterm), PetscAbsReal(rterm))); 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 2580 d4 = 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()
|