From ada9989e96923e1c84445228e207c764c9c5e894 Mon Sep 17 00:00:00 2001 From: Bob Myhill Date: Wed, 31 Jul 2024 16:45:57 +0100 Subject: [PATCH] update tests --- tests/composite_viscous_outputs.cc | 6 +- tests/composite_viscous_outputs/screen-output | 24 +++---- tests/composite_viscous_outputs_isostress.cc | 2 +- tests/composite_viscous_outputs_limited.cc | 9 +-- .../screen-output | 24 +++---- tests/composite_viscous_outputs_no_peierls.cc | 6 +- .../screen-output | 24 +++---- tests/composite_viscous_outputs_phases.cc | 9 +-- .../screen-output | 68 +++++++++---------- 9 files changed, 87 insertions(+), 85 deletions(-) diff --git a/tests/composite_viscous_outputs.cc b/tests/composite_viscous_outputs.cc index 6f8f65301af..8fc7cae9dd3 100644 --- a/tests/composite_viscous_outputs.cc +++ b/tests/composite_viscous_outputs.cc @@ -107,7 +107,7 @@ void f(const aspect::SimulatorAccess &simulator_access, strain_rate[2][1] = 0.; strain_rate[2][2] = 0.; - std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max)" << std::endl; + std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max)" << std::endl; // Loop through strain rates, tracking whether there is a discrepancy in // the decomposed strain rates. @@ -120,7 +120,7 @@ void f(const aspect::SimulatorAccess &simulator_access, double disl_stress; double prls_stress; double drpr_stress; - std::vector partial_strain_rates(5, 0.); + std::vector partial_strain_rates(6, 0.); for (unsigned int i=0; i <= 10; i++) { @@ -134,7 +134,7 @@ void f(const aspect::SimulatorAccess &simulator_access, // of the max viscosity dashpot from the total strain rate // The creep stress is then calculated by subtracting the stress running // through the strain rate limiter from the total stress - creep_strain_rate = total_strain_rate - partial_strain_rates[4]; + creep_strain_rate = total_strain_rate - partial_strain_rates[4] - partial_strain_rates[5]; creep_stress = 2.*(viscosity*total_strain_rate - lim_visc*creep_strain_rate); // Print the output diff --git a/tests/composite_viscous_outputs/screen-output b/tests/composite_viscous_outputs/screen-output index 77b2499d232..25b78b8b328 100644 --- a/tests/composite_viscous_outputs/screen-output +++ b/tests/composite_viscous_outputs/screen-output @@ -2,18 +2,18 @@ Loading shared library <./libcomposite_viscous_outputs.debug.so> * Connecting signals -temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max) -1000 3.45602e+19 6.89205e+08 1e-11 1.29846e-06 0.00363514 0.0654305 0.930933 3.45602e-09 -1100 2.98863e+19 5.95725e+08 1e-11 7.23285e-05 0.835897 0.163395 0.0006365 2.98863e-09 -1200 7.70568e+18 1.52114e+08 1e-11 0.000594407 0.999399 6.44893e-06 1.44515e-33 7.70568e-10 -1300 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 3.19443e-09 1.32277e-59 2.39282e-10 -1400 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 1.26589e-11 5.56077e-82 9.18169e-11 -1500 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 2.09885e-13 1.4634e-101 4.32084e-11 -1600 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 8.74429e-15 3.2006e-119 2.47246e-11 -1700 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 5.42942e-16 3.38184e-136 1.67378e-11 -1800 1.28447e+17 568948 1e-11 0.749975 0.250025 2.2682e-17 6.38432e-155 1.28447e-11 -1900 1.09563e+17 191256 1e-11 0.962698 0.0373021 2.59166e-19 1.35594e-178 1.09563e-11 -2000 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 2.26225e-21 4.97674e-204 1.02965e-11 +temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max) +1000 3.45602e+19 6.89205e+08 1e-11 1.29846e-06 0.00363514 0.0654305 0.930933 0 3.45602e-09 +1100 2.98863e+19 5.95725e+08 1e-11 7.23285e-05 0.835897 0.163395 0.0006365 0 2.98863e-09 +1200 7.70568e+18 1.52114e+08 1e-11 0.000594407 0.999399 6.44893e-06 1.44515e-33 0 7.70568e-10 +1300 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 3.19443e-09 1.32277e-59 0 2.39282e-10 +1400 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 1.26589e-11 5.56077e-82 0 9.18169e-11 +1500 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 2.09885e-13 1.4634e-101 0 4.32084e-11 +1600 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 8.74429e-15 3.2006e-119 0 2.47246e-11 +1700 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 5.42942e-16 3.38184e-136 0 1.67378e-11 +1800 1.28447e+17 568948 1e-11 0.749975 0.250025 2.2682e-17 6.38432e-155 0 1.28447e-11 +1900 1.09563e+17 191256 1e-11 0.962698 0.0373021 2.59166e-19 1.35594e-178 0 1.09563e-11 +2000 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 2.26225e-21 4.97674e-204 0 1.02965e-11 OK Number of active cells: 100 (on 1 levels) Number of degrees of freedom: 5,534 (3,969+242+1,323) diff --git a/tests/composite_viscous_outputs_isostress.cc b/tests/composite_viscous_outputs_isostress.cc index a8073a5e2c7..608280d62f5 100644 --- a/tests/composite_viscous_outputs_isostress.cc +++ b/tests/composite_viscous_outputs_isostress.cc @@ -136,7 +136,7 @@ void f(const aspect::SimulatorAccess &simulator_access, // of the max viscosity dashpot from the total strain rate // The creep stress is then calculated by subtracting the stress running // through the strain rate limiter from the total stress - creep_strain_rate = total_strain_rate - partial_strain_rates[5]; + creep_strain_rate = total_strain_rate - partial_strain_rates[4] - partial_strain_rates[5]; creep_stress = 2.*(viscosity*total_strain_rate - lim_visc*creep_strain_rate); // Print the output diff --git a/tests/composite_viscous_outputs_limited.cc b/tests/composite_viscous_outputs_limited.cc index d7ec481b536..e6dff709bcb 100644 --- a/tests/composite_viscous_outputs_limited.cc +++ b/tests/composite_viscous_outputs_limited.cc @@ -109,6 +109,7 @@ void f(const aspect::SimulatorAccess &simulator_access, double temperature; const double pressure = 1.e9; const double grain_size = 1.e-3; + const double inverse_kelvin_viscosity = 0.; SymmetricTensor<2,dim> strain_rate; strain_rate[0][0] = -1e-11; strain_rate[0][1] = 0.; @@ -117,7 +118,7 @@ void f(const aspect::SimulatorAccess &simulator_access, strain_rate[2][1] = 0.; strain_rate[2][2] = 0.; - std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max)" << std::endl; + std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max)" << std::endl; // Loop through strain rates, tracking whether there is a discrepancy in // the decomposed strain rates. @@ -130,21 +131,21 @@ void f(const aspect::SimulatorAccess &simulator_access, double disl_stress; double prls_stress; double drpr_stress; - std::vector partial_strain_rates(5, 0.); + std::vector partial_strain_rates(6, 0.); for (unsigned int i=0; i <= 10; i++) { temperature = 1000. + i*100.; // Compute the viscosity - viscosity = composite_creep->compute_viscosity(pressure, temperature, grain_size, volume_fractions, strain_rate, partial_strain_rates); + viscosity = composite_creep->compute_viscosity(pressure, temperature, grain_size, volume_fractions, strain_rate, inverse_kelvin_viscosity, partial_strain_rates); total_strain_rate = std::accumulate(partial_strain_rates.begin(), partial_strain_rates.end(), 0.); // The creep strain rate is calculated by subtracting the strain rate // of the max viscosity dashpot from the total strain rate // The creep stress is then calculated by subtracting the stress running // through the strain rate limiter from the total stress - creep_strain_rate = total_strain_rate - partial_strain_rates[4]; + creep_strain_rate = total_strain_rate - partial_strain_rates[4] - partial_strain_rates[5]; creep_stress = 2.*(viscosity*total_strain_rate - lim_visc*creep_strain_rate); // Print the output diff --git a/tests/composite_viscous_outputs_limited/screen-output b/tests/composite_viscous_outputs_limited/screen-output index 3c63ad3cc86..1415da04a77 100644 --- a/tests/composite_viscous_outputs_limited/screen-output +++ b/tests/composite_viscous_outputs_limited/screen-output @@ -2,18 +2,18 @@ Loading shared library <./libcomposite_viscous_outputs_limited.debug.so> * Connecting signals -temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max) -1000 3.99999e+18 7.99998e+07 1e-11 1.5072e-07 1.93692e-06 1.58471e-11 1.60671e-47 0.999998 -1100 3.997e+18 7.99201e+07 1e-11 9.7033e-06 0.000739244 7.95399e-10 1.52864e-47 0.999251 -1200 3.70523e+18 7.21395e+07 1e-11 0.000281896 0.0734102 8.42438e-09 9.12346e-50 0.926308 -1300 2.33707e+18 3.56552e+07 1e-11 0.00262872 0.413104 3.88781e-10 4.546e-65 0.584267 -1400 1.5335e+18 1.42267e+07 1e-11 0.0130073 0.603618 4.18203e-12 5.08959e-85 0.383375 -1500 1.22332e+18 5.95532e+06 1e-11 0.0482676 0.645901 9.24051e-14 6.26253e-104 0.305831 -1600 1.09939e+18 2.65043e+06 1e-11 0.144969 0.580183 4.10546e-15 1.64922e-121 0.274848 -1700 1.04457e+18 1.1886e+06 1e-11 0.350471 0.388386 2.29166e-16 6.35863e-139 0.261143 -1800 1.01759e+18 469097 1e-11 0.618354 0.127248 6.34559e-18 4.11763e-159 0.254398 -1900 1.00547e+18 145858 1e-11 0.734184 0.0144489 4.62654e-20 1.77001e-184 0.251367 -2000 1.00167e+18 44514.1 1e-11 0.748313 0.00126937 3.89551e-22 2.99484e-210 0.250417 +temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max) +1000 3.99999e+18 7.99998e+07 1e-11 1.5072e-07 1.93692e-06 1.58471e-11 1.60671e-47 0 0.999998 +1100 3.997e+18 7.99201e+07 1e-11 9.7033e-06 0.000739244 7.95399e-10 1.52864e-47 0 0.999251 +1200 3.70523e+18 7.21395e+07 1e-11 0.000281896 0.0734102 8.42438e-09 9.12346e-50 0 0.926308 +1300 2.33707e+18 3.56552e+07 1e-11 0.00262872 0.413104 3.88781e-10 4.546e-65 0 0.584267 +1400 1.5335e+18 1.42267e+07 1e-11 0.0130073 0.603618 4.18203e-12 5.08959e-85 0 0.383375 +1500 1.22332e+18 5.95532e+06 1e-11 0.0482676 0.645901 9.24051e-14 6.26253e-104 0 0.305831 +1600 1.09939e+18 2.65043e+06 1e-11 0.144969 0.580183 4.10546e-15 1.64922e-121 0 0.274848 +1700 1.04457e+18 1.1886e+06 1e-11 0.350471 0.388386 2.29166e-16 6.35863e-139 0 0.261143 +1800 1.01759e+18 469097 1e-11 0.618354 0.127248 6.34559e-18 4.11763e-159 0 0.254398 +1900 1.00547e+18 145858 1e-11 0.734184 0.0144489 4.62654e-20 1.77001e-184 0 0.251367 +2000 1.00167e+18 44514.1 1e-11 0.748313 0.00126937 3.89551e-22 2.99484e-210 0 0.250417 OK Number of active cells: 100 (on 1 levels) Number of degrees of freedom: 5,534 (3,969+242+1,323) diff --git a/tests/composite_viscous_outputs_no_peierls.cc b/tests/composite_viscous_outputs_no_peierls.cc index e9bc6b0a158..5c139f5c610 100644 --- a/tests/composite_viscous_outputs_no_peierls.cc +++ b/tests/composite_viscous_outputs_no_peierls.cc @@ -101,7 +101,7 @@ void f(const aspect::SimulatorAccess &simulator_access, strain_rate[2][1] = 0.; strain_rate[2][2] = 0.; - std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max)" << std::endl; + std::cout << "temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max)" << std::endl; // Loop through strain rates, tracking whether there is a discrepancy in // the decomposed strain rates. @@ -113,7 +113,7 @@ void f(const aspect::SimulatorAccess &simulator_access, double diff_stress; double disl_stress; double drpr_stress; - std::vector partial_strain_rates(5, 0.); + std::vector partial_strain_rates(6, 0.); for (unsigned int i=0; i <= 10; i++) { @@ -127,7 +127,7 @@ void f(const aspect::SimulatorAccess &simulator_access, // of the max viscosity dashpot from the total strain rate // The creep stress is then calculated by subtracting the stress running // through the strain rate limiter from the total stress - creep_strain_rate = total_strain_rate - partial_strain_rates[4]; + creep_strain_rate = total_strain_rate - partial_strain_rates[4] - partial_strain_rates[5]; creep_stress = 2.*(viscosity*total_strain_rate - lim_visc*creep_strain_rate); // Print the output diff --git a/tests/composite_viscous_outputs_no_peierls/screen-output b/tests/composite_viscous_outputs_no_peierls/screen-output index e1c4155b475..ab3c2caa4bb 100644 --- a/tests/composite_viscous_outputs_no_peierls/screen-output +++ b/tests/composite_viscous_outputs_no_peierls/screen-output @@ -2,18 +2,18 @@ Loading shared library <./libcomposite_viscous_outputs_no_peierls.debug.so> * Connecting signals -temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max) -1000 3.46071e+19 6.90142e+08 1e-11 1.30023e-06 0.00365246 0 0.996346 3.46071e-09 -1100 3.13843e+19 6.25686e+08 1e-11 7.59662e-05 0.992522 0 0.00740179 3.13843e-09 -1200 7.7057e+18 1.52114e+08 1e-11 0.000594408 0.999406 0 1.44528e-33 7.7057e-10 -1300 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 0 1.32277e-59 2.39282e-10 -1400 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 0 5.56077e-82 9.18169e-11 -1500 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 0 1.4634e-101 4.32084e-11 -1600 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 0 3.2006e-119 2.47246e-11 -1700 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 0 3.38184e-136 1.67378e-11 -1800 1.28447e+17 568948 1e-11 0.749975 0.250025 0 6.38432e-155 1.28447e-11 -1900 1.09563e+17 191256 1e-11 0.962698 0.0373021 0 1.35594e-178 1.09563e-11 -2000 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 0 4.97674e-204 1.02965e-11 +temperature (K) eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max) +1000 3.46071e+19 6.90142e+08 1e-11 1.30023e-06 0.00365246 0 0.996346 0 3.46071e-09 +1100 3.13843e+19 6.25686e+08 1e-11 7.59662e-05 0.992522 0 0.00740179 0 3.13843e-09 +1200 7.7057e+18 1.52114e+08 1e-11 0.000594408 0.999406 0 1.44528e-33 0 7.7057e-10 +1300 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 0 1.32277e-59 0 2.39282e-10 +1400 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 0 5.56077e-82 0 9.18169e-11 +1500 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 0 1.4634e-101 0 4.32084e-11 +1600 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 0 3.2006e-119 0 2.47246e-11 +1700 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 0 3.38184e-136 0 1.67378e-11 +1800 1.28447e+17 568948 1e-11 0.749975 0.250025 0 6.38432e-155 0 1.28447e-11 +1900 1.09563e+17 191256 1e-11 0.962698 0.0373021 0 1.35594e-178 0 1.09563e-11 +2000 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 0 4.97674e-204 0 1.02965e-11 OK Number of active cells: 100 (on 1 levels) Number of degrees of freedom: 5,534 (3,969+242+1,323) diff --git a/tests/composite_viscous_outputs_phases.cc b/tests/composite_viscous_outputs_phases.cc index 3650747cfab..e98f6e1ff03 100644 --- a/tests/composite_viscous_outputs_phases.cc +++ b/tests/composite_viscous_outputs_phases.cc @@ -125,6 +125,7 @@ void f(const aspect::SimulatorAccess &simulator_access, double temperature; double pressure; const double grain_size = 1.e-3; + const double inverse_kelvin_viscosity = 0.; SymmetricTensor<2,dim> strain_rate; strain_rate[0][0] = -1e-11; strain_rate[0][1] = 0.; @@ -133,7 +134,7 @@ void f(const aspect::SimulatorAccess &simulator_access, strain_rate[2][1] = 0.; strain_rate[2][2] = 0.; - std::cout << "temperature (K) phase transition progress eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max)" << std::endl; + std::cout << "temperature (K) phase transition progress eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max)" << std::endl; // Loop through strain rates, tracking whether there is a discrepancy in // the decomposed strain rates. @@ -146,7 +147,7 @@ void f(const aspect::SimulatorAccess &simulator_access, double disl_stress; double prls_stress; double drpr_stress; - std::vector partial_strain_rates(5, 0.); + std::vector partial_strain_rates(6, 0.); for (unsigned int i=0; i <= 2; i++) { @@ -176,14 +177,14 @@ void f(const aspect::SimulatorAccess &simulator_access, Rheology::DruckerPragerParameters p = drucker_prager_power->compute_drucker_prager_parameters(composition, phase_function_values, n_phase_transitions_per_composition); // Compute the viscosity - viscosity = composite_creep->compute_viscosity(pressure, temperature, grain_size, volume_fractions, strain_rate, partial_strain_rates, phase_function_values, n_phase_transitions_per_composition); + viscosity = composite_creep->compute_viscosity(pressure, temperature, grain_size, volume_fractions, strain_rate, inverse_kelvin_viscosity, partial_strain_rates, phase_function_values, n_phase_transitions_per_composition); total_strain_rate = std::accumulate(partial_strain_rates.begin(), partial_strain_rates.end(), 0.); // The creep strain rate is calculated by subtracting the strain rate // of the max viscosity dashpot from the total strain rate // The creep stress is then calculated by subtracting the stress running // through the strain rate limiter from the total stress - creep_strain_rate = total_strain_rate - partial_strain_rates[4]; + creep_strain_rate = total_strain_rate - partial_strain_rates[4] - partial_strain_rates[5]; creep_stress = 2.*(viscosity*total_strain_rate - lim_visc*creep_strain_rate); // Print the output diff --git a/tests/composite_viscous_outputs_phases/screen-output b/tests/composite_viscous_outputs_phases/screen-output index 6b0f9cf0338..6d9def9739b 100644 --- a/tests/composite_viscous_outputs_phases/screen-output +++ b/tests/composite_viscous_outputs_phases/screen-output @@ -2,43 +2,43 @@ Loading shared library <./libcomposite_viscous_outputs_phases.debug.so> * Connecting signals -temperature (K) phase transition progress eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, max) +temperature (K) phase transition progress eta (Pas) creep stress (Pa) edot_ii (/s) edot_ii fractions (diff, disl, prls, drpr, kel, max) pressure: 1 GPa -1000 0.0179862 4.49858e+19 8.97717e+08 1e-11 1.6913e-06 0.00916831 0.99083 5.3662e-13 4.49858e-09 -1100 0.0179862 2.98905e+19 5.9581e+08 1e-11 7.23388e-05 0.836312 0.163616 6.73098e-22 2.98905e-09 -1200 0.0179862 7.70568e+18 1.52114e+08 1e-11 0.000594407 0.999399 6.44893e-06 1.51744e-51 7.70568e-10 -1300 0.0179862 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 3.19443e-09 1.38894e-77 2.39282e-10 -1400 0.0179862 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 1.26589e-11 5.83892e-100 9.18169e-11 -1500 0.0179862 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 2.09885e-13 1.5366e-119 4.32084e-11 -1600 0.0179862 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 8.74429e-15 3.3607e-137 2.47246e-11 -1700 0.0179862 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 5.42942e-16 3.551e-154 1.67378e-11 -1800 0.0179862 1.28447e+17 568948 1e-11 0.749975 0.250025 2.2682e-17 6.70367e-173 1.28447e-11 -1900 0.0179862 1.09563e+17 191256 1e-11 0.962698 0.0373021 2.59166e-19 1.42377e-196 1.09563e-11 -2000 0.0179862 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 2.26225e-21 5.22568e-222 1.02965e-11 +1000 0.0179862 4.49858e+19 8.97717e+08 1e-11 1.6913e-06 0.00916831 0.99083 5.3662e-13 0 4.49858e-09 +1100 0.0179862 2.98905e+19 5.9581e+08 1e-11 7.23388e-05 0.836312 0.163616 6.73098e-22 0 2.98905e-09 +1200 0.0179862 7.70568e+18 1.52114e+08 1e-11 0.000594407 0.999399 6.44893e-06 1.51744e-51 0 7.70568e-10 +1300 0.0179862 2.39282e+18 4.58564e+07 1e-11 0.00338081 0.996619 3.19443e-09 1.38894e-77 0 2.39282e-10 +1400 0.0179862 9.18168e+17 1.63634e+07 1e-11 0.0149609 0.985039 1.26589e-11 5.83892e-100 0 9.18169e-11 +1500 0.0179862 4.32083e+17 6.64167e+06 1e-11 0.0538305 0.94617 2.09885e-13 1.5366e-119 0 4.32084e-11 +1600 0.0179862 2.47246e+17 2.94493e+06 1e-11 0.161077 0.838923 8.74429e-15 3.3607e-137 0 2.47246e-11 +1700 0.0179862 1.67379e+17 1.34758e+06 1e-11 0.397345 0.602655 5.42942e-16 3.551e-154 0 1.67378e-11 +1800 0.0179862 1.28447e+17 568948 1e-11 0.749975 0.250025 2.2682e-17 6.70367e-173 0 1.28447e-11 +1900 0.0179862 1.09563e+17 191256 1e-11 0.962698 0.0373021 2.59166e-19 1.42377e-196 0 1.09563e-11 +2000 0.0179862 1.02964e+17 59280.1 1e-11 0.99654 0.00345962 2.26225e-21 5.22568e-222 0 1.02965e-11 pressure: 3 GPa -1000 0.5 5.3691e+19 1.07182e+09 1e-11 4.76869e-07 0.000587744 0.99515 0.00426168 5.3691e-09 -1100 0.5 4.29097e+19 8.56193e+08 1e-11 2.79907e-05 0.139287 0.860685 5.65005e-08 4.29097e-09 -1200 0.5 1.70541e+19 3.39081e+08 1e-11 0.000398001 0.998688 0.000914421 4.35181e-28 1.70541e-09 -1300 0.5 4.9077e+18 9.6154e+07 1e-11 0.00233578 0.997664 1.07007e-07 1.86999e-55 4.9077e-10 -1400 0.5 1.72879e+18 3.25758e+07 1e-11 0.0106234 0.989377 1.4219e-10 5.86483e-79 1.72879e-10 -1500 0.5 7.33471e+17 1.26694e+07 1e-11 0.0392318 0.960768 1.03038e-12 1.82512e-99 7.3347e-11 -1600 0.5 3.72285e+17 5.4457e+06 1e-11 0.120854 0.879146 2.33681e-14 8.43651e-118 3.72285e-11 -1700 0.5 2.2329e+17 2.46581e+06 1e-11 0.311075 0.688925 9.90838e-16 5.26533e-135 2.2329e-11 -1800 0.5 1.53962e+17 1.07924e+06 1e-11 0.638071 0.361929 3.94829e-17 6.01736e-153 1.53962e-11 -1900 0.5 1.19579e+17 391585 1e-11 0.922158 0.0778423 5.25907e-19 5.8183e-175 1.19579e-11 -2000 0.5 1.06072e+17 121442 1e-11 0.992096 0.0079044 3.18167e-21 2.19751e-200 1.06071e-11 +1000 0.5 5.3691e+19 1.07182e+09 1e-11 4.76869e-07 0.000587744 0.99515 0.00426168 0 5.3691e-09 +1100 0.5 4.29097e+19 8.56193e+08 1e-11 2.79907e-05 0.139287 0.860685 5.65005e-08 0 4.29097e-09 +1200 0.5 1.70541e+19 3.39081e+08 1e-11 0.000398001 0.998688 0.000914421 4.35181e-28 0 1.70541e-09 +1300 0.5 4.9077e+18 9.6154e+07 1e-11 0.00233578 0.997664 1.07007e-07 1.86999e-55 0 4.9077e-10 +1400 0.5 1.72879e+18 3.25758e+07 1e-11 0.0106234 0.989377 1.4219e-10 5.86483e-79 0 1.72879e-10 +1500 0.5 7.33471e+17 1.26694e+07 1e-11 0.0392318 0.960768 1.03038e-12 1.82512e-99 0 7.3347e-11 +1600 0.5 3.72285e+17 5.4457e+06 1e-11 0.120854 0.879146 2.33681e-14 8.43651e-118 0 3.72285e-11 +1700 0.5 2.2329e+17 2.46581e+06 1e-11 0.311075 0.688925 9.90838e-16 5.26533e-135 0 2.2329e-11 +1800 0.5 1.53962e+17 1.07924e+06 1e-11 0.638071 0.361929 3.94829e-17 6.01736e-153 0 1.53962e-11 +1900 0.5 1.19579e+17 391585 1e-11 0.922158 0.0778423 5.25907e-19 5.8183e-175 0 1.19579e-11 +2000 0.5 1.06072e+17 121442 1e-11 0.992096 0.0079044 3.18167e-21 2.19751e-200 0 1.06071e-11 pressure: 5 GPa -1000 0.982014 4.06599e+19 8.11198e+08 1e-11 8.52316e-08 7.64121e-06 0.00625559 0.993737 4.06599e-09 -1100 0.982014 4.0594e+19 8.0988e+08 1e-11 7.12921e-06 0.00536781 0.0784945 0.916131 4.0594e-09 -1200 0.982014 3.57677e+19 7.13354e+08 1e-11 0.000251506 0.814985 0.183155 0.0016077 3.57677e-09 -1300 0.982014 1.01801e+19 2.01602e+08 1e-11 0.00161362 0.998378 7.89782e-06 5.82951e-31 1.01801e-09 -1400 0.982014 3.34138e+18 6.48275e+07 1e-11 0.00754067 0.992459 3.15133e-09 1.3456e-55 3.34138e-10 -1500 0.982014 1.30692e+18 2.41383e+07 1e-11 0.0285575 0.971443 9.10293e-12 4.74544e-77 1.30691e-10 -1600 0.982014 6.01678e+17 1.00336e+07 1e-11 0.0903463 0.909654 1.03399e-13 4.10849e-96 6.01679e-11 -1700 0.982014 3.23232e+17 4.46465e+06 1e-11 0.240983 0.759017 2.73137e-15 1.07226e-113 3.23232e-11 -1800 0.982014 1.99404e+17 1.98807e+06 1e-11 0.527179 0.472821 9.07293e-17 2.90092e-131 1.99404e-11 -1900 0.982014 1.38804e+17 776080 1e-11 0.855044 0.144956 1.43049e-18 1.08706e-151 1.38805e-11 -2000 0.982014 1.12372e+17 247431 1e-11 0.982283 0.0177172 7.09191e-21 1.63545e-176 1.12372e-11 +1000 0.982014 4.06599e+19 8.11198e+08 1e-11 8.52316e-08 7.64121e-06 0.00625559 0.993737 0 4.06599e-09 +1100 0.982014 4.0594e+19 8.0988e+08 1e-11 7.12921e-06 0.00536781 0.0784945 0.916131 0 4.0594e-09 +1200 0.982014 3.57677e+19 7.13354e+08 1e-11 0.000251506 0.814985 0.183155 0.0016077 0 3.57677e-09 +1300 0.982014 1.01801e+19 2.01602e+08 1e-11 0.00161362 0.998378 7.89782e-06 5.82951e-31 0 1.01801e-09 +1400 0.982014 3.34138e+18 6.48275e+07 1e-11 0.00754067 0.992459 3.15133e-09 1.3456e-55 0 3.34138e-10 +1500 0.982014 1.30692e+18 2.41383e+07 1e-11 0.0285575 0.971443 9.10293e-12 4.74544e-77 0 1.30691e-10 +1600 0.982014 6.01678e+17 1.00336e+07 1e-11 0.0903463 0.909654 1.03399e-13 4.10849e-96 0 6.01679e-11 +1700 0.982014 3.23232e+17 4.46465e+06 1e-11 0.240983 0.759017 2.73137e-15 1.07226e-113 0 3.23232e-11 +1800 0.982014 1.99404e+17 1.98807e+06 1e-11 0.527179 0.472821 9.07293e-17 2.90092e-131 0 1.99404e-11 +1900 0.982014 1.38804e+17 776080 1e-11 0.855044 0.144956 1.43049e-18 1.08706e-151 0 1.38805e-11 +2000 0.982014 1.12372e+17 247431 1e-11 0.982283 0.0177172 7.09191e-21 1.63545e-176 0 1.12372e-11 OK Number of active cells: 100 (on 1 levels) Number of degrees of freedom: 5,534 (3,969+242+1,323)