From 18f9e131072c15daa9b45334112babbba911b191 Mon Sep 17 00:00:00 2001
From: Emily Boudreaux <thomas@boudreauxmail.com>
Date: Mon, 14 Jul 2025 14:54:22 -0400
Subject: [PATCH] test(sandbox-test): updated sandbox

---
 tests/graphnet_sandbox/main.cpp | 64 ++++++++++-----------------------
 1 file changed, 18 insertions(+), 46 deletions(-)

diff --git a/tests/graphnet_sandbox/main.cpp b/tests/graphnet_sandbox/main.cpp
index 175a947f..08bf973f 100644
--- a/tests/graphnet_sandbox/main.cpp
+++ b/tests/graphnet_sandbox/main.cpp
@@ -27,12 +27,10 @@
 
 #include "gridfire/partition/composite/partition_composite.h"
 
-// Keep a copy of the previous handler
 static std::terminate_handler g_previousHandler = nullptr;
 
 void measure_execution_time(const std::function<void()>& callback, const std::string& name)
 {
-    // variable names in camelCase
     const auto startTime   = std::chrono::steady_clock::now();
     callback();
     const auto endTime     = std::chrono::steady_clock::now();
@@ -43,14 +41,7 @@ void measure_execution_time(const std::function<void()>& callback, const std::st
 
 void quill_terminate_handler()
 {
-    // 1. Stop the Quill backend (flushes all sinks and joins thread)
     quill::Backend::stop();
-
-    // 2. Optionally, flush individual loggers:
-    // for (auto logger : quill::get_all_loggers())
-    //   logger->flush_log();
-
-    // 3. Chain to the previous terminate handler (or abort if none)
     if (g_previousHandler)
         g_previousHandler();
     else
@@ -64,8 +55,8 @@ int main() {
     LOG_DEBUG(logger, "Starting Adaptive Engine View Example...");
 
     using namespace gridfire;
-    const std::vector<double> comp = {0.708, 0.0, 2.94e-5, 0.276, 0.003, 0.0011, 9.62e-3, 1.62e-3, 5.16e-4};
-    const std::vector<std::string> symbols = {"H-1", "H-2", "He-3", "He-4", "C-12", "N-14", "O-16", "Ne-20", "Mg-24"};
+    const std::vector<double> comp = {0.708, 2.94e-5, 0.276, 0.003, 0.0011, 9.62e-3, 1.62e-3, 5.16e-4};
+    const std::vector<std::string> symbols = {"H-1", "He-3", "He-4", "C-12", "N-14", "O-16", "Ne-20", "Mg-24"};
 
 
     fourdst::composition::Composition composition;
@@ -77,44 +68,25 @@ int main() {
         BasePartitionType::RauscherThielemann,
         BasePartitionType::GroundState
     });
-    GraphEngine ReaclibEngine(composition, partitionFunction);
-
-
 
     NetIn netIn;
     netIn.composition = composition;
-    netIn.temperature = 1.5e7;
-    netIn.density = 1e2;
-    netIn.energy = 0.0;
+    netIn.temperature = 1.5e6;
+    netIn.density = 1.5e3;
+    netIn.energy = 0;
+    netIn.tMax = 3e17;
+    netIn.dt0 = 1e-12;
 
-    netIn.tMax = 1; // 1 year in seconds
-    netIn.dt0 = 1e12;
+    GraphEngine ReaclibEngine(composition, partitionFunction, NetworkBuildDepth::SecondOrder);
+    auto primedReport = ReaclibEngine.primeEngine(netIn);
+    std::cout << primedReport << std::endl;
+    std::cout << "Initial Composition\n";
+    for (const auto& [symbol, entry] : netIn.composition) {
+        std::cout << "\t" << symbol << ": " << entry.mass_fraction() << "\n";
+    }
 
-    NetOut netOut;
-    netIn.dt0 = 1e-15;
-
-    const auto primedComposition = primeNetwork(netIn, ReaclibEngine);
-    netIn.composition = primedComposition;
-    netIn.composition.finalize(true);
-    MultiscalePartitioningEngineView multiscaleEngine(ReaclibEngine);
-    multiscaleEngine.equilibrateNetwork(netIn, 1e12);
-
-
-    // // std::cout << ReaclibEngine.getPartitionFunction().type() << std::endl;
-    // ReaclibEngine.setPrecomputation(true);
-    // ReaclibEngine.setUseReverseReactions(false);
-    // // // AdaptiveEngineView adaptiveEngine(ReaclibEngine);
-    // io::SimpleReactionListFileParser parser{};
-    // FileDefinedEngineView approx8EngineView(ReaclibEngine, "approx8.net", parser);
-    // approx8EngineView.setScreeningModel(screening::ScreeningType::WEAK);
-    // solver::DirectNetworkSolver solver(approx8EngineView);
-    // netOut = solver.evaluate(netIn);
-    // std::cout << netOut.composition << std::endl;
-
-
-    // measure_execution_time([&]() {
-    //     netOut = solver.evaluate(netIn);
-    // }, "Approx8 Network Evaluation (Precomputation)");
-    // ReaclibEngine.setPrecomputation(false);
-    // std::cout << "Precomputation H-1: " << netOut.composition.getMassFraction("H-1") << " in " << netOut.num_steps << " steps." << std::endl;
+    std::cout << "\nPrimed Composition\n";
+    for (const auto& [symbol, entry] : primedReport.primedComposition) {
+        std::cout << "\t" << symbol << ": " << entry.mass_fraction() << "\n";
+    }
 }
\ No newline at end of file