Lines Matching refs:G
149 const PetscReal G = PetscRealPart(constants[0]); in f0_quadratic_linear_u() local
154 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_quadratic_linear_u()
157 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * t; in f0_quadratic_linear_u()
158 f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * t; in f0_quadratic_linear_u()
231 const PetscReal G = PetscRealPart(constants[0]); in f0_quadratic_trig_u() local
236 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_quadratic_trig_u()
239 for (d = 0; d < dim - 1; ++d) f0[d] -= 2.0 * G - alpha * PetscCosReal(t); in f0_quadratic_trig_u()
240 f0[dim - 1] -= 2.0 * lambda + 4.0 * G - alpha * PetscCosReal(t); in f0_quadratic_trig_u()
344 const PetscReal G = PetscRealPart(constants[0]); in f0_trig_linear_u() local
349 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f0_trig_linear_u()
352 …cSqr(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[d]) * (2. * G + lambda) + 2.0 * (G + lambda) … in f0_trig_linear_u()
353 …f0[dim - 1] += PetscSqr(2. * PETSC_PI) * PetscSinReal(2. * PETSC_PI * x[dim - 1]) * (2. * G + lamb… in f0_trig_linear_u()
383 PetscScalar G = param->mu; /* Pa */ in terzaghi_drainage_pressure() local
386 …PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) … in terzaghi_drainage_pressure()
387 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) … in terzaghi_drainage_pressure()
389 u[0] = ((P_0 * eta) / (G * S)); in terzaghi_drainage_pressure()
404 PetscScalar G = param->mu; /* Pa */ in terzaghi_initial_u() local
407 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_initial_u()
411 u[1] = ((P_0 * L * (1.0 - 2.0 * nu_u)) / (2.0 * G * (1.0 - nu_u))) * (1.0 - zstar); in terzaghi_initial_u()
424 PetscScalar G = param->mu; /* Pa */ in terzaghi_initial_eps() local
426 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_initial_eps()
428 u[0] = -(P_0 * (1.0 - 2.0 * nu_u)) / (2.0 * G * (1.0 - nu_u)); in terzaghi_initial_eps()
445 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_u() local
452 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_u()
453 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_2d_u()
454 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_u()
465 …(P_0 * L * (1.0 - 2.0 * nu_u)) / (2.0 * G * (1.0 - nu_u))) * (1.0 - zstar) + ((P_0 * L * (nu_u - n… in terzaghi_2d_u()
482 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_eps() local
489 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in terzaghi_2d_eps()
490 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_2d_eps()
491 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_eps()
501 …] = -((P_0 * L * (1.0 - 2.0 * nu_u)) / (2.0 * G * (1.0 - nu_u) * L)) + ((P_0 * L * (nu_u - nu)) / … in terzaghi_2d_eps()
519 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p() local
526 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p()
527 … PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_p()
534 …/ G) - S) <= 1.0e-10, PETSC_COMM_SELF, PETSC_ERR_PLIB, "S %g != check %g", (double)PetscAbsScalar(… in terzaghi_2d_p()
539 u[0] = ((P_0 * eta) / (G * S)) * F1; /* Pa */ in terzaghi_2d_p()
557 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_u_t() local
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()
565 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_2d_u_t()
566 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_u_t()
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()
594 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_eps_t() local
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()
602 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in terzaghi_2d_eps_t()
603 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_eps_t()
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()
628 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p_t() local
634 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p_t()
635 … PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_p_t()
638 u[0] = -((P_0 * eta) / (G * S)) * PetscSqr(0 * PETSC_PI) * c / PetscSqr(2.0 * L); /* Pa / s */ in terzaghi_2d_p_t()
643 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p_t() local
650 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p_t()
651 … PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in terzaghi_2d_p_t()
658 …/ G) - S) <= 1.0e-10, PETSC_COMM_SELF, PETSC_ERR_PLIB, "S %g != check %g", (double)PetscAbsScalar(… in terzaghi_2d_p_t()
663 u[0] = ((P_0 * eta) / (G * S)) * F1_t; /* Pa / s */ in terzaghi_2d_p_t()
679 PetscScalar G = param->mu; /* Pa */ in mandel_drainage_pressure() local
686 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in mandel_drainage_pressure()
688 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in mandel_drainage_pressure()
717 PetscScalar G = param->mu; /* Pa */ in mandel_initial_u() local
724 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in mandel_initial_u()
725 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in mandel_initial_u()
726 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in mandel_initial_u()
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()
737 u[1] = (-1 * (P_0 * (1.0 - nu)) / (2 * G * a) + (P_0 * (1 - nu_u)) / (G * a) * A_s) * x[1]; in mandel_initial_u()
752 PetscScalar G = param->mu; /* Pa */ in mandel_initial_eps() local
759 … PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in mandel_initial_eps()
760 … PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in mandel_initial_eps()
775 …u[0] = (P_0 / G) * eps_A + ((P_0 * nu) / (2.0 * G * a)) - eps_B / (G * a) - (P_0 * (1 - nu)) / (2 … in mandel_initial_eps()
795 PetscScalar G = param->mu; in mandel_2d_u() local
801 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_u()
802 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_u()
805 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha … in mandel_2d_u()
817 u[0] = ((F * nu) / (2.0 * G * a) - (F * nu_u) / (G * a) * A_x) * x[0] + F / G * B_x; in mandel_2d_u()
818 u[1] = (-1 * (F * (1.0 - nu)) / (2 * G * a) + (F * (1 - nu_u)) / (G * a) * A_x) * x[1]; in mandel_2d_u()
838 PetscScalar G = param->mu; in mandel_2d_eps() local
844 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_eps()
845 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_eps()
851 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha * (1.0… in mandel_2d_eps()
868 …u[0] = (F / G) * eps_A + ((F * nu) / (2.0 * G * a)) - eps_B / (G * a) - (F * (1 - nu)) / (2 * G * … in mandel_2d_eps()
889 PetscScalar G = param->mu; in mandel_2d_p() local
895 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_p()
896 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_p()
901 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha * (… in mandel_2d_p()
930 PetscScalar G = param->mu; in mandel_2d_u_t() local
934 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_u_t()
935 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_u_t()
938 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha … in mandel_2d_u_t()
951 u[0] = (F / G) * A_s_t - ((F * nu_u * x[0]) / (G * a)) * B_s_t; in mandel_2d_u_t()
952 u[1] = ((F * x[1] * (1 - nu_u)) / (G * a)) * B_s_t; in mandel_2d_u_t()
970 PetscScalar G = param->mu; in mandel_2d_eps_t() local
976 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_eps_t()
977 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_eps_t()
983 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha * (1.0… in mandel_2d_eps_t()
998 u[0] = (F / G) * eps_As - ((F * nu_u) / (G * a)) * eps_Bs + ((F * (1 - nu_u)) / (G * a)) * eps_Cs; in mandel_2d_eps_t()
1014 PetscScalar G = param->mu; in mandel_2d_p_t() local
1018 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandel_2d_p_t()
1019 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandel_2d_p_t()
1059 PetscScalar G = param->mu; /* Pa */ in cryer_initial_u() local
1062 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in cryer_initial_u()
1064 PetscScalar u_0 = -P_0 * R_0 * (1. - 2. * nu_u) / (2. * G * (1. + nu_u)); /* Cheng (7.407) */ in cryer_initial_u()
1082 PetscScalar G = param->mu; /* Pa */ in cryer_initial_eps() local
1085 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in cryer_initial_eps()
1087 PetscScalar u_0 = -P_0 * R_0 * (1. - 2. * nu_u) / (2. * G * (1. + nu_u)); /* Cheng (7.407) */ in cryer_initial_eps()
1111 PetscScalar G = param->mu; /* Pa */ in cryer_3d_u() local
1118 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) … in cryer_3d_u()
1119 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) … in cryer_3d_u()
1120 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) … in cryer_3d_u()
1122 …PetscScalar u_inf = -P_0 * R_0 * (1. - 2. * nu) / (2. * G * (1. + nu)); /* m, Cheng (7.388)… in cryer_3d_u()
1161 PetscScalar G = param->mu; /* Pa */ in cryer_3d_eps() local
1168 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) … in cryer_3d_eps()
1169 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) … in cryer_3d_eps()
1170 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) … in cryer_3d_eps()
1172 …PetscScalar u_inf = -P_0 * R_0 * (1. - 2. * nu) / (2. * G * (1. + nu)); /* m, Cheng (7.388)… in cryer_3d_eps()
1213 PetscScalar G = param->mu; /* Pa */ in cryer_3d_p() local
1220 …PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in cryer_3d_p()
1221 …PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in cryer_3d_p()
1222 …PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in cryer_3d_p()
1223 …PetscScalar S = (3.0 * K_u + 4.0 * G) / (M * (3.0 * K_d + 4.0 * G)); /* Pa^{-1}, Cheng (B.14) */ in cryer_3d_p()
1289 PetscScalar G = constants[0];
1297 PetscScalar nu = (3.0*K_d - 2.0*G) / (2.0*(3.0*K_d + G));
1298 PetscScalar nu_u = (3.0*K_u - 2.0*G) / (2.0*(3.0*K_u + G));
1300 …PetscReal c = PetscRealPart(((2.0*kappa*G) * (1.0 - nu) * (nu_u - nu)) / (alpha*alpha * (1.0 - 2…
1385 const PetscReal G = PetscRealPart(constants[0]); in f1_u() local
1390 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in f1_u()
1394 for (d = 0; d < dim; ++d) f1[c * dim + d] -= G * (u_x[c * dim + d] + u_x[d * dim + c]); in f1_u()
1449 const PetscReal G = PetscRealPart(constants[0]); in g2_ue() local
1454 const PetscReal lambda = K_d - (2.0 * G) / 3.0; in g2_ue()
1472 const PetscReal G = PetscRealPart(constants[0]); in g3_uu() local
1477 g3[((c * Nc + c) * dim + d) * dim + d] -= G; in g3_uu()
1478 g3[((c * Nc + d) * dim + d) * dim + c] -= G; in g3_uu()
1527 PetscScalar G = param->mu; in mandelZeros() local
1533 PetscScalar nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); in mandelZeros()
1534 PetscScalar nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); in mandelZeros()
1569 PetscReal G = PetscRealPart(param->mu); /* Pa */ in cryerZeros() local
1573 PetscReal nu = (3.0 * K_d - 2.0 * G) / (2.0 * (3.0 * K_d + G)); /* -, Cheng (B.8) */ in cryerZeros()
1574 PetscReal nu_u = (3.0 * K_u - 2.0 * G) / (2.0 * (3.0 * K_u + G)); /* -, Cheng (B.9) */ in cryerZeros()