blob: dc018da40882d035dc520bff114106a5525461ca (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
|
#include <algorithm>
#include <iostream>
#include <pcm_device.h>
#define ZERO_PULSE_THRESHOLD (1e-10)
using namespace std;
void PCMElement::on_input_changed()
{
// setting the transmission for the initial state
update_transmission_local();
if (specsGlobalConfig.verbose_component_initialization)
{
cout << name() << ":" << endl;
cout << "Emelt = " << m_meltEnergy*1e9 << " nJ" << endl;
cout << "Nstates = " << m_nStates << "" << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p_in.get_interface()))->name();
cout << " --> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p_out.get_interface()))->name();
cout << endl;
cout << endl;
}
m_last_pulse_power = 0;
bool in_window = false;
vector<pulse_sample_t> vector_of_samples;
uint32_t cur_wavelength_id;
m_memory_in[0] = 0;
double total_in_power = 0;
while (true) {
// Wait next input change
wait();
// Read signal and store field in memory for that wavelength
auto s = p_in->read();
cur_wavelength_id = s.m_wavelength_id;
m_memory_in[cur_wavelength_id] = s.m_field;
// Summing powers of all wavelengths
total_in_power = 0;
for (auto lambdaID_field : m_memory_in)
{
total_in_power += norm(lambdaID_field.second);
}
// Storage of optical state on the input
const sc_time &now = sc_time_stamp();
bool last_was_zero = m_last_pulse_power < ZERO_PULSE_THRESHOLD;
bool current_is_zero = total_in_power < ZERO_PULSE_THRESHOLD;
bool rising_edge = last_was_zero && !current_is_zero;
if (!rising_edge)
{
// Enter here when:
// - a signal is received on top of another one (pulse intersect)
// - a signal ends (power falls to 0)
// Record the event
vector_of_samples.emplace_back(now.to_seconds(), total_in_power);
}
else
{
// Enter here when:
// - a signal starts (power rises from 0 to non-0 value)
// This case is necessarily the first rise
// Will unroll events and check if there was a phase change before according to
// the vector and advance the state accordingly to get the transmission
// cout << "\t\t first rise detected" << endl;
in_window = phase_change(vector_of_samples, true);
if (!in_window)
vector_of_samples.clear();
// Record the current event
vector_of_samples.emplace_back(now.to_seconds(), total_in_power);
}
m_last_pulse_power = total_in_power;
// Write attenuated signal to output
s *= m_transmission;
s.getNewId();
m_out_writer.delayedWrite(s, SC_ZERO_TIME);
}
}
bool PCMElement::phase_change(const vector<pulse_sample_t> &samples, const bool &local)
{
// Here should decide between sending the
// sim off to external or resolving inside
// Should set new values for m_state and m_transmission
double energy_absorbed = 0;
double dur_pulse = 0;
// will return true if we shouldn't erase the vector
bool in_influence_window = false;
if (samples.empty())
return true;
if (sc_time_stamp().to_seconds() - samples.back().first < influence_time_ns*1e-9)
return true;
if (local)
{
// Integrating power along each duration
// pulse 0 stays at pow[0] for tim[1]-tim[0] seconds
// the last sample is zero in power, just here to get the final time
for (unsigned int i = 0; i < samples.size()-1; i++)
{
dur_pulse = samples[i+1].first - samples[i].first;
energy_absorbed += samples[i].second * dur_pulse;
}
// if enough energy and not saturated
if ((energy_absorbed > m_meltEnergy) && (m_state < m_nStates))
{
// cout << sc_time_stamp() << ": \n\t" << name() << " phase changed with "
// << energy_absorbed*1e6 << " uJ" << endl;
m_state++;
// cout << "\tnew state: " << m_state << endl;
update_transmission_local();
}
}
else // external phase change simulation
{
// send the simulation to external software
cout << "TODO: implement handles to external sim"
<< endl;
}
return in_influence_window;
}
void PCMElement::update_transmission_local()
{
double sp = 0.85;
double ep = 0.95;
double speed = 3;
m_transmission = sp+(ep-sp)*tanh(speed*m_state/m_nStates);
// cout << "\tnew trans_power: " << m_transmission << "- W/W" << endl;
// For field need to do a square root
m_transmission = sqrt(m_transmission);
}
|
