Lines Matching refs:PetscScalar
47 PetscScalar mu; /* shear modulus */
48 PetscScalar K_u; /* undrained bulk modulus */
49 PetscScalar alpha; /* Biot effective stress coefficient */
50 PetscScalar M; /* Biot modulus */
51 PetscScalar k; /* (isotropic) permeability */
52 PetscScalar mu_f; /* fluid dynamic viscosity */
53 PetscScalar P_0; /* magnitude of vertical stress */
71 …e zero(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in zero()
113 …atic_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in quadratic_u()
121 …ar_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in linear_eps()
127 …near_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in linear_linear_p()
137 …ar_p_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in linear_linear_p_t()
147 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_quadratic_linear_u()
161 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_quadratic_linear_p()
209 …trig_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in linear_trig_p()
219 …ig_p_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in linear_trig_p_t()
229 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_quadratic_trig_u()
243 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_quadratic_trig_p()
304 …trig_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in trig_u()
312 …ig_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in trig_eps()
322 …near_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in trig_linear_p()
332 …ar_p_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in trig_linear_p_t()
342 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_trig_linear_u()
356 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_trig_linear_p()
373 …essure(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_drainage_pressure()
380 PetscScalar alpha = param->alpha; /* - */ in terzaghi_drainage_pressure()
381 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_drainage_pressure()
382 PetscScalar M = param->M; /* Pa */ in terzaghi_drainage_pressure()
383 PetscScalar G = param->mu; /* Pa */ in terzaghi_drainage_pressure()
384 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_drainage_pressure()
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()
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()
396 …tial_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_initial_u()
403 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_initial_u()
404 PetscScalar G = param->mu; /* Pa */ in terzaghi_initial_u()
405 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_initial_u()
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()
416 …al_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_initial_eps()
423 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_initial_eps()
424 PetscScalar G = param->mu; /* Pa */ in terzaghi_initial_eps()
425 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_initial_eps()
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()
433 …i_2d_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_u()
442 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_u()
443 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_u()
444 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_u()
445 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_u()
446 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_u()
447 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_u()
451 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_u()
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()
455 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_u()
459 PetscScalar F2 = 0.0; in terzaghi_2d_u()
470 …2d_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_eps()
479 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_eps()
480 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_eps()
481 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_eps()
482 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_eps()
483 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_eps()
484 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_eps()
488 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_eps()
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()
492 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_eps()
496 PetscScalar F2_z = 0.0; in terzaghi_2d_eps()
507 …i_2d_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_p()
516 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p()
517 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_p()
518 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_p()
519 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p()
520 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_p()
521 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_p()
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()
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()
528 … PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_p()
532 PetscScalar F1 = 0.0; in terzaghi_2d_p()
544 …2d_u_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_u_t()
554 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_u_t()
555 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_u_t()
556 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_u_t()
557 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_u_t()
558 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_u_t()
559 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_u_t()
563 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_u_t()
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()
567 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_u_t()
571 PetscScalar F2_t = 0.0; in terzaghi_2d_u_t()
582 …_eps_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_eps_t()
591 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_eps_t()
592 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_eps_t()
593 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_eps_t()
594 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_eps_t()
595 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_eps_t()
596 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_eps_t()
600 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in terzaghi_2d_eps_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()
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()
604 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_eps_t()
608 PetscScalar F2_zt = 0.0; in terzaghi_2d_eps_t()
618 …2d_p_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in terzaghi_2d_p_t()
625 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p_t()
626 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_p_t()
627 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_p_t()
628 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p_t()
629 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_p_t()
630 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_p_t()
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()
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()
636 … PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_p_t()
640 PetscScalar alpha = param->alpha; /* - */ in terzaghi_2d_p_t()
641 PetscScalar K_u = param->K_u; /* Pa */ in terzaghi_2d_p_t()
642 PetscScalar M = param->M; /* Pa */ in terzaghi_2d_p_t()
643 PetscScalar G = param->mu; /* Pa */ in terzaghi_2d_p_t()
644 PetscScalar P_0 = param->P_0; /* Pa */ in terzaghi_2d_p_t()
645 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in terzaghi_2d_p_t()
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()
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()
652 … PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in terzaghi_2d_p_t()
656 PetscScalar F1_t = 0.0; in terzaghi_2d_p_t()
669 …essure(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_drainage_pressure()
676 PetscScalar alpha = param->alpha; /* - */ in mandel_drainage_pressure()
677 PetscScalar K_u = param->K_u; /* Pa */ in mandel_drainage_pressure()
678 PetscScalar M = param->M; /* Pa */ in mandel_drainage_pressure()
679 PetscScalar G = param->mu; /* Pa */ in mandel_drainage_pressure()
680 PetscScalar P_0 = param->P_0; /* Pa */ in mandel_drainage_pressure()
681 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in mandel_drainage_pressure()
685 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_drainage_pressure()
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()
687 …PetscScalar B = alpha * M / K_u; /* -, Cheng (B.12) */ 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()
689 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in mandel_drainage_pressure()
691 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u)); in mandel_drainage_pressure()
707 …tial_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_initial_u()
714 PetscScalar alpha = param->alpha; /* - */ in mandel_initial_u()
715 PetscScalar K_u = param->K_u; /* Pa */ in mandel_initial_u()
716 PetscScalar M = param->M; /* Pa */ in mandel_initial_u()
717 PetscScalar G = param->mu; /* Pa */ in mandel_initial_u()
718 PetscScalar P_0 = param->P_0; /* Pa */ in mandel_initial_u()
719 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in mandel_initial_u()
723 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_initial_u()
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()
742 …al_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_initial_eps()
749 PetscScalar alpha = param->alpha; /* - */ in mandel_initial_eps()
750 PetscScalar K_u = param->K_u; /* Pa */ in mandel_initial_eps()
751 PetscScalar M = param->M; /* Pa */ in mandel_initial_eps()
752 PetscScalar G = param->mu; /* Pa */ in mandel_initial_eps()
753 PetscScalar P_0 = param->P_0; /* Pa */ in mandel_initial_eps()
754 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in mandel_initial_eps()
758 … PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) */ in mandel_initial_eps()
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()
781 …l_2d_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_u()
792 PetscScalar alpha = param->alpha; in mandel_2d_u()
793 PetscScalar K_u = param->K_u; in mandel_2d_u()
794 PetscScalar M = param->M; in mandel_2d_u()
795 PetscScalar G = param->mu; in mandel_2d_u()
796 PetscScalar k = param->k; in mandel_2d_u()
797 PetscScalar mu_f = param->mu_f; in mandel_2d_u()
798 PetscScalar F = param->P_0; in mandel_2d_u()
800 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_u()
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()
803 PetscScalar kappa = k / mu_f; in mandel_2d_u()
808 PetscScalar A_x = 0.0; in mandel_2d_u()
809 PetscScalar B_x = 0.0; in mandel_2d_u()
824 …2d_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_eps()
835 PetscScalar alpha = param->alpha; in mandel_2d_eps()
836 PetscScalar K_u = param->K_u; in mandel_2d_eps()
837 PetscScalar M = param->M; in mandel_2d_eps()
838 PetscScalar G = param->mu; in mandel_2d_eps()
839 PetscScalar k = param->k; in mandel_2d_eps()
840 PetscScalar mu_f = param->mu_f; in mandel_2d_eps()
841 PetscScalar F = param->P_0; in mandel_2d_eps()
843 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_eps()
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()
846 PetscScalar kappa = k / mu_f; in mandel_2d_eps()
874 …l_2d_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_p()
886 PetscScalar alpha = param->alpha; in mandel_2d_p()
887 PetscScalar K_u = param->K_u; in mandel_2d_p()
888 PetscScalar M = param->M; in mandel_2d_p()
889 PetscScalar G = param->mu; in mandel_2d_p()
890 PetscScalar k = param->k; in mandel_2d_p()
891 PetscScalar mu_f = param->mu_f; in mandel_2d_p()
892 PetscScalar F = param->P_0; in mandel_2d_p()
894 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_p()
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()
897 PetscScalar kappa = k / mu_f; in mandel_2d_p()
898 PetscScalar B = (alpha * M) / (K_d + alpha * alpha * M); in mandel_2d_p()
902 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u)); in mandel_2d_p()
918 …2d_u_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_u_t()
927 PetscScalar alpha = param->alpha; in mandel_2d_u_t()
928 PetscScalar K_u = param->K_u; in mandel_2d_u_t()
929 PetscScalar M = param->M; in mandel_2d_u_t()
930 PetscScalar G = param->mu; in mandel_2d_u_t()
931 PetscScalar F = param->P_0; in mandel_2d_u_t()
933 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_u_t()
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()
936 PetscScalar kappa = param->k / param->mu_f; in mandel_2d_u_t()
941 PetscScalar A_s_t = 0.0; in mandel_2d_u_t()
942 PetscScalar B_s_t = 0.0; in mandel_2d_u_t()
958 …_eps_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_eps_t()
967 PetscScalar alpha = param->alpha; in mandel_2d_eps_t()
968 PetscScalar K_u = param->K_u; in mandel_2d_eps_t()
969 PetscScalar M = param->M; in mandel_2d_eps_t()
970 PetscScalar G = param->mu; in mandel_2d_eps_t()
971 PetscScalar k = param->k; in mandel_2d_eps_t()
972 PetscScalar mu_f = param->mu_f; in mandel_2d_eps_t()
973 PetscScalar F = param->P_0; in mandel_2d_eps_t()
975 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_eps_t()
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()
978 PetscScalar kappa = k / mu_f; in mandel_2d_eps_t()
986 PetscScalar eps_As = 0.0; in mandel_2d_eps_t()
987 PetscScalar eps_Bs = 0.0; in mandel_2d_eps_t()
988 PetscScalar eps_Cs = 0.0; in mandel_2d_eps_t()
1003 …2d_p_t(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in mandel_2d_p_t()
1011 PetscScalar alpha = param->alpha; in mandel_2d_p_t()
1012 PetscScalar K_u = param->K_u; in mandel_2d_p_t()
1013 PetscScalar M = param->M; in mandel_2d_p_t()
1014 PetscScalar G = param->mu; in mandel_2d_p_t()
1015 PetscScalar F = param->P_0; in mandel_2d_p_t()
1017 PetscScalar K_d = K_u - alpha * alpha * M; in mandel_2d_p_t()
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()
1031 …essure(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_drainage_pressure()
1038 PetscScalar alpha = param->alpha; /* - */ in cryer_drainage_pressure()
1039 PetscScalar K_u = param->K_u; /* Pa */ in cryer_drainage_pressure()
1040 PetscScalar M = param->M; /* Pa */ in cryer_drainage_pressure()
1041 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_drainage_pressure()
1042 PetscScalar B = alpha * M / K_u; /* -, Cheng (B.12) */ in cryer_drainage_pressure()
1051 …tial_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_initial_u()
1058 PetscScalar K_u = param->K_u; /* Pa */ in cryer_initial_u()
1059 PetscScalar G = param->mu; /* Pa */ in cryer_initial_u()
1060 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_initial_u()
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()
1074 …al_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_initial_eps()
1081 PetscScalar K_u = param->K_u; /* Pa */ in cryer_initial_eps()
1082 PetscScalar G = param->mu; /* Pa */ in cryer_initial_eps()
1083 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_initial_eps()
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()
1099 …r_3d_u(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_3d_u()
1108 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_u()
1109 PetscScalar K_u = param->K_u; /* Pa */ in cryer_3d_u()
1110 PetscScalar M = param->M; /* Pa */ in cryer_3d_u()
1111 PetscScalar G = param->mu; /* Pa */ in cryer_3d_u()
1112 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_3d_u()
1113 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in cryer_3d_u()
1117 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) … in cryer_3d_u()
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()
1121 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) … 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()
1128 PetscScalar u_sc; in cryer_3d_u()
1149 …3d_eps(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_3d_eps()
1158 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_eps()
1159 PetscScalar K_u = param->K_u; /* Pa */ in cryer_3d_eps()
1160 PetscScalar M = param->M; /* Pa */ in cryer_3d_eps()
1161 PetscScalar G = param->mu; /* Pa */ in cryer_3d_eps()
1162 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_3d_eps()
1163 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in cryer_3d_eps()
1167 …PetscScalar K_d = K_u - alpha * alpha * M; /* Pa, Cheng (B.5) … in cryer_3d_eps()
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()
1171 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) … 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()
1201 …r_3d_p(PetscInt dim, PetscReal time, const PetscReal x[], PetscInt Nc, PetscScalar *u, PetscCtx ct… in cryer_3d_p()
1210 PetscScalar alpha = param->alpha; /* - */ in cryer_3d_p()
1211 PetscScalar K_u = param->K_u; /* Pa */ in cryer_3d_p()
1212 PetscScalar M = param->M; /* Pa */ in cryer_3d_p()
1213 PetscScalar G = param->mu; /* Pa */ in cryer_3d_p()
1214 PetscScalar P_0 = param->P_0; /* Pa */ in cryer_3d_p()
1216 PetscScalar kappa = param->k / param->mu_f; /* m^2 / (Pa s) */ in cryer_3d_p()
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()
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()
1224 …PetscScalar c = kappa / S; /* m^2 / s, Cheng (B.16) */ in cryer_3d_p()
1243 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_terzaghi_bd_u()
1253 …etscInt uOff[], const PetscInt uOff_x[], const PetscScalar u[], const PetscScalar u_t[], const Pet…
1254 …etscInt aOff[], const PetscInt aOff_x[], const PetscScalar a[], const PetscScalar a_t[], const Pet…
1255 …al x[], const PetscReal n[], PetscInt numConstants, const PetscScalar constants[], PetscScalar f0[…
1289 PetscScalar G = constants[0];
1290 PetscScalar K_u = constants[1];
1291 PetscScalar alpha = constants[2];
1292 PetscScalar M = constants[3];
1293 PetscScalar kappa = constants[4];
1294 PetscScalar F = constants[5];
1296 PetscScalar K_d = K_u - alpha*alpha*M;
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));
1301 PetscScalar B = (3.0 * (nu_u - nu)) / ( alpha * (1.0 - 2.0*nu) * (1.0 + nu_u));
1302 PetscScalar A1 = 3.0 / (B * (1.0 + nu_u));
1303 PetscScalar A2 = (alpha * (1.0 - 2.0*nu)) / (1.0 - nu);
1331 PetscScalar A_x = 0.0;
1332 PetscScalar B_x = 0.0;
1342 PetscScalar sigma_zz = -1.0*(F/aL) - ((2.0*F)/aL) * (A2/A1) * A_x + ((2.0*F)/aL) * B_x;
1352 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_cryer_bd_u()
1362 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_epsilon()
1371 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f0_p()
1382 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f1_u()
1401 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in f1_p()
1418 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g0_ee()
1424 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g0_pe()
1432 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g0_pp()
1440 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g1_eu()
1447 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g2_ue()
1460 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g2_up()
1469 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g3_uu()
1484 …PetscScalar u[], const PetscScalar u_t[], const PetscScalar u_x[], const PetscInt aOff[], const Pe… in g3_pp()
1524 PetscScalar alpha = param->alpha; in mandelZeros()
1525 PetscScalar K_u = param->K_u; in mandelZeros()
1526 PetscScalar M = param->M; in mandelZeros()
1527 PetscScalar G = param->mu; in mandelZeros()
1532 PetscScalar K_d = K_u - alpha * alpha * M; in mandelZeros()
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()
1656 PetscScalar K_d = p->K_u - p->alpha * p->alpha * p->M; in SetupParameters()
1657 PetscScalar nu_u = (3.0 * p->K_u - 2.0 * p->mu) / (2.0 * (3.0 * p->K_u + p->mu)); in SetupParameters()
1658 PetscScalar nu = (3.0 * K_d - 2.0 * p->mu) / (2.0 * (3.0 * K_d + p->mu)); in SetupParameters()
1659 PetscScalar S = (3.0 * p->K_u + 4.0 * p->mu) / (p->M * (3.0 * K_d + 4.0 * p->mu)); in SetupParameters()
1712 …PetscErrorCode (*exact[3])(PetscInt, PetscReal, const PetscReal[], PetscInt, PetscScalar *, void *… in SetupPrimalProblem()
1713 …PetscErrorCode (*exact_t[3])(PetscInt, PetscReal, const PetscReal[], PetscInt, PetscScalar *, void… in SetupPrimalProblem()
1898 PetscScalar constants[6]; in SetupPrimalProblem()
1963 …PetscErrorCode (*funcs[3])(PetscInt, PetscReal, const PetscReal[], PetscInt, PetscScalar *, void *… in SetInitialConditions()
2027 …PetscErrorCode (**exacts)(PetscInt, PetscReal, const PetscReal x[], PetscInt, PetscScalar *u, Pets… in SolutionMonitor()