From e8781cfdd74f8e0a01f0c03ea0b47d37dc866453 Mon Sep 17 00:00:00 2001 From: EstherLerouzic Date: Thu, 27 Apr 2023 16:57:00 +0200 Subject: [PATCH] Feat: apply roadm-path loss Signed-off-by: EstherLerouzic Change-Id: I87e27d5b653bdce814f43e4c9c1183fb51fbcc1e --- gnpy/core/elements.py | 15 ++++++++------- tests/test_roadm_restrictions.py | 20 ++++++++++++-------- 2 files changed, 20 insertions(+), 15 deletions(-) diff --git a/gnpy/core/elements.py b/gnpy/core/elements.py index a78a5b733..8255984ed 100644 --- a/gnpy/core/elements.py +++ b/gnpy/core/elements.py @@ -374,13 +374,14 @@ def get_per_degree_power(self, degree, spectral_info): def propagate(self, spectral_info, degree, from_degree): """Equalization targets are read from topology file if defined and completed with default definition of the library. - If the input power is lower than the target one, use the input power instead because - a ROADM doesn't amplify, it can only attenuate. - There is no difference for add or express : the same target is applied. For the moment - propagates operates with spectral info carriers all having the same source or destination. + If the input power is lower than the target one, use the input power minus the ROADM loss + if is exists, because a ROADM doesn't amplify, it can only attenuate. + There is no difference for add or express : the same target is applied. + For the moment propagate operates with spectral info carriers all having the same source or destination. """ - # TODO maybe add a minimum loss for the ROADM - + # apply min ROADM loss if it exists + roadm_maxloss_db = self.get_roadm_path(from_degree, degree).impairment.maxloss + spectral_info.apply_attenuation_db(roadm_maxloss_db) # find the target power for the reference carrier ref_per_degree_pch = self.get_per_degree_ref_power(degree) # find the target powers for each signal carrier @@ -391,7 +392,7 @@ def propagate(self, spectral_info, degree, from_degree): # the target power out configured for this ROADM degree's egress. Since ROADM does not amplify, # the power out of the ROADM for the ref channel is the min value between target power and input power. # (TODO add a minimum loss for the ROADM crossing) - self.ref_pch_out_dbm = min(self.ref_pch_in_dbm[from_degree], ref_per_degree_pch) + self.ref_pch_out_dbm = min(self.ref_pch_in_dbm[from_degree] - roadm_maxloss_db, ref_per_degree_pch) # Definition of effective_loss: # Optical power of carriers are equalized by the ROADM, so that the experienced loss is not the same for # different carriers. effective_loss records the loss for the reference carrier. diff --git a/tests/test_roadm_restrictions.py b/tests/test_roadm_restrictions.py index 18dab8682..ebe76fa56 100644 --- a/tests/test_roadm_restrictions.py +++ b/tests/test_roadm_restrictions.py @@ -213,9 +213,11 @@ def test_restrictions(restrictions, equipment): raise AssertionError() +@pytest.mark.parametrize('roadm_type_variety, roadm_b_maxloss', [('default', 0), + ('example_detailed_impairments', 16.5)]) @pytest.mark.parametrize('power_dbm', [0, +1, -2]) @pytest.mark.parametrize('prev_node_type, effective_pch_out_db', [('edfa', -20.0), ('fused', -22.0)]) -def test_roadm_target_power(prev_node_type, effective_pch_out_db, power_dbm): +def test_roadm_target_power(prev_node_type, effective_pch_out_db, power_dbm, roadm_type_variety, roadm_b_maxloss): """Check that egress power of roadm is equal to target power if input power is greater than target power else, that it is equal to input power. Use a simple two hops A-B-C topology for the test where the prev_node in ROADM B is either an amplifier or a fused, so that the target @@ -226,6 +228,8 @@ def test_roadm_target_power(prev_node_type, effective_pch_out_db, power_dbm): json_network = load_json(TEST_DIR / 'data/twohops_roadm_power_test.json') prev_node = next(n for n in json_network['elements'] if n['uid'] == 'west edfa in node B to ila2') json_network['elements'].remove(prev_node) + roadm_b = next(element for element in json_network['elements'] if element['uid'] == 'roadm node B') + roadm_b['type_variety'] = roadm_type_variety if prev_node_type == 'edfa': prev_node = {'uid': 'west edfa in node B to ila2', 'type': 'Edfa'} elif prev_node_type == 'fused': @@ -269,10 +273,9 @@ def test_roadm_target_power(prev_node_type, effective_pch_out_db, power_dbm): # if previous was an EDFA, power level at ROADM input is enough for the ROADM to apply its # target power (as specified in equipment ie -20 dBm) # if it is a Fused, the input power to the ROADM is smaller than the target power, and the - # ROADM cannot apply this target. In this case, it is assumed that the ROADM has 0 dB loss - # so the output power will be the same as the input power, which for this particular case - # corresponds to -22dBm + power_dbm - # next step (for ROADM modelling) will be to apply a minimum loss for ROADMs ! + # ROADM cannot apply this target. If the ROADM has 0 dB loss the output power will be the same + # as the input power, which for this particular case corresponds to -22dBm + power_dbm. + # If ROADM has a minimum losss, then output power will be -22dBm + power_dbm - ROADM loss. ! if prev_node_type == 'edfa': # edfa prev_node sets input power to roadm to a high enough value: # check that target power is correctly set in the ROADM @@ -283,9 +286,10 @@ def test_roadm_target_power(prev_node_type, effective_pch_out_db, power_dbm): # fused prev_node does not reamplify power after fiber propagation, so input power # to roadm is low. # check that target power correctly reports power_dbm from previous propagation - assert_allclose(el.ref_pch_out_dbm, effective_pch_out_db + power_dbm, rtol=1e-3) - # Check that egress power of roadm is not equalized: power out is the same as power in. - assert_allclose(power_out_roadm, power_in_roadm, rtol=1e-3) + assert_allclose(el.ref_pch_out_dbm, effective_pch_out_db + power_dbm - roadm_b_maxloss, rtol=1e-3) + # Check that egress power of roadm is not equalized: + # power out is the same as power in minus the ROADM loss. + assert_allclose(power_out_roadm, power_in_roadm / db2lin(roadm_b_maxloss), rtol=1e-3) assert effective_pch_out_db + power_dbm ==\ pytest.approx(lin2db(min(power_in_roadm) * 1e3), rel=1e-3) else: