Lines Matching refs:alpha
49 PetscScalar alpha; /* Biot effective stress coefficient */ member
151 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_linear_u() local
153 const PetscReal K_d = K_u - alpha * alpha * M; 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()
163 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_linear_p() local
169 f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0; in f0_quadratic_linear_p()
233 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_trig_u() local
235 const PetscReal K_d = K_u - alpha * alpha * M; 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()
245 const PetscReal alpha = PetscRealPart(constants[2]); in f0_quadratic_trig_p() local
252 f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0; in f0_quadratic_trig_p()
346 const PetscReal alpha = PetscRealPart(constants[2]); in f0_trig_linear_u() local
348 const PetscReal K_d = K_u - alpha * alpha * M; in f0_trig_linear_u()
352 …TSC_PI * x[d]) * (2. * G + lambda) + 2.0 * (G + lambda) - 2. * PETSC_PI * alpha * PetscSinReal(2. … in f0_trig_linear_u()
353 …SinReal(2. * PETSC_PI * x[dim - 1]) * (2. * G + lambda) - 2. * PETSC_PI * alpha * PetscSinReal(2. … in f0_trig_linear_u()
358 const PetscReal alpha = PetscRealPart(constants[2]); in f0_trig_linear_p() local
365 f0[0] += u_t ? alpha * u_t[uOff[1]] : 0.0; in f0_trig_linear_p()
380 PetscScalar alpha = param->alpha; /* - */ in terzaghi_drainage_pressure() local
385 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) … in terzaghi_drainage_pressure()
386 …PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) … in terzaghi_drainage_pressure()
442 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_u() local
451 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_u()
479 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_eps() local
488 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_eps()
516 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p() local
525 … PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_p()
526 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p()
534 …alpha * eta) / G) - S) <= 1.0e-10, PETSC_COMM_SELF, PETSC_ERR_PLIB, "S %g != check %g", (double)Pe… in terzaghi_2d_p()
554 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_u_t() local
563 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_u_t()
591 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_eps_t() local
600 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_eps_t()
625 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p_t() local
633 … PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_p_t()
634 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p_t()
640 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p_t() local
649 … PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_p_t()
650 … PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in terzaghi_2d_p_t()
658 …alpha * eta) / G) - S) <= 1.0e-10, PETSC_COMM_SELF, PETSC_ERR_PLIB, "S %g != check %g", (double)Pe… in terzaghi_2d_p_t()
676 PetscScalar alpha = param->alpha; /* - */ in mandel_drainage_pressure() local
685 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_drainage_pressure()
687 …PetscScalar B = alpha * M / K_u; /* -, Cheng (B.12) */ in mandel_drainage_pressure()
714 PetscScalar alpha = param->alpha; /* - */ in mandel_initial_u() local
723 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_initial_u()
749 PetscScalar alpha = param->alpha; /* - */ in mandel_initial_eps() local
758 … PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_initial_eps()
792 PetscScalar alpha = param->alpha; in mandel_2d_u() local
800 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_u()
805 … = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha * (1.0 - 2.0 *… in mandel_2d_u()
835 PetscScalar alpha = param->alpha; in mandel_2d_eps() local
843 PetscScalar K_d = K_u - alpha * alpha * M; 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()
886 PetscScalar alpha = param->alpha; in mandel_2d_p() local
894 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_p()
898 PetscScalar B = (alpha * M) / (K_d + alpha * alpha * M); in mandel_2d_p()
901 …PetscReal c = PetscRealPart(((2.0 * kappa * G) * (1.0 - nu) * (nu_u - nu)) / (alpha * alpha * (… in mandel_2d_p()
927 PetscScalar alpha = param->alpha; in mandel_2d_u_t() local
933 PetscScalar K_d = K_u - alpha * alpha * M; 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()
967 PetscScalar alpha = param->alpha; in mandel_2d_eps_t() local
975 PetscScalar K_d = K_u - alpha * alpha * M; 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()
1011 PetscScalar alpha = param->alpha; in mandel_2d_p_t() local
1017 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_p_t()
1025 u[0] = ((2.0 * F * (-2.0 * nu + 3.0 * nu_u)) / (3.0 * a * alpha * (1.0 - 2.0 * nu))); in mandel_2d_p_t()
1038 PetscScalar alpha = param->alpha; /* - */ in cryer_drainage_pressure() local
1042 PetscScalar B = alpha * M / K_u; /* -, Cheng (B.12) */ in cryer_drainage_pressure()
1108 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_u() local
1117 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) … in cryer_3d_u()
1158 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_eps() local
1167 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) … in cryer_3d_eps()
1210 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_p() local
1219 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in cryer_3d_p()
1220 …PetscScalar eta = (3.0 * alpha * G) / (3.0 * K_d + 4.0 * G); /* -, Cheng (B.11) */ in cryer_3d_p()
1291 PetscScalar alpha = constants[2];
1296 PetscScalar K_d = K_u - alpha*alpha*M;
1300 …PetscReal c = PetscRealPart(((2.0*kappa*G) * (1.0 - nu) * (nu_u - nu)) / (alpha*alpha * (1.0 - 2…
1301 PetscScalar B = (3.0 * (nu_u - nu)) / ( alpha * (1.0 - 2.0*nu) * (1.0 + nu_u));
1303 PetscScalar A2 = (alpha * (1.0 - 2.0*nu)) / (1.0 - nu);
1373 const PetscReal alpha = PetscRealPart(constants[2]); in f0_p() local
1376 f0[0] += alpha * u_t[uOff[1]]; in f0_p()
1387 const PetscReal alpha = PetscRealPart(constants[2]); in f1_u() local
1389 const PetscReal K_d = K_u - alpha * alpha * M; in f1_u()
1396 f1[c * dim + c] += alpha * u[uOff[2]]; in f1_u()
1426 const PetscReal alpha = PetscRealPart(constants[2]); in g0_pe() local
1428 g0[0] = u_tShift * alpha; in g0_pe()
1451 const PetscReal alpha = PetscRealPart(constants[2]); in g2_ue() local
1453 const PetscReal K_d = K_u - alpha * alpha * M; in g2_ue()
1462 const PetscReal alpha = PetscRealPart(constants[2]); in g2_up() local
1465 for (d = 0; d < dim; ++d) g2[d * dim + d] += alpha; in g2_up()
1524 PetscScalar alpha = param->alpha; in mandelZeros() local
1532 PetscScalar K_d = K_u - alpha * alpha * M; in mandelZeros()
1566 PetscReal alpha = PetscRealPart(param->alpha); /* - */ in cryerZeros() local
1572 PetscReal K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in cryerZeros()
1621 …PetscCall(PetscBagRegisterScalar(bag, &p->alpha, 0.6, "alpha", "Biot Effective Stress Coefficient,… in SetupParameters()
1630 …PetscCall(PetscBagRegisterScalar(bag, &p->alpha, 0.6, "alpha", "Biot Effective Stress Coefficient,… in SetupParameters()
1639 …PetscCall(PetscBagRegisterScalar(bag, &p->alpha, 0.6, "alpha", "Biot Effective Stress Coefficient,… in SetupParameters()
1648 …PetscCall(PetscBagRegisterScalar(bag, &p->alpha, 1.0, "alpha", "Biot Effective Stress Coefficient,… in SetupParameters()
1656 PetscScalar K_d = p->K_u - p->alpha * p->alpha * p->M; in SetupParameters()
1901 constants[2] = param->alpha; /* Biot effective stress coefficient, - */ in SetupPrimalProblem()