From 1fbfbc5134806c80cc2fe2b61623fe7530fc34b8 Mon Sep 17 00:00:00 2001 From: Piyush Sharda Date: Mon, 26 Aug 2024 23:56:49 +0200 Subject: [PATCH] add checks at intermediate densities for high-z test --- .github/workflows/burn_cell_metal_chem.yml | 2 +- unit_test/burn_cell_metal_chem/burn_cell.H | 5 +++++ .../reference_solution_1_z10.out | 18 ++++++++++++++++++ 3 files changed, 24 insertions(+), 1 deletion(-) diff --git a/.github/workflows/burn_cell_metal_chem.yml b/.github/workflows/burn_cell_metal_chem.yml index b188353d5..419b08aa5 100644 --- a/.github/workflows/burn_cell_metal_chem.yml +++ b/.github/workflows/burn_cell_metal_chem.yml @@ -80,7 +80,7 @@ jobs: run: | set -e cd unit_test/burn_cell_metal_chem - line_numbers=(4461 4463 4465 4467 4468 4471 4472 4475 4481 4486) + line_numbers=(1612 1614 4479 4481 4483 4485 4486 4489 4490 4493 4499 4504) threshold=0.01 error_found=false for line_number in "${line_numbers[@]}"; do diff --git a/unit_test/burn_cell_metal_chem/burn_cell.H b/unit_test/burn_cell_metal_chem/burn_cell.H index 31f0e4a49..296467fc3 100644 --- a/unit_test/burn_cell_metal_chem/burn_cell.H +++ b/unit_test/burn_cell_metal_chem/burn_cell.H @@ -344,6 +344,11 @@ auto burn_cell_c() -> int { state_over_time << std::endl; std::cout << "step " << n << " done with dt = " << dt << std::endl; + // print some output at intermediate densities for testing purposes if redshift > 0 + if (state.rho > 1e-18 && state.rho < 1.2e-18 && redshift > 0) { + std::cout << "Tgas = " << state.T << std::endl; + } + } state_over_time.close(); diff --git a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out index 722332478..d67ac10d0 100644 --- a/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out +++ b/unit_test/burn_cell_metal_chem/reference_solution_1_z10.out @@ -1585,23 +1585,41 @@ step 1537 done with dt = 2.137708209e+10 step 1538 done with dt = 2.127099153e+10 step 1539 done with dt = 2.116542749e+10 step 1540 done with dt = 2.106038734e+10 +Tgas = 31.53250536 step 1541 done with dt = 2.09558685e+10 +Tgas = 31.52704064 step 1542 done with dt = 2.085186837e+10 +Tgas = 31.52147312 step 1543 done with dt = 2.074838438e+10 +Tgas = 31.51580941 step 1544 done with dt = 2.064541397e+10 +Tgas = 31.51005587 step 1545 done with dt = 2.054295459e+10 +Tgas = 31.50421855 step 1546 done with dt = 2.04410037e+10 +Tgas = 31.49830335 step 1547 done with dt = 2.033955877e+10 +Tgas = 31.49231585 step 1548 done with dt = 2.023861731e+10 +Tgas = 31.48626145 step 1549 done with dt = 2.01381768e+10 +Tgas = 31.48014529 step 1550 done with dt = 2.003823477e+10 +Tgas = 31.47397231 step 1551 done with dt = 1.993878873e+10 +Tgas = 31.46774723 step 1552 done with dt = 1.983983623e+10 +Tgas = 31.46147455 step 1553 done with dt = 1.974137481e+10 +Tgas = 31.4551586 step 1554 done with dt = 1.964340205e+10 +Tgas = 31.44880348 step 1555 done with dt = 1.954591551e+10 +Tgas = 31.44241076 step 1556 done with dt = 1.944891278e+10 +Tgas = 31.43598915 step 1557 done with dt = 1.935239145e+10 +Tgas = 31.42953383 step 1558 done with dt = 1.925634915e+10 step 1559 done with dt = 1.91607835e+10 step 1560 done with dt = 1.906569212e+10