From 1d4cccd3a20d14e82f578f3e4ea2cf1c9279b21c Mon Sep 17 00:00:00 2001 From: Martin Kronbichler Date: Thu, 9 May 2024 17:26:48 +0200 Subject: [PATCH] Test step-67 and step-76 tutorial programs --- tests/examples/step-67.diff | 20 ++++++++++++++++++++ tests/examples/step-67.output | 11 +++++++++++ tests/examples/step-76.diff | 20 ++++++++++++++++++++ tests/examples/step-76.output | 16 ++++++++++++++++ 4 files changed, 67 insertions(+) create mode 100644 tests/examples/step-67.diff create mode 100644 tests/examples/step-67.output create mode 100644 tests/examples/step-76.diff create mode 100644 tests/examples/step-76.output diff --git a/tests/examples/step-67.diff b/tests/examples/step-67.diff new file mode 100644 index 0000000000..41c7f399c6 --- /dev/null +++ b/tests/examples/step-67.diff @@ -0,0 +1,20 @@ +75c75 +< constexpr unsigned int n_global_refinements = 3; +--- +> constexpr unsigned int n_global_refinements = 2; +82c82 +< constexpr double final_time = testcase == 0 ? 10 : 2.0; +--- +> constexpr double final_time = testcase == 0 ? 5 : 2.0; +2196,2198d2195 +< const unsigned int n_vect_number = VectorizedArray::size(); +< const unsigned int n_vect_bits = 8 * sizeof(Number) * n_vect_number; +< +2202,2206d2198 +< pcout << "Vectorization over " << n_vect_number << ' ' +< << (std::is_same_v ? "doubles" : "floats") << " = " +< << n_vect_bits << " bits (" +< << Utilities::System::get_current_vectorization_level() << ')' +< << std::endl; +2284d2275 +< timer.print_wall_time_statistics(MPI_COMM_WORLD); diff --git a/tests/examples/step-67.output b/tests/examples/step-67.output new file mode 100644 index 0000000000..14c7d95f1d --- /dev/null +++ b/tests/examples/step-67.output @@ -0,0 +1,11 @@ +Running with 1 MPI processes +Number of degrees of freedom: 36864 ( = 4 [vars] x 256 [cells] x 36 [dofs/cell/var] ) +Time step size: 0.0137935, minimal h: 0.625, initial transport scaling: 0.205622 + +Time: 0, dt: 0.014, error rho: 1.867e-05, rho * u: 7.97e-05, energy: 0.0001589 +Time: 1.01, dt: 0.014, error rho: 5.715e-05, rho * u: 0.0001503, energy: 0.0002395 +Time: 2, dt: 0.014, error rho: 7.119e-05, rho * u: 0.0001799, energy: 0.0002681 +Time: 3.01, dt: 0.014, error rho: 7.981e-05, rho * u: 0.0001989, energy: 0.0002824 +Time: 4, dt: 0.014, error rho: 9.041e-05, rho * u: 0.0002087, energy: 0.0002955 +Time: 5.01, dt: 0.014, error rho: 7.695e-05, rho * u: 0.000158, energy: 0.0002245 + diff --git a/tests/examples/step-76.diff b/tests/examples/step-76.diff new file mode 100644 index 0000000000..4c7d92d47b --- /dev/null +++ b/tests/examples/step-76.diff @@ -0,0 +1,20 @@ +67c67 +< constexpr unsigned int fe_degree = 5; +--- +> constexpr unsigned int fe_degree = 3; +89c89 +< constexpr double final_time = testcase == 0 ? 10 : 2.0; +--- +> constexpr double final_time = testcase == 0 ? 10 : 0.5; +1519,1521d1518 +< const unsigned int n_vect_number = VectorizedArrayType::size(); +< const unsigned int n_vect_bits = 8 * sizeof(Number) * n_vect_number; +< +1525,1529d1521 +< pcout << "Vectorization over " << n_vect_number << ' ' +< << (std::is_same_v ? "doubles" : "floats") << " = " +< << n_vect_bits << " bits (" +< << Utilities::System::get_current_vectorization_level() << ')' +< << std::endl; +1593d1584 +< timer.print_wall_time_statistics(MPI_COMM_WORLD); diff --git a/tests/examples/step-76.output b/tests/examples/step-76.output new file mode 100644 index 0000000000..dac4117ab4 --- /dev/null +++ b/tests/examples/step-76.output @@ -0,0 +1,16 @@ +Running with 1 MPI processes +Number of degrees of freedom: 102400 ( = 4 [vars] x 1600 [cells] x 16 [dofs/cell/var] ) +Time step size: 0.000636798, minimal h: 0.0075, initial transport scaling: 0.00441186 + +Time: 0, dt: 0.00064, norm rho: 3.914e-16, rho * u: 1.431e-16, energy: 8.334e-16 +Time: 0.0503, dt: 0.00054, norm rho: 0.0207, rho * u: 0.03801, energy: 0.08745 +Time: 0.1, dt: 0.00051, norm rho: 0.02202, rho * u: 0.04509, energy: 0.08908 +Time: 0.15, dt: 0.00049, norm rho: 0.02256, rho * u: 0.04585, energy: 0.08944 +Time: 0.2, dt: 0.0005, norm rho: 0.02052, rho * u: 0.04354, energy: 0.08191 +Time: 0.25, dt: 0.00051, norm rho: 0.01688, rho * u: 0.04195, energy: 0.06839 +Time: 0.3, dt: 0.00051, norm rho: 0.01649, rho * u: 0.04002, energy: 0.0658 +Time: 0.35, dt: 0.00049, norm rho: 0.01768, rho * u: 0.04255, energy: 0.07025 +Time: 0.4, dt: 0.0005, norm rho: 0.01599, rho * u: 0.04491, energy: 0.06383 +Time: 0.45, dt: 0.0005, norm rho: 0.01523, rho * u: 0.04347, energy: 0.06164 +Time: 0.5, dt: 0.0005, norm rho: 0.01567, rho * u: 0.04187, energy: 0.06332 + -- 2.39.5