Merge pull request #482 from lucafedeli88/qed_evolve_optical_depth

QED: add evolution of the optical depth according to Breit-Wheeler and Quantum Synchrotron processes

1 files changed, 29 insertions, 0 deletions

diff --git a/Source/QED/QuantumSyncEngineWrapper.cpp b/Source/QED/QuantumSyncEngineWrapper.cpp
index b2630dc4d..ffafec761 100644
--- a/Source/QED/QuantumSyncEngineWrapper.cpp
+++ b/Source/QED/QuantumSyncEngineWrapper.cpp
@@ -1,6 +1,7 @@
 #include "QuantumSyncEngineWrapper.H"
 
 #include "QedTableParserHelperFunctions.H"
+#include "QuantumSyncDummyTable.H"
 
 #include <utility>
 
@@ -139,6 +140,28 @@ QuantumSynchrotronEngine::init_lookup_tables_from_raw_data (
     return true;
 }
 
+void QuantumSynchrotronEngine::init_dummy_tables()
+{
+    m_innards.ctrl = QedUtils::QuantumSyncEngineInnardsDummy.ctrl;
+    m_innards.KKfunc_coords.assign(
+        QedUtils::QuantumSyncEngineInnardsDummy.KKfunc_coords.begin(),
+        QedUtils::QuantumSyncEngineInnardsDummy.KKfunc_coords.end());
+    m_innards.KKfunc_data.assign(
+        QedUtils::QuantumSyncEngineInnardsDummy.KKfunc_data.begin(),
+        QedUtils::QuantumSyncEngineInnardsDummy.KKfunc_data.end());
+    m_innards.cum_distrib_coords_1.assign(
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_coords_1.begin(),
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_coords_1.end());
+    m_innards.cum_distrib_coords_2.assign(
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_coords_2.begin(),
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_coords_2.end());
+    m_innards.cum_distrib_data.assign(
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_data.begin(),
+        QedUtils::QuantumSyncEngineInnardsDummy.cum_distrib_data.end());
+
+    m_lookup_tables_initialized = true;
+}
+
 Vector<char> QuantumSynchrotronEngine::export_lookup_tables_data () const
 {
     Vector<char> res{};
@@ -175,6 +198,12 @@ QuantumSynchrotronEngine::get_default_ctrl() const
     return PicsarQuantumSynchrotronCtrl();
 }
 
+const PicsarQuantumSynchrotronCtrl&
+QuantumSynchrotronEngine::get_ref_ctrl() const
+{
+    return m_innards.ctrl;
+}
+
 void QuantumSynchrotronEngine::compute_lookup_tables (
     PicsarQuantumSynchrotronCtrl ctrl)
 {