1 files changed, 140 insertions, 0 deletions
diff --git a/src/devices/value_list_source.h b/src/devices/value_list_source.h
new file mode 100644
index 0000000..a7ddb6e
--- /dev/null
+++ b/src/devices/value_list_source.h
@@ -0,0 +1,140 @@
+#pragma once
+
+#include <systemc.h>
+#include <fstream>
+#include <algorithm>
+#include <vector>
+#include <utility>
+
+#include "strutils.h"
+#include "optical_output_port.h"
+#include "optical_signal.h"
+#include "specs.h"
+#include "spx_module.h"
+
+using std::vector;
+using std::pair;
+
+class VLSource : public spx_module {
+public:
+    typedef pair<double, spx::oa_value_type> time_value_pair_type;
+    // Ports
+    spx::oa_port_out_type p_out;
+
+    // Timed ports writers
+    OpticalOutputPort m_out_writer;
+
+    // Signal to emit
+    vector<time_value_pair_type> m_values_queue;
+
+    // Source emission control
+    spx::ed_signal_type enable;
+
+    // Processes
+    void runner();
+
+    VLSource(sc_module_name name, const vector<time_value_pair_type> &values = {})
+    : spx_module(name)
+    , m_out_writer((string(this->name()) + "_out_writer").c_str(), p_out)
+    , m_values_queue(values)
+    {
+        sortValues();
+        enable = sc_logic(0);
+
+        SC_HAS_PROCESS(VLSource);
+        SC_THREAD(runner);
+    }
+
+    void setValues(const vector<time_value_pair_type> &values)
+    {
+        m_values_queue = values;
+        sortValues();
+    }
+
+    void setValues(const string &filename)
+    {
+        std::ifstream file(filename);
+
+        if (!file.is_open()) {
+            cerr << "Error: Cannot open the file." << endl;
+            exit(1);
+        }
+
+        string line;
+
+        getline(file, line);
+
+        int lineNumber = 0;
+        while (getline(file, line)) {
+            ++lineNumber;
+            istringstream s(line);
+            string field;
+
+            double t, P, WL;
+
+            // Skip if line is empty or starts with ;
+            if (line.empty())
+            {
+                cout << "Skipping empty line." << endl;
+                continue;
+            }
+            // Skip if line starts with ;
+            if (line[0] == ';')
+            {
+                cout << "Skipping commented line: " << line << endl;
+                continue;
+            }
+
+            // Read first field (time)
+            getline(s, field, ',');
+            strutils::trim(field);
+
+            // Skip header line if it exists
+            if(lineNumber == 1 && field == "time")
+            {
+                cout << "Skipping header line: " << line << endl;
+                continue;
+            }
+
+            try {
+                t = stod(field);
+            } catch (std::invalid_argument const& ex) {
+                cerr << "Invalid value for time: \"" << field << "\"" << endl;
+                exit(1);
+            }
+
+            // Read second field (P)
+            getline(s, field, ',');
+            try {
+                P = stod(field);
+            } catch (std::invalid_argument const& ex) {
+                cerr << "Invalid value for power: \"" << field << "\"" << endl;
+                exit(1);
+            }
+
+            // Read wavelength (WL)
+            std::getline(s, field);
+            try {
+                WL = stod(field);
+            } catch (std::invalid_argument const& ex) {
+                cerr << "Invalid value for wavelength: \"" << field << "\"" << endl;
+                exit(1);
+            }
+
+            cout << t << ": " << P << "W @" << WL*1e9 << "nm" << endl;
+            m_values_queue.emplace_back(t, OpticalSignal(sqrt(P), WL));
+        }
+
+        file.close();
+        sortValues();
+    }
+
+    void sortValues()
+    {
+        auto cmp = [](const time_value_pair_type& p1, const time_value_pair_type& p2) {
+            return p1.first < p2.first;
+        };
+        // Stable sort will keep initial ordering for values with identical time
+        std::stable_sort(m_values_queue.begin(), m_values_queue.end(), cmp);
+    }
+};