| /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 150 d3 = (xc - xt); in FormObjectiveLocal() 155 d3 /= hy; in FormObjectiveLocal() 159 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormObjectiveLocal() 168 d3 = (user->left[j + 1] - user->left[j + 2]) / hy; in FormObjectiveLocal() 170 ft += PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormObjectiveLocal() 176 d3 = (user->bottom[i + 1] - x[0][i]) / hy; in FormObjectiveLocal() 177 ft += PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); 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 257 d3 = (xc - xt); 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 208 d3 = (xc - xt); in FormGradient() 217 df3dxc = d3 * hxdhy; in FormGradient() 218 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormGradient() 224 d3 /= hy; in FormGradient() 233 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormGradient() 234 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); 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 346 d3 = (xc - xt) / hy; in FormJacobian() 355 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); 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 182 d3 = (xc - xt); in FormFunctionGradient() 191 df3dxc = d3 * hxdhy; in FormFunctionGradient() 192 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient() 198 d3 *= rhy; in FormFunctionGradient() 207 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormFunctionGradient() 208 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient() 226 d3 = (user->left[j + 1] - user->left[j + 2]) * rhy; in FormFunctionGradient() 228 ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient() 233 d3 = (user->bottom[i + 1] - x[i]) * rhy; in FormFunctionGradient() [all …]
|
| H A D | minsurf2.c | 170 PetscReal f2, f4, d1, d2, d3, d4, xc, xl, xr, xt, xb; in FormFunction() local 218 d3 = (xc - xt); in FormFunction() 223 d3 *= rhy; in FormFunction() 227 f4 = PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction() 236 d3 = (user->left[j - ys + 1] - user->left[j - ys + 2]) * rhy; in FormFunction() 238 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); in FormFunction() 244 d3 = (user->bottom[i - xs + 1] - x[0][i]) * rhy; in FormFunction() 245 ft = ft + PetscSqrtReal(1.0 + d3 * d3 + d2 * d2); 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 362 d3 = (xc - xt); 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 196 d3 = (xc - xt); in FormConstraints() 205 df3dxc = d3 * hxdhy; in FormConstraints() 206 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormConstraints() 212 d3 /= hy; in FormConstraints() 221 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormConstraints() 222 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); 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 321 d3 = (xc - xt) / hy; in FormJacobian() 330 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormJacobian() [all …]
|
| /petsc/src/tao/bound/tutorials/ |
| H A D | plate2f.F90 | 72 PetscReal f1, f2, f3, f4, f5, f6, d1, d2, d3 161 d3 = xc - xt 170 df3dxc = d3*hxdhy 171 df4dxc = d2*hydhx + d3*hxdhy 177 d3 = d3*rhy 186 f3 = sqrt(1.0 + d3*d3 + d8*d8) 187 f4 = sqrt(1.0 + d3*d3 + d2*d2) 207 d3 = (left_v(1 + j - ys + 1) - left_v(1 + j - ys + 2))*rhy 209 ft = ft + sqrt(1.0 + d3*d3 + d2*d2) 216 d3 = (bottom_v(1 + i - xs + 1) - x_v(1 + i - gxs))*rhy [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 273 d3 = (xc - xt); in FormFunctionGradient() 282 df3dxc = d3 * hxdhy; in FormFunctionGradient() 283 df4dxc = (d2 * hydhx + d3 * hxdhy); in FormFunctionGradient() 289 d3 *= rhy; in FormFunctionGradient() 298 f3 = PetscSqrtScalar(1.0 + d3 * d3 + d8 * d8); in FormFunctionGradient() 299 f4 = PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient() 319 d3 = (left[j - ys + 1] - left[j - ys + 2]) * rhy; in FormFunctionGradient() 321 ft = ft + PetscSqrtScalar(1.0 + d3 * d3 + d2 * d2); in FormFunctionGradient() 327 d3 = (bottom[i - xs + 1] - x[i - gxs]) * rhy; in FormFunctionGradient() [all …]
|
| /petsc/src/dm/ftn-mod/ |
| H A D | petscdmcomposite.h90 | 39 subroutine DMCompositeGetEntries3(dm, d1, d2, d3, ierr) 43 DM d1, d2, d3 46 subroutine DMCompositeGetEntries4(dm, d1, d2, d3, d4, ierr) 50 DM d1, d2, d4, d3
|
| /petsc/src/ts/utils/dmplexlandau/kokkos/ |
| H A D | landau.kokkos.cxx | 614 PetscInt fieldB, d2, d3, f_off_r, grid_r, ipidx_g, nip_loc_r, loc_Nf_r; in LandauKokkosJacobian() local 646 for (d3 = 0; d3 < 2; ++d3) { in LandauKokkosJacobian() 648 ggg.gg2[d2] += Uk[d2][d3] * temp1[d3]; in LandauKokkosJacobian() 650 ggg.gg3[d2][d3] += Ud[d2][d3] * temp2; in LandauKokkosJacobian() 655 for (d3 = 0; d3 < 3; ++d3) { in LandauKokkosJacobian() 657 ggg.gg2[d2] += U[d2][d3]*temp1[d3]; in LandauKokkosJacobian() 659 ggg.gg3[d2][d3] += U[d2][d3]*temp2; in LandauKokkosJacobian() 667 PetscInt d2, d3; in LandauKokkosJacobian() local 670 …for (d3 = 0; d3 < dim; d3++) gg3(d2, d3, fieldA, myQi) = -gg_temp.gg3[d2][d3] * d_alpha[fieldA + f… in LandauKokkosJacobian() 676 int d, d2, d3, dp; in LandauKokkosJacobian() local [all …]
|
| /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,... 52 8.261d3,7.03d3,6.005d3,5.147d3,4.427d3,3.82d3,3.307d3,2.872d3];
|
| H A D | g_dfovec_1d.m | 42 ….963d4, 1.637d4, 1.372d4, 1.154d4, 9.744d3, 8.261d3, 7.03d3, 6.005d3, 5.147d3, 4.427d3, 3.82d3, 3.…
|
| /petsc/share/petsc/saws/js/ |
| H A D | tree.js | 66 var canvas = d3.select("#tree").append("svg") 73 var tree = d3.layout.tree() 113 var diagonal = d3.svg.diagonal()
|
| /petsc/src/ksp/ksp/impls/bcgs/pipebcgs/ |
| H A D | pipebcgs.c | 18 PetscScalar rho, rhoold, alpha, beta, omega = 0.0, d1, d2, d3; in KSPSolve_PIPEBCGS() local 148 …PetscCall(VecDotBegin(Z, RP, &d3)); … in KSPSolve_PIPEBCGS() 158 PetscCall(VecDotEnd(Z, RP, &d3)); in KSPSolve_PIPEBCGS() 160 …PetscCheck(d2 + beta * d1 - beta * omega * d3 != 0.0, PetscObjectComm((PetscObject)ksp), PETSC_ERR… in KSPSolve_PIPEBCGS() 163 …alpha = rho / (d2 + beta * d1 - beta * omega * d3); /* alpha <- rho / (d2 + beta * d1 - beta * ome… in KSPSolve_PIPEBCGS()
|
| /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 106 d3 = d[3]; in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 143 m3 = uik[3] = -(d3 * u0 + d9 * u1 + d15 * u2 + d21 * u3 + d27 * u4 + d33 * u5); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 150 m9 = uik[9] = -(d3 * u6 + d9 * u7 + d15 * u8 + d21 * u9 + d27 * u10 + d33 * u11); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 157 m15 = uik[15] = -(d3 * u12 + d9 * u13 + d15 * u14 + d21 * u15 + d27 * u16 + d33 * u17); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 164 m21 = uik[21] = -(d3 * u18 + d9 * u19 + d15 * u20 + d21 * u21 + d27 * u22 + d33 * u23); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 171 m27 = uik[27] = -(d3 * u24 + d9 * u25 + d15 * u26 + d21 * u27 + d27 * u28 + d33 * u29); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering() 178 m33 = uik[33] = -(d3 * u30 + d9 * u31 + d15 * u32 + d21 * u33 + d27 * u34 + d33 * u35); in MatCholeskyFactorNumeric_SeqSBAIJ_6_NaturalOrdering()
|
| /petsc/src/ts/utils/dmplexlandau/ |
| H A D | plexland.c | 337 PetscInt d, d2, dp, d3, IPf_idx; in LandauFormJacobian_Internal() local 384 for (d3 = 0; d3 < 2; ++d3) { in LandauFormJacobian_Internal() 386 gg2_temp[d2] += Uk[d2][d3] * temp1[d3]; in LandauFormJacobian_Internal() 388 gg3_temp[d2][d3] += Ud[d2][d3] * temp2; in LandauFormJacobian_Internal() 393 for (d3 = 0; d3 < 3; ++d3) { in LandauFormJacobian_Internal() 395 gg2_temp[d2] += U[d2][d3] * temp1[d3]; in LandauFormJacobian_Internal() 397 gg3_temp[d2][d3] += U[d2][d3] * temp2; in LandauFormJacobian_Internal() 410 …for (d3 = 0; d3 < LANDAU_DIM; d3++) gg3[fieldA][d2][d3] = -gg3_temp[d2][d3] * nu_alpha[fieldA + f_… in LandauFormJacobian_Internal() 422 for (d3 = 0; d3 < dim; ++d3) { in LandauFormJacobian_Internal() 423 …for (dp = 0; dp < dim; ++dp) g3[fieldA][d][d2] += invJj[d * dim + d3] * gg3[fieldA][d3][dp] * invJ… in LandauFormJacobian_Internal()
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | wedge_cylinder.cas | 240 (13 (7 d3 de 3 0)
|
| /petsc/src/dm/dt/interface/ |
| H A D | dt.c | 2431 PetscReal lterm, rterm, maxTerm = 0.0, d1, d2, d3, d4; in PetscDTTanhSinhIntegrate() local 2461 d3 = PetscLog10Real(maxTerm) - p; 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 2578 d3 = mpfr_get_d(tmp, MPFR_RNDN) - digits; in PetscDTTanhSinhIntegrateMPFR() 2581 d = PetscAbsInt(PetscMin(0, PetscMax(PetscMax(PetscMax(PetscSqr(d1) / d2, 2 * d1), d3), d4))); in PetscDTTanhSinhIntegrateMPFR()
|