| /petsc/src/ts/tutorials/ |
| H A D | ex53.c | 452 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_u() local
465 …* (1.0 - nu_u))) * (1.0 - zstar) + ((P_0 * L * (nu_u - nu)) / (2.0 * G * (1.0 - nu_u) * (1.0 - nu)… in terzaghi_2d_u()
489 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_eps() local
501 … / (2.0 * G * (1.0 - nu_u) * L)) + ((P_0 * L * (nu_u - nu)) / (2.0 * G * (1.0 - nu_u) * (1.0 - nu)… in terzaghi_2d_eps()
564 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_u_t() local
577 u[1] = ((P_0 * L * (nu_u - nu)) / (2.0 * G * (1.0 - nu_u) * (1.0 - nu))) * F2_t; /* m / s */ in terzaghi_2d_u_t()
601 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_eps_t() local
613 u[0] = ((P_0 * L * (nu_u - nu)) / (2.0 * G * (1.0 - nu_u) * (1.0 - nu))) * F2_zt; /* 1 / s */ in terzaghi_2d_eps_t()
724 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in mandel_initial_u() local
736 u[0] = ((P_0 * nu) / (2.0 * G * a) - (P_0 * nu_u) / (G * a) * A_s) * x[0] + P_0 / G * B_s; in mandel_initial_u()
[all …]
|
| H A D | ex76.c | 87 PetscReal nu; /* Kinematic viscosity */ member
177 const PetscReal nu = PetscRealPart(constants[NU]); in f0_quadratic_v() local
179 … * X[1]) + 2 * X[0] * X[0] * X[0] + 4 * X[0] * X[0] * X[1] - 2 * X[0] * X[1] * X[1] - 4.0 * nu + 2; in f0_quadratic_v()
180 … * X[1]) + 4 * X[0] * X[1] * X[1] + 2 * X[0] * X[0] * X[1] - 2 * X[1] * X[1] * X[1] - 4.0 * nu + 2; in f0_quadratic_v()
304 const PetscReal nu = PetscRealPart(constants[NU]); in f0_cubic_v() local
306 … * X[1] * X[1] - 6 * X[0] * X[0] * X[1] * X[1] * X[1] - (6 * X[0] + 6 * X[1]) * nu + 3 * X[0] + 1); in f0_cubic_v()
307 …* X[1] * X[1] - 6 * X[0] * X[1] * X[1] * X[1] * X[1] - (12 * X[0] - 6 * X[1]) * nu + 3 * X[1] + 1); in f0_cubic_v()
376 const PetscReal nu = PetscRealPart(constants[NU]); in f0_cubic_trig_v() local
378 … X[0] * X[1] * X[1] - 6 * X[0] * X[0] * X[1] * X[1] * X[1] - (6 * X[0] + 6 * X[1]) * nu + 3 * X[0]; in f0_cubic_trig_v()
379 …X[1] * X[1] * X[1] - 6 * X[0] * X[1] * X[1] * X[1] * X[1] - (12 * X[0] - 6 * X[1]) * nu + 3 * X[1]; in f0_cubic_trig_v()
[all …]
|
| H A D | ex77.c | 48 PetscReal nu; /* Kinematic viscosity */ member
193 const PetscReal nu = PetscRealPart(constants[0]); in f1_v() local
198 for (d = 0; d < dim; ++d) f1[c * dim + d] = nu * (u_x[c * dim + d] + u_x[d * dim + c]); in f1_v()
252 const PetscReal nu = PetscRealPart(constants[0]); in g3_vu() local
258 g3[((c * Nc + c) * dim + d) * dim + d] += nu; // gradU in g3_vu()
259 g3[((c * Nc + d) * dim + d) * dim + c] += nu; // gradU transpose in g3_vu()
328 PetscCall(PetscBagRegisterReal(bag, &p->nu, 1.0, "nu", "Kinematic viscosity")); in SetupParameters()
387 constants[0] = param->nu; in SetupProblem()
|
| H A D | ex11.h | 364 PetscReal nu, p; in EulerFlux() local
369 nu = DotDIMReal(x->ru, n); in EulerFlux()
370 f->r = nu; /* A rho u */ in EulerFlux()
371 nu /= x->r; /* A u */ in EulerFlux()
372 for (i = 0; i < DIM; i++) f->ru[i] = nu * x->ru[i] + n[i] * p; /* r u^2 + p */ in EulerFlux()
373 f->E = nu * (x->E + p); /* u(e+p) */ in EulerFlux()
|
| /petsc/src/ksp/ksp/impls/cg/pipeprcg/ |
| H A D | pipeprcg.c | 42 …PetscScalar alpha = 0.0, beta = 0.0, nu = 0.0, nu_old = 0.0, mudelgam[3], *mu_p, *delta_p, *ga… in KSPSolve_PIPEPRCG() local
91 PetscCall(VecDotBegin(RT, R, &nu)); in KSPSolve_PIPEPRCG()
95 PetscCall(VecDotEnd(RT, R, &nu)); /* nu <- (rt,r) */ in KSPSolve_PIPEPRCG()
115 dp = PetscSqrtReal(PetscAbsScalar(nu)); in KSPSolve_PIPEPRCG()
131 alpha = nu / *mu_p; in KSPSolve_PIPEPRCG()
132 nu_old = nu; in KSPSolve_PIPEPRCG()
133 nu = nu_old - 2. * alpha * (*delta_p) + (alpha * alpha) * (*gamma_p); in KSPSolve_PIPEPRCG()
134 beta = nu / nu_old; in KSPSolve_PIPEPRCG()
146 PetscCall(VecDotBegin(RT, R, &nu)); in KSPSolve_PIPEPRCG()
166 PetscCall(VecDotEnd(RT, R, &nu)); in KSPSolve_PIPEPRCG()
|
| /petsc/src/ksp/ksp/tutorials/ |
| H A D | ex29.c | 37 PetscReal nu; member
63 user.nu = 0.1; in main()
64 …PetscCall(PetscOptionsReal("-nu", "The width of the Gaussian source", "ex29.c", user.nu, &user.nu,… in main()
120 …)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScalar(-((PetscReal)j * Hy) * ((PetscReal)j *… in ComputeRHS()
122 …)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScalar(-((PetscReal)j * Hy) * ((PetscReal)j *… in ComputeRHS()
|
| H A D | ex32.c | 42 PetscScalar nu; member
65 user.nu = 0.1; in main()
66 …PetscCall(PetscOptionsScalar("-nu", "The width of the Gaussian source", "ex29.c", 0.1, &user.nu, N… in main()
99 …tscReal)i + 0.5) * Hx) / user->nu) * PetscExpScalar(-(((PetscReal)j + 0.5) * Hy) * (((PetscReal)j … in ComputeRHS()
|
| H A D | ex36.cxx | 81 PetscReal nu; member
231 return PetscExpScalar(-(xx + yy + zz) / user->nu); in ExactSolution()
246 return (4.0 / user->nu * duxyz * duxyz - 6.0) * exact / user->nu; in ForcingFunction()
528 user->nu = 0.05; in InitializeOptions()
546 …-nu", "The width of the Gaussian source (for -problem 1)", "ex36.cxx", user->nu, &user->nu, NULL)); in InitializeOptions()
|
| H A D | ex35.cxx | 98 PetscReal nu; member
254 …return user->nu * sin(PETSC_PI * coords[0] / user->bounds[1]) * sin(PETSC_PI * coords[1] / user->b… in ComputeForcingFunction()
256 …ref) * (coords[0] - user->xref) + (coords[1] - user->yref) * (coords[1] - user->yref)) / user->nu); in ComputeForcingFunction()
572 user->nu = 0.05; in InitializeOptions()
588 …-nu", "The width of the Gaussian source (for -problem 1)", "ex35.cxx", user->nu, &user->nu, NULL)); in InitializeOptions()
|
| H A D | ex49.c | 84 PetscScalar E, nu, fx, fy; member
91 PetscScalar nu[GAUSS_POINTS]; member
405 …lPart(_coefficients[j][i].E[p]), (double)PetscRealPart(_coefficients[j][i].nu[p]), (double)PetscRe… in DMDAViewCoefficientsGnuplot2d()
417 …oid FormStressOperatorQ1(PetscScalar Ke[], PetscScalar coords[], PetscScalar E[], PetscScalar nu[]) in FormStressOperatorQ1() argument
450 prop_nu = nu[p]; in FormStressOperatorQ1()
601 prop_nu = props[ej][ei].nu; in AssembleA_Elasticity()
843 element_props[j][i].nu[p] = opts_nu; in solve_elasticity_2d()
870 element_props[j][i].nu[p] = opts_nu0; in solve_elasticity_2d()
873 element_props[j][i].nu[p] = opts_nu1; in solve_elasticity_2d()
908 element_props[j][i].nu[p] = values_nu[index]; in solve_elasticity_2d()
[all …]
|
| H A D | ex73.c | 40 PetscReal nu; member
57 user->nu = 0.1; in UserContextCreate()
58 …(PetscOptionsReal("-nu", "The width of the Gaussian source", "ex29.c", user->nu, &user->nu, NULL)); in UserContextCreate()
956 …)i * Hx) * ((PetscReal)i * Hx) / user->nu) * PetscExpScalar(-((PetscReal)j * Hy) * ((PetscReal)j *… in ComputeRHS_DMDA()
|
| /petsc/src/snes/tutorials/ |
| H A D | ex71.c | 31 PetscReal nu; /* Kinematic viscosity */ member
72 PetscReal nu = param->nu; in quadratic_u() local
78 for (d = 1; d < dim; ++d) u[0] += Delta / (fac * 2.0 * nu) * X[d] * (1.0 - X[d]); in quadratic_u()
107 const PetscReal nu = PetscRealPart(constants[1]); in f1_u() local
114 f1[c * dim + d] = nu * u_x[c * dim + d]; in f1_u()
158 const PetscReal nu = PetscRealPart(constants[1]); in g3_uu() local
163 for (d = 0; d < dim; ++d) g3[((c * Nc + c) * dim + d) * dim + d] = nu; in g3_uu()
178 PetscCall(PetscBagRegisterReal(bag, &p->nu, 1.0, "nu", "Kinematic viscosity")); in SetupParameters()
249 constants[1] = param->nu; in SetupProblem()
|
| H A D | ex64.c | 135 PetscScalar nu; member
144 user->nu = 0.3; in ProcessOptions()
149 PetscCall(PetscOptionsScalar("-nu", "Poisson ratio", NULL, user->nu, &user->nu, NULL)); in ProcessOptions()
195 constants[0] = user->E / (2.0 * (1.0 + user->nu)); in SetupEqn()
196 constants[1] = user->nu * user->E / ((1.0 + user->nu) * (1.0 - 2.0 * user->nu)); in SetupEqn()
201 …Object)dm), "E = %g, nu = %g\n", (double)PetscRealPart(user->E), (double)PetscRealPart(user->nu))); in SetupEqn()
|
| H A D | ex76.c | 35 PetscReal nu; /* Kinematic viscosity */ member
105 const PetscReal nu = PetscRealPart(constants[0]); in f0_quadratic_v() local
110 …f0[0] -= (2 * X[0] * X[0] * X[0] + 4 * X[0] * X[0] * X[1] - 2 * X[0] * X[1] * X[1] - 4.0 * nu + 1); in f0_quadratic_v()
111 …f0[1] -= (4 * X[0] * X[1] * X[1] + 2 * X[0] * X[0] * X[1] - 2 * X[1] * X[1] * X[1] - 4.0 * nu + 1); in f0_quadratic_v()
165 const PetscReal nu = PetscRealPart(constants[0]); in f0_cubic_v() local
170 …X[0] * X[1] * X[1] - 6 * X[0] * X[0] * X[1] * X[1] * X[1] - (6 * X[0] + 6 * X[1]) * nu + 3 * X[0]); in f0_cubic_v()
171 …[0] * X[0] * X[1] - 6 * X[0] * X[1] * X[1] * X[1] * X[1] - (12 * X[0] - 6 * X[1]) * nu + 3 * X[1]); in f0_cubic_v()
191 const PetscReal nu = PetscRealPart(constants[0]); in f1_v() local
197 f1[c * dim + d] = nu * (u_x[c * dim + d] + u_x[d * dim + c]); in f1_v()
251 const PetscReal nu = PetscRealPart(constants[0]); in g3_vu() local
[all …]
|
| H A D | ex12.c | 314 PetscScalar nu = 0.0; in f1_analytic_nonlinear_u() local
316 for (d = 0; d < dim; ++d) nu += u_x[d] * u_x[d]; in f1_analytic_nonlinear_u()
317 for (d = 0; d < dim; ++d) f1[d] = 0.5 * nu * u_x[d]; in f1_analytic_nonlinear_u()
330 PetscScalar nu = 0.0; in g3_analytic_nonlinear_uu() local
332 for (d = 0; d < dim; ++d) nu += u_x[d] * u_x[d]; in g3_analytic_nonlinear_uu()
334 g3[d * dim + d] = 0.5 * nu; in g3_analytic_nonlinear_uu()
667 Vec nu; in SetupMaterial() local
672 PetscCall(DMCreateLocalVector(dmAux, &nu)); in SetupMaterial()
673 PetscCall(PetscObjectSetName((PetscObject)nu, "Coefficient")); in SetupMaterial()
674 PetscCall(DMProjectFunctionLocal(dmAux, 0.0, matFuncs, ctx, INSERT_ALL_VALUES, nu)); in SetupMaterial()
[all …]
|
| /petsc/src/snes/tutorials/output/ |
| H A D | ex64_2d_p2_p1_p2_fetidp.out | 1 E = 1., nu = 0.3
|
| H A D | ex64_2d_p2_p1_p2.out | 1 E = 1., nu = 0.3
|
| /petsc/src/ksp/ksp/tests/ |
| H A D | ex10.c | 198 PetscReal E, nu; /* Physical constants. */ variable
317 nu = 0.3; in paulsetup20()
460 temp = (1.0 + nu) * (1.0 - 2.0 * nu); in paulintegrate20()
462 C[0][0] = temp * (1.0 - nu); in paulintegrate20()
465 C[3][3] = temp * (0.5 - nu); in paulintegrate20()
468 C[0][1] = temp * nu; in paulintegrate20()
|
| H A D | ex11.c | 79 PetscInt lnu, lnp, nu, np, i, start_u, end_u, start_p, end_p; in LoadTestMatrices() local
88 PetscCall(VecGetSize(f, &nu)); in LoadTestMatrices()
|
| /petsc/src/snes/impls/ngmres/ |
| H A D | ngmresfunc.c | 47 PetscReal nu; in SNESNGMRESFormCombinedSolution_Private() local
52 nu = fMnorm * fMnorm; in SNESNGMRESFormCombinedSolution_Private()
68 for (i = 0; i < l; i++) beta[i] = nu - xi[i]; in SNESNGMRESFormCombinedSolution_Private()
72 for (i = 0; i < l; i++) H(i, j) = Q(i, j) - xi[i] - xi[j] + nu; in SNESNGMRESFormCombinedSolution_Private()
|
| /petsc/include/petsc/private/ |
| H A D | petschpddm.h | 25 PetscInt nu; /* number of local deflation vectors */ member
|
| /petsc/src/ksp/pc/impls/fieldsplit/ |
| H A D | fieldsplit.c | 1866 PetscReal lowbnd, nu, beta; in PCApply_FieldSplit_GKB() local
1886 nu = jac->gkbnu; in PCApply_FieldSplit_GKB()
1892 nu = 1; in PCApply_FieldSplit_GKB()
1901 PetscCall(VecAXPBY(work1, 1.0 / nu, -1.0, ilinkD->x)); /* c = b - B'*x */ in PCApply_FieldSplit_GKB()
1906 beta = PetscSqrtReal(nu) * beta; in PCApply_FieldSplit_GKB()
1907 PetscCall(VecAXPBY(v, nu / beta, 0.0, work1)); /* v = nu/beta *c */ in PCApply_FieldSplit_GKB()
1936 PetscCall(VecAXPBY(v, nu, -alpha, work1)); in PCApply_FieldSplit_GKB()
1938 beta = beta / PetscSqrtReal(nu); in PCApply_FieldSplit_GKB()
3285 PetscErrorCode PCFieldSplitSetGKBNu(PC pc, PetscReal nu) in PCFieldSplitSetGKBNu() argument
3289 PetscValidLogicalCollectiveReal(pc, nu, 2); in PCFieldSplitSetGKBNu()
[all …]
|
| /petsc/src/ksp/pc/impls/hpddm/ |
| H A D | pchpddm.cxx | 354 data->levels[i - 1]->nu = 0; in PCSetFromOptions_HPDDM()
357 … computed by SLEPc", "EPSSetDimensions", data->levels[i - 1]->nu, &data->levels[i - 1]->nu, nullpt… in PCSetFromOptions_HPDDM()
370 …PetscCheck(data->levels[0]->nu == 0 || nsv == 0, PETSC_COMM_SELF, PETSC_ERR_ARG_WRONG, "Cannot sup… in PCSetFromOptions_HPDDM()
371 if (data->levels[0]->nu == 0) { /* -eps_nev has not been used, so nu is 0 */ in PCSetFromOptions_HPDDM()
372 data->levels[0]->nu = nsv; /* nu may still be 0 if -svd_nsv has not been used */ in PCSetFromOptions_HPDDM()
378 … if (data->levels[0]->nu == 0 || nsv == 0) { /* if neither -eps_nev nor -svd_nsv has been used */ in PCSetFromOptions_HPDDM()
385 …PetscCheck(data->levels[0]->nu || PetscAbsReal(data->levels[i - 1]->threshold + static_cast<PetscR… in PCSetFromOptions_HPDDM()
403 …if (data->levels[i - 1]->threshold <= PetscReal() && data->levels[i - 1]->nu <= 0 && !(data->defla… in PCSetFromOptions_HPDDM()
584 PetscCall(PetscViewerASCIIPrintf(viewer, " %" PetscInt_FMT, data->levels[i - 1]->nu)); in PCView_HPDDM()
|
| /petsc/doc/tutorials/physics/ |
| H A D | guide_to_stokes.md | 42 …we divide by the constant density, we would have the kinematic viscosity $\nu$ and a force per uni…
|
| /petsc/share/petsc/datafiles/meshes/ |
| H A D | testcase3D.cas | 2313 (dpm/spray-suite/ssd-nu-t 1.)
4679 (reacting-channel/ud-nu-num 3.66)
4680 (reacting-channel/use-ud-nu? #f)
4681 (reacting-channel/nu-entry-len-correction? #f)
8839 …nu): ~a")) (k/residual-tolerance 0.001 real #f 0.001 0.0001 ("Residual-tolerance for turb…
|