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 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282
|
from __future__ import annotations
import logging
import math
import pytest
from pint import OffsetUnitCalculusError, Unit, UnitRegistry
from pint.facets.plain.unit import UnitsContainer
from pint.testsuite import QuantityTestCase, helpers
@pytest.fixture(scope="module")
def module_registry_auto_offset():
return UnitRegistry(autoconvert_offset_to_baseunit=True)
# TODO: do not subclass from QuantityTestCase
class TestLogarithmicQuantity(QuantityTestCase):
def test_log_quantity_creation(self, caplog):
# Following Quantity Creation Pattern
for args in (
(4.2, "dBm"),
(4.2, UnitsContainer(decibelmilliwatt=1)),
(4.2, self.ureg.dBm),
):
x = self.Q_(*args)
assert x.magnitude == 4.2
assert x.units == UnitsContainer(decibelmilliwatt=1)
x = self.Q_(self.Q_(4.2, "dBm"))
assert x.magnitude == 4.2
assert x.units == UnitsContainer(decibelmilliwatt=1)
x = self.Q_(4.2, UnitsContainer(decibelmilliwatt=1))
y = self.Q_(x)
assert x.magnitude == y.magnitude
assert x.units == y.units
assert x is not y
# Using multiplications for dB units requires autoconversion to baseunits
new_reg = UnitRegistry(autoconvert_offset_to_baseunit=True)
x = new_reg.Quantity("4.2 * dBm")
assert x.magnitude == 4.2
assert x.units == UnitsContainer(decibelmilliwatt=1)
with caplog.at_level(logging.DEBUG):
assert "wally" not in caplog.text
assert 4.2 * new_reg.dBm == new_reg.Quantity(4.2, 2 * new_reg.dBm)
assert len(caplog.records) == 1
def test_log_convert(self):
# # 1 dB = 1/10 * bel
# helpers.assert_quantity_almost_equal(self.Q_(1.0, "dB").to("dimensionless"), self.Q_(1, "bell") / 10)
# # Uncomment Bell unit in default_en.txt
# ## Test dB to dB units octave - decade
# 1 decade = log2(10) octave
helpers.assert_quantity_almost_equal(
self.Q_(1.0, "decade"), self.Q_(math.log2(10), "octave")
)
# ## Test dB to dB units dBm - dBu
# 0 dBm = 1mW = 1e3 uW = 30 dBu
helpers.assert_quantity_almost_equal(
self.Q_(0.0, "dBm"), self.Q_(29.999999999999996, "dBu"), atol=1e-7
)
# ## Test dB to dB units dBm - dBW
# 0 dBW = 1W = 1e3 mW = 30 dBm
helpers.assert_quantity_almost_equal(
self.Q_(0.0, "dBW"), self.Q_(29.999999999999996, "dBm"), atol=1e-7
)
def test_mix_regular_log_units(self):
# Test regular-logarithmic mixed definition, such as dB/km or dB/cm
# Multiplications and divisions with a mix of Logarithmic Units and regular Units is normally not possible.
# The reason is that dB are considered by pint like offset units.
# Multiplications and divisions that involve offset units are badly defined, so pint raises an error
with pytest.raises(OffsetUnitCalculusError):
(-10.0 * self.ureg.dB) / (1 * self.module_registry.cm)
# However, if the flag autoconvert_offset_to_baseunit=True is given to UnitRegistry, then pint converts the unit to plain.
# With this flag on multiplications and divisions are now possible:
new_reg = UnitRegistry(autoconvert_offset_to_baseunit=True)
helpers.assert_quantity_almost_equal(
-10 * new_reg.dB / new_reg.cm, 0.1 / new_reg.cm
)
log_unit_names = [
"decibelwatt",
"dBW",
"decibelmilliwatt",
"dBm",
"decibelmicrowatt",
"dBu",
"decibel",
"dB",
"decade",
"octave",
"oct",
]
@pytest.mark.parametrize("unit_name", log_unit_names)
def test_unit_by_attribute(module_registry, unit_name):
"""Can the logarithmic units be accessed by attribute lookups?"""
unit = getattr(module_registry, unit_name)
assert isinstance(unit, Unit)
@pytest.mark.parametrize("unit_name", log_unit_names)
def test_unit_parsing(module_registry, unit_name):
"""Can the logarithmic units be understood by the parser?"""
unit = module_registry.parse_units(unit_name)
assert isinstance(unit, Unit)
@pytest.mark.parametrize("mag", [1.0, 4.2])
@pytest.mark.parametrize("unit_name", log_unit_names)
def test_quantity_by_constructor(module_registry, unit_name, mag):
"""Can Quantity() objects be constructed using logarithmic units?"""
q = module_registry.Quantity(mag, unit_name)
assert q.magnitude == pytest.approx(mag)
assert q.units == getattr(module_registry, unit_name)
@pytest.mark.parametrize("mag", [1.0, 4.2])
@pytest.mark.parametrize("unit_name", log_unit_names)
def test_quantity_by_multiplication(module_registry_auto_offset, unit_name, mag):
"""Test that logarithmic units can be defined with multiplication
Requires setting `autoconvert_offset_to_baseunit` to True
"""
unit = getattr(module_registry_auto_offset, unit_name)
q = mag * unit
assert q.magnitude == pytest.approx(mag)
assert q.units == unit
@pytest.mark.parametrize(
"unit1,unit2",
[
("decibelwatt", "dBW"),
("decibelmilliwatt", "dBm"),
("decibelmicrowatt", "dBu"),
("decibel", "dB"),
("octave", "oct"),
],
)
def test_unit_equivalence(module_registry, unit1, unit2):
"""Are certain pairs of units equivalent?"""
assert getattr(module_registry, unit1) == getattr(module_registry, unit2)
@pytest.mark.parametrize(
"db_value,scalar",
[
(0.0, 1.0), # 0 dB == 1x
(-10.0, 0.1), # -10 dB == 0.1x
(10.0, 10.0),
(30.0, 1e3),
(60.0, 1e6),
],
)
def test_db_conversion(module_registry, db_value, scalar):
"""Test that a dB value can be converted to a scalar and back."""
Q_ = module_registry.Quantity
assert Q_(db_value, "dB").to("dimensionless").magnitude == pytest.approx(scalar)
assert Q_(scalar, "dimensionless").to("dB").magnitude == pytest.approx(db_value)
@pytest.mark.parametrize(
"octave,scalar",
[
(2.0, 4.0), # 2 octave == 4x
(1.0, 2.0), # 1 octave == 2x
(0.0, 1.0),
(-1.0, 0.5),
(-2.0, 0.25),
],
)
def test_octave_conversion(module_registry, octave, scalar):
"""Test that an octave can be converted to a scalar and back."""
Q_ = module_registry.Quantity
assert Q_(octave, "octave").to("dimensionless").magnitude == pytest.approx(scalar)
assert Q_(scalar, "dimensionless").to("octave").magnitude == pytest.approx(octave)
@pytest.mark.parametrize(
"decade,scalar",
[
(2.0, 100.0), # 2 decades == 100x
(1.0, 10.0), # 1 octave == 2x
(0.0, 1.0),
(-1.0, 0.1),
(-2.0, 0.01),
],
)
def test_decade_conversion(module_registry, decade, scalar):
"""Test that a decade can be converted to a scalar and back."""
Q_ = module_registry.Quantity
assert Q_(decade, "decade").to("dimensionless").magnitude == pytest.approx(scalar)
assert Q_(scalar, "dimensionless").to("decade").magnitude == pytest.approx(decade)
@pytest.mark.parametrize(
"dbm_value,mw_value",
[
(0.0, 1.0), # 0.0 dBm == 1.0 mW
(10.0, 10.0),
(20.0, 100.0),
(-10.0, 0.1),
(-20.0, 0.01),
],
)
def test_dbm_mw_conversion(module_registry, dbm_value, mw_value):
"""Test dBm values can convert to mW and back."""
Q_ = module_registry.Quantity
assert Q_(dbm_value, "dBm").to("mW").magnitude == pytest.approx(mw_value)
assert Q_(mw_value, "mW").to("dBm").magnitude == pytest.approx(dbm_value)
@pytest.mark.xfail
def test_compound_log_unit_multiply_definition(module_registry_auto_offset):
"""Check that compound log units can be defined using multiply."""
Q_ = module_registry_auto_offset.Quantity
canonical_def = Q_(-161, "dBm") / module_registry_auto_offset.Hz
mult_def = -161 * module_registry_auto_offset("dBm/Hz")
assert mult_def == canonical_def
@pytest.mark.xfail
def test_compound_log_unit_quantity_definition(module_registry_auto_offset):
"""Check that compound log units can be defined using ``Quantity()``."""
Q_ = module_registry_auto_offset.Quantity
canonical_def = Q_(-161, "dBm") / module_registry_auto_offset.Hz
quantity_def = Q_(-161, "dBm/Hz")
assert quantity_def == canonical_def
def test_compound_log_unit_parse_definition(module_registry_auto_offset):
Q_ = module_registry_auto_offset.Quantity
canonical_def = Q_(-161, "dBm") / module_registry_auto_offset.Hz
parse_def = module_registry_auto_offset("-161 dBm/Hz")
assert parse_def == canonical_def
def test_compound_log_unit_parse_expr(module_registry_auto_offset):
"""Check that compound log units can be defined using ``parse_expression()``."""
Q_ = module_registry_auto_offset.Quantity
canonical_def = Q_(-161, "dBm") / module_registry_auto_offset.Hz
parse_def = module_registry_auto_offset.parse_expression("-161 dBm/Hz")
assert canonical_def == parse_def
@pytest.mark.xfail
def test_dbm_db_addition(module_registry_auto_offset):
"""Test a dB value can be added to a dBm and the answer is correct."""
power = (5 * module_registry_auto_offset.dBm) + (
10 * module_registry_auto_offset.dB
)
assert power.to("dBm").magnitude == pytest.approx(15)
@pytest.mark.xfail
@pytest.mark.parametrize(
"freq1,octaves,freq2",
[
(100, 2.0, 400),
(50, 1.0, 100),
(200, 0.0, 200),
], # noqa: E231
)
def test_frequency_octave_addition(module_registry_auto_offset, freq1, octaves, freq2):
"""Test an Octave can be added to a frequency correctly"""
freq1 = freq1 * module_registry_auto_offset.Hz
shift = octaves * module_registry_auto_offset.Octave
new_freq = freq1 + shift
assert new_freq.units == freq1.units
assert new_freq.magnitude == pytest.approx(freq2)
|