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
|
#include "specs.h"
#include "crossing.h"
using namespace std;
void CrossingUni::on_input_changed()
{
// If it's NAN, it's because it was not specified and thus the linear should be zero (-inf dB)
const double crosstalk_field_lin = (isnan(m_crosstalk_power_dB)) ? 0 : pow(10.0, m_crosstalk_power_dB / 20);
const double transmission_field_lin = pow(10.0, -m_attenuation_power_dB / 20); // due to attenuations
// Pre-calculate S-parameters
OpticalSignal::field_type S13, S14, S23, S24;
// the second part relates to power that is crossed over (not transmitted)
S13 = polar(transmission_field_lin, 0.0); // From in1 to out1
S24 = polar(transmission_field_lin, 0.0); // From in2 to out2
S14 = polar(crosstalk_field_lin, 0.0); // From in1 to out2 (crosstalk)
S23 = polar(crosstalk_field_lin, 0.0); // From in2 to out1 (crosstalk)
if (specsGlobalConfig.verbose_component_initialization)
{
cout << name() << ":" << endl;
cout << "transmission_power = " << norm(S13) << " W/W" << endl;
cout << "crosstalk_power = " << norm(S14)<< " W/W" << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p_in1.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p_out1.get_interface()))->name();
cout << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p_in2.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p_out2.get_interface()))->name();
cout << endl;
cout << endl;
}
while (true) {
// Wait for a new input signal
wait();
// Read current inputs
auto s1 = p_in1->read();
auto s2 = p_in2->read();
// Apply S-parameters
auto s3 = s1 * S13 + s2 * S23;
auto s4 = s1 * S14 + s2 * S24;
// Get new IDs for signal
s3.getNewId();
s4.getNewId();
// Write to ouput port after delay
m_out1_writer.delayedWrite(s3, sc_time(0, SC_NS));
m_out2_writer.delayedWrite(s4, sc_time(0, SC_NS));
}
}
void CrossingBi::on_input_changed()
{
// If it's NAN, it's because it was not specified and thus the linear should be zero (-inf dB)
const double crosstalk_field_lin = (isnan(m_crosstalk_power_dB)) ? 0 : pow(10.0, m_crosstalk_power_dB / 20);
const double transmission_field_lin = pow(10.0, -m_attenuation_power_dB / 20); // due to attenuations
// Pre-calculate S-parameters
OpticalSignal::field_type through = polar(transmission_field_lin, 0.0);
OpticalSignal::field_type cross = polar(crosstalk_field_lin, 0.0);
OpticalSignal::field_type S02, S03, S12, S13;
// the second part relates to power that is crossed over (not transmitted)
S02 = through; // From in1 to out1
S13 = through; // From in2 to out2
S03 = cross; // From in1 to out2 (crosstalk)
S12 = cross; // From in2 to out1 (crosstalk)
if (specsGlobalConfig.verbose_component_initialization)
{
cout << name() << ":" << endl;
cout << "transmission_power = " << norm(through) << " W/W" << endl;
cout << "crosstalk_power = " << norm(cross)<< " W/W" << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p0_in.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p2_out.get_interface()))->name();
cout << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p1_in.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p3_out.get_interface()))->name();
cout << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p2_in.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p0_out.get_interface()))->name();
cout << endl;
cout << (dynamic_cast<spx::oa_signal_type *>(p3_in.get_interface()))->name();
cout << " --- ---> ";
cout << (dynamic_cast<spx::oa_signal_type *>(p1_out.get_interface()))->name();
cout << endl;
cout << endl;
}
while (true) {
// Wait for a new input signal
wait();
// Read current inputs
auto s0_in = p0_in->read();
auto s1_in = p1_in->read();
auto s2_in = p2_in->read();
auto s3_in = p3_in->read();
// Apply S-parameters
auto s0_out = s2_in * S02 + s3_in * S03;
auto s1_out = s2_in * S12 + s3_in * S13;
auto s2_out = s0_in * S02 + s1_in * S12;
auto s3_out = s0_in * S03 + s1_in * S13;
// Get new IDs for signal
s0_out.getNewId();
s1_out.getNewId();
s2_out.getNewId();
s3_out.getNewId();
// Write to ouput port after delay
m_p0_out_writer.delayedWrite(s0_out, sc_time(0, SC_NS));
m_p1_out_writer.delayedWrite(s1_out, sc_time(0, SC_NS));
m_p2_out_writer.delayedWrite(s2_out, sc_time(0, SC_NS));
m_p3_out_writer.delayedWrite(s3_out, sc_time(0, SC_NS));
}
}
|