gwpy 3.0.8__py3-none-any.whl → 3.0.10__py3-none-any.whl

gwpy/_version.py

@@ -12,5 +12,5 @@ __version__: str
 __version_tuple__: VERSION_TUPLE
 version_tuple: VERSION_TUPLE
 
-__version__ = version = '3.0.8'
-__version_tuple__ = version_tuple = (3, 0, 8)
+__version__ = version = '3.0.10'
+__version_tuple__ = version_tuple = (3, 0, 10)

gwpy/astro/range.py

@@ -25,7 +25,7 @@ import warnings
 from functools import wraps
 from math import pi
 
-from scipy.integrate import trapz
+from scipy.integrate import trapezoid
 from scipy.interpolate import interp1d
 
 from astropy import (
@@ -238,7 +238,7 @@ def sensemon_range(psd, snr=8, mass1=1.4, mass2=1.4, fmin=None, fmax=None,
     integrand = sensemon_range_psd(psd[frange], snr=snr, mass1=mass1,
                                    mass2=mass2, horizon=horizon)
     return (units.Quantity(
-        trapz(integrand.value, f.value[frange]),
+        trapezoid(integrand.value, f.value[frange]),
         unit=integrand.unit * units.Hertz,
     ) ** (1/2.)).to('Mpc')
 
@@ -528,7 +528,7 @@ def burst_range(psd, snr=8, energy=1e-2, fmin=100, fmax=500):
     # calculate integrand and integrate
     integrand = burst_range_spectrum(
         psd[frange], snr=snr, energy=energy) ** 3
-    out = trapz(integrand.value, f[frange])
+    out = trapezoid(integrand.value, f[frange])
     # normalize and return
     return (units.Quantity(
         out / (fmax - fmin),

gwpy/astro/tests/test_range.py

@@ -24,7 +24,7 @@ from unittest import mock
 import pytest
 
 from astropy import units
-from scipy.integrate import trapz
+from scipy.integrate import trapezoid
 
 from ... import astro
 from ...testing import utils
@@ -71,7 +71,7 @@ def test_sensemon_range_psd(psd):
     r = astro.sensemon_range_psd(psd[frange])
     assert isinstance(r, FrequencySeries)
     utils.assert_quantity_almost_equal(
-        trapz(r, r.frequencies) ** (1/2.),
+        trapezoid(r, r.frequencies) ** (1/2.),
         TEST_RESULTS['sensemon_range'],
     )
     assert r.f0.value > 0
@@ -103,9 +103,8 @@ def test_inspiral_range_psd(psd):
     frange = (psd.frequencies.value < 4096)
     r = astro.inspiral_range_psd(psd[frange])
     assert isinstance(r, FrequencySeries)
-    print(trapz(r, r.frequencies) ** (1/2.))
     utils.assert_quantity_almost_equal(
-        trapz(r, r.frequencies) ** (1/2.),
+        trapezoid(r, r.frequencies) ** (1/2.),
         TEST_RESULTS['inspiral_range'],
     )
     assert r.f0.value > 0
@@ -129,7 +128,7 @@ def test_burst_range_spectrum(psd):
     r = astro.burst_range_spectrum(psd[frange])
     assert isinstance(r, FrequencySeries)
     utils.assert_quantity_almost_equal(
-        (trapz(r**3, f[frange]) / (400 * units.Hz)) ** (1/3.),
+        (trapezoid(r**3, f[frange]) / (400 * units.Hz)) ** (1/3.),
         TEST_RESULTS['burst_range'],
     )
     assert r.f0.value > 0

gwpy/cli/gwpy_plot.py

@@ -55,7 +55,7 @@ Examples:
     $ gwpy-plot coherencegram --chan H1:GDS-CALIB_STRAIN H1:PEM-CS_ACC_PSL_PERISCOPE_X_DQ --start 1126260017 --duration 600
 
 Written by {__author__}.
-Report bugs to https://github.com/gwpy/gwpy/issues/.
+Report bugs to https://gitlab.com/gwpy/gwpy/-/issues/.
 """  # noqa: E501
 
 

gwpy/cli/qtransform.py

@@ -215,7 +215,7 @@ class Qtransform(Spectrogram):
             outseg = Segment(gps, gps).protract(args.plot[self.plot_num])
 
             # use the precomputed ASD as the whitener if needed
-            if self.qxfrm_args.get("whiten"):
+            if self.qxfrm_args.get("whiten") is True:
                 self.qxfrm_args["whiten"] = asd
 
             # This section tries to optimize the amount of data that is

gwpy/cli/tests/base.py

@@ -19,6 +19,7 @@
 """Unit tests for :mod:`gwpy.cli`
 """
 
+import warnings
 from argparse import ArgumentParser
 from unittest import mock
 
@@ -237,7 +238,16 @@ class _TestCliProduct(object):
         xmin = min(series.xspan[0] for series in data)
         xmax = max(series.xspan[1] for series in data)
 
-        plotprod.set_plot_properties()
+        # ignore warnings from matplotlib about having no labels
+        # (because we have cut some corners in preparing this test)
+        with warnings.catch_warnings():
+            warnings.filterwarnings(
+                "ignore",
+                message="No artists with labels found to put in legend",
+                category=UserWarning,
+            )
+            plotprod.set_plot_properties()
+
         ax = plotprod.ax
 
         ymin, ymax = ax.get_ylim()

gwpy/detector/units.py

@@ -167,6 +167,7 @@ for unit, aliases in [
     # GW observatories like to record 'time' as the unit
     (units.Unit('second'), (
         'time',
+        'time (s)',
         'time [s]',
         'Time [sec]',
         'Time (sec)',
@@ -186,7 +187,7 @@ units.def_unit(['NONE', 'undef'], namespace=_ns,
 
 # other dimenionless units
 units.def_unit('strain', namespace=_ns)
-units.def_unit('coherence', namespace=_ns)
+units.def_unit('coherence', represents=units.dimensionless_unscaled, namespace=_ns)
 
 # alias for 'second' but with prefices
 units.def_unit((['sec'], ['sec']), represents=units.second, prefixes=True,

gwpy/frequencyseries/tests/test_hist.py

@@ -77,7 +77,8 @@ class TestSpectralVariance(_TestArray2D):
         utils.assert_array_equal(array.yindex, array.bins[:-1])
 
     def test_transpose(self, array):
-        return NotImplemented
+        with pytest.raises(NotImplementedError):
+            array.T
 
     # -- test utilities -------------------------
 
@@ -110,10 +111,10 @@ class TestSpectralVariance(_TestArray2D):
         pytest.skip("float precision test not supported for SpectralVariance")
 
     def test_is_compatible_yindex(self, array):
-        return NotImplemented
+        pytest.skip(f"not implemented for {type(self).__name__}")
 
     def test_is_compatible_error_yindex(self, array):
-        return NotImplemented
+        pytest.skip(f"not implemented for {type(self).__name__}")
 
     def test_plot(self, array):
         with rc_context(rc={'text.usetex': False}):

gwpy/io/datafind.py

@@ -65,6 +65,7 @@ GRB_TYPE = re.compile(r'^(?!.*_GRB\d{6}([A-Z])?$)')  # *_GRBYYMMDD{A}
 HIGH_PRIORITY_TYPE = re.compile("({})".format("|".join((
     r'\A[A-Z]\d_HOFT_C\d\d(_T\d{7}_v\d)?\Z',  # X1_HOFT_CXY
     r'\AV1Online\Z',
+    r'\AHoftOnline\Z',
     r'\AV1O[0-9]+([A-Z]+)?Repro[0-9]+[A-Z]+\Z',  # V1OXReproXY
 ))))
 LOW_PRIORITY_TYPE = re.compile("({})".format("|".join((

gwpy/io/ffldatafind.py

@@ -28,8 +28,9 @@ __author__ = "Duncan Macleod <duncan.macleod@ligo.org>"
 
 import os
 import re
-from warnings import warn
+from collections import defaultdict
 from functools import lru_cache
+from warnings import warn
 
 from ligo.segments import (
     segment,
@@ -37,8 +38,10 @@ from ligo.segments import (
 )
 
 from .cache import (
+    _CacheEntry,
     _iter_cache,
     cache_segments,
+    file_segment,
     read_cache_entry,
 )
 
@@ -123,10 +126,10 @@ def _find_ffl_files(basedir=None):
 def _find_ffls(basedir=None):
     """Find all readable FFL files.
     """
-    ffls = {}
+    ffls = defaultdict(list)
     for path in _find_ffl_files(basedir=basedir):
         try:
-            ffls[_get_site_tag(path)] = path
+            ffls[_get_site_tag(path)].append(path)
         except (
             OSError,  # file is empty (or cannot be read at all)
             AttributeError,  # last entry didn't match _SITE_REGEX
@@ -135,8 +138,8 @@ def _find_ffls(basedir=None):
     return ffls
 
 
-def _ffl_path(site, tag, basedir=None):
-    """Return the path of the FFL file for a given site and tag.
+def _ffl_paths(site, tag, basedir=None):
+    """Return the paths of all FFL files for a given site and tag.
     """
     try:
         return _find_ffls(basedir=basedir)[(site, tag)]
@@ -147,15 +150,18 @@ def _ffl_path(site, tag, basedir=None):
 
 
 @lru_cache()
-def _read_ffl(site, tag, basedir=None):
-    """Read an FFL file as a list of `CacheEntry` objects
+def _read_ffls(site, tag, basedir=None):
+    """Read all FFL files for a given site and tag
+    as a list of `CacheEntry` objects.
     """
-    ffl = _ffl_path(site, tag, basedir=basedir)
-    with open(ffl, "r") as fobj:
-        return [
-            type(entry)(site, tag, entry.segment, entry.path)
-            for entry in _iter_cache(fobj, gpstype=float)
-        ]
+    entries = []
+    for ffl in _ffl_paths(site, tag, basedir=basedir):
+        with open(ffl, "r") as fobj:
+            entries.extend(
+                _CacheEntry(site, tag, entry.segment, entry.path)
+                for entry in _iter_cache(fobj, gpstype=float)
+            )
+    return entries
 
 
 def _handle_error(action, message):
@@ -261,7 +267,7 @@ def find_urls(
     span = segment(gpsstart, gpsend)
 
     cache = [
-        e for e in _read_ffl(site, tag) if (
+        e for e in _read_ffls(site, tag) if (
             e.observatory == site
             and e.description == tag
             and e.segment.intersects(span)
@@ -305,11 +311,16 @@ def find_latest(site, tag, on_missing="warn"):
         for a specific site and tag.
     """
     try:
-        fflfile = _ffl_path(site, tag)
+        fflfiles = _ffl_paths(site, tag)
     except ValueError:  # no readable FFL file
         urls = []
     else:
-        urls = [read_cache_entry(_read_last_line(fflfile), gpstype=float)]
+        urls = [
+            read_cache_entry(_read_last_line(fflfile), gpstype=float)
+            for fflfile in fflfiles
+        ]
+        if urls:  # if multiple, find the latest one
+            urls = sorted(urls, key=file_segment)[-1:]
 
     if not urls:
         _handle_error(on_missing, "No files found")

gwpy/io/tests/test_ffldatafind.py

@@ -96,9 +96,11 @@ def test_is_ffl_file(path, result):
 # -- test ffl UI ------------
 
 FFLS = {
-    "test.ffl": [
+    "a/test.ffl": [
         "/tmp/X-test-0-1.gwf 0 1 0 0",
         "/tmp/X-test-1-1.gwf 1 1 0 0",
+    ],
+    "b/test.ffl": [
         "/tmp/X-test-2-1.gwf 2 1 0 0",
     ],
     "test2.ffl": [
@@ -114,7 +116,11 @@ FFLS = {
     "test-empty.ffl": [],
     "test-bad.ffl": ["badness"],
 }
-TEST_URLS = [x.split()[0] for x in FFLS["test.ffl"]]
+TEST_URLS = [
+    x.split()[0]
+    for key in ("a/test.ffl", "b/test.ffl")
+    for x in FFLS[key]
+]
 
 
 @pytest.fixture(autouse=True)
@@ -124,6 +130,7 @@ def mock_ffl(tmp_path):
     """
     for path, lines in FFLS.items():
         ffl = tmp_path / path
+        ffl.parent.mkdir(parents=True, exist_ok=True)
         ffl.write_text("\n".join(lines))
     with mock.patch.dict(
         "os.environ",
@@ -190,7 +197,7 @@ def test_find_latest():
     assert ffldatafind.find_latest(
         "X",
         "test",
-    ) == sorted(x.split()[0] for x in FFLS["test.ffl"])[-1:]
+    ) == sorted(x.split()[0] for x in FFLS["b/test.ffl"])[-1:]
 
 
 @pytest.mark.parametrize(("on_missing", "ctx"), (

gwpy/plot/plot.py

@@ -27,7 +27,12 @@ from itertools import zip_longest
 
 import numpy
 
-from matplotlib import (figure, get_backend, _pylab_helpers)
+from matplotlib import (
+    _pylab_helpers,
+    backends,
+    figure,
+    get_backend,
+)
 from matplotlib.artist import setp
 from matplotlib.gridspec import GridSpec
 from matplotlib.ticker import LogFormatterSciNotation
@@ -53,7 +58,12 @@ iterable_types = (list, tuple, KeysView, ValuesView,)
 def interactive_backend():
     """Returns `True` if the current backend is interactive
     """
-    from matplotlib.rcsetup import interactive_bk
+    try:
+        interactive_bk = backends.backend_registry.list_builtin(
+            backends.BackendFilter.INTERACTIVE,
+        )
+    except AttributeError:  # matplotlib < 3.9.0
+        from matplotlib.rcsetup import interactive_bk
     return get_backend() in interactive_bk
 
 

gwpy/plot/tests/test_segments.py

@@ -136,7 +136,7 @@ class TestSegmentAxes(_TestAxes):
     # -- disable tests from upstream
 
     def test_imshow(self):
-        return NotImplemented
+        pytest.skip(f"not implemented for {type(self).__name__}")
 
 
 def test_segmentrectangle():

gwpy/plot/tex.py

@@ -39,10 +39,12 @@ def _test_usetex():
     from matplotlib import (pyplot, rc_context)
     with rc_context({"text.usetex": True}):
         fig = pyplot.figure()
-        ax = fig.gca()
-        ax.set_xlabel(r"\LaTeX")
-        fig.canvas.draw()
-    pyplot.close(fig)
+        try:
+            ax = fig.gca()
+            ax.set_xlabel(r"\LaTeX")
+            fig.canvas.draw()
+        finally:
+            pyplot.close(fig)
 
 
 @lru_cache(maxsize=None)

gwpy/segments/flag.py

@@ -453,18 +453,20 @@ class DataQualityFlag(object):
         Examples
         --------
         >>> from gwpy.segments import DataQualityFlag
-        >>> print(DataQualityFlag.fetch_open_data('H1_DATA', 'Jan 1 2010',
-        ...                                       'Jan 2 2010'))
+        >>> print(DataQualityFlag.fetch_open_data(
+        ...     "H1_DATA",
+        ...     "Sep 14 2015",
+        ...     "Sep 15 2015",
+        ... ))
         <DataQualityFlag('H1:DATA',
-                         known=[[946339215 ... 946425615)],
-                         active=[[946340946 ... 946351800)
-                                 [946356479 ... 946360620)
-                                 [946362652 ... 946369150)
-                                 [946372854 ... 946382630)
-                                 [946395595 ... 946396751)
-                                 [946400173 ... 946404977)
-                                 [946412312 ... 946413577)
-                                 [946415770 ... 946422986)],
+                         known=[[1126224017 ... 1126310417)]
+                         active=[[1126251604 ... 1126252133)
+                                 [1126252291 ... 1126274322)
+                                 [1126276234 ... 1126281754)
+                               ...
+                                 [1126308670 ... 1126309577)
+                                 [1126309637 ... 1126309817)
+                                 [1126309877 ... 1126310417)]
                          description=None)>
         """
         start = to_gps(start).gpsSeconds

gwpy/signal/filter_design.py

@@ -30,6 +30,7 @@ from scipy import signal
 
 from astropy.units import (Unit, Quantity)
 
+from ..types.array import COPY_IF_NEEDED
 from .window import (get_window, planck)
 
 __author__ = "Duncan Macleod <duncan.macleod@ligo.org>"
@@ -88,8 +89,8 @@ def _design_fir(wp, ws, sample_rate, gpass, gstop, window='hamming', **kwargs):
     if wp[0] > ws[0]:
         kwargs.setdefault('pass_zero', False)
     if ws.shape == (1,):
-        kwargs.setdefault('width', ws - wp)
-    kwargs.setdefault('nyq', sample_rate/2.)
+        kwargs.setdefault("width", ws.item() - wp.item())
+    kwargs.setdefault('fs', sample_rate)
     return signal.firwin(nt, wp, window=window, **kwargs)
 
 
@@ -281,9 +282,9 @@ def bilinear_zpk(zeros, poles, gain, fs=1.0, unit='Hz'):
     zpk : `tuple`
         digital version of input zpk
     """
-    zeros = numpy.array(zeros, dtype=float, copy=False)
+    zeros = numpy.array(zeros, dtype=float, copy=COPY_IF_NEEDED)
     zeros = zeros[numpy.isfinite(zeros)]
-    poles = numpy.array(poles, dtype=float, copy=False)
+    poles = numpy.array(poles, dtype=float, copy=COPY_IF_NEEDED)
     gain = gain
 
     # convert from Hz to rad/s if needed

gwpy/signal/tests/test_coherence.py

@@ -21,6 +21,8 @@
 
 __author__ = "Alex Southgate <alex.southgate@ligo.org>"
 
+import warnings
+
 import numpy as np
 import pytest
 import scipy.signal as sig
@@ -121,7 +123,11 @@ def test_coherence_resample(unequal_fs_series_data):
     # the first coherence val coh12 is broken intentionally since
     # secondarr data should not have fs_1, instead fs_2
     coh12 = spectral.coherence(first, second, segmentlength=seglen)
-    coh13 = spectral.coherence(first, third, segmentlength=seglen)
+    with pytest.warns(
+        UserWarning,
+        match="Sampling frequencies are unequal",
+    ):
+        coh13 = spectral.coherence(first, third, segmentlength=seglen)
 
     # get the frequency at minimum coherence, this should be the extra
     # component in secondarr
@@ -136,7 +142,7 @@ def test_coherence_resample(unequal_fs_series_data):
     assert not (4 <= maxf12.value <= 6)
 
 
-def test_coherence_resample_arg(series_data):
+def test_coherence_resample_downsample(series_data):
     """Ensure warning is raised by unequal sampling frequencies.
     """
     firstarr, secondarr, seglen = series_data
@@ -145,15 +151,30 @@ def test_coherence_resample_arg(series_data):
     first = TimeSeries(firstarr, sample_rate=f_s)
     second = TimeSeries(secondarr, sample_rate=f_s * 2.32)
 
-    with pytest.warns(UserWarning, match="Sampling frequencies are unequal"):
-        spectral.coherence(first, second, segmentlength=seglen)
+    with pytest.warns(
+        UserWarning,
+        match="Sampling frequencies are unequal",
+    ):
+        coh1 = spectral.coherence(first, second, segmentlength=seglen)
 
+    # check that forcibly disabling downsample results in an error
     with pytest.raises(ValueError):
-        spectral.coherence(first, second, segmentlength=seglen,
-                           downsample=False)
-
-    coh1 = spectral.coherence(first, second, segmentlength=seglen)
-    coh2 = spectral.coherence(first, second, segmentlength=seglen,
-                              downsample=True)
+        spectral.coherence(
+            first,
+            second,
+            segmentlength=seglen,
+            downsample=False,
+        )
+
+    # but that accepting downsampling gives you the same result as
+    # doing nothing (but doesn't emit a warning)
+    with warnings.catch_warnings():
+        warnings.simplefilter("error", DeprecationWarning)
+        coh2 = spectral.coherence(
+            first,
+            second,
+            segmentlength=seglen,
+            downsample=True,
+        )
 
     assert all(np.array(coh1.data) == np.array(coh2.data))

gwpy/signal/tests/test_filter_design.py

@@ -50,7 +50,7 @@ LOWPASS_IIR_100HZ = signal.iirdesign(
     analog=False, ftype='cheby1', output='zpk',
 )
 LOWPASS_FIR_100HZ = signal.firwin(
-    30, 100, window='hamming', width=50., nyq=512.,
+    30, 100, window='hamming', width=50., fs=1024.,
 )
 
 HIGHPASS_IIR_100HZ = signal.iirdesign(
@@ -60,7 +60,7 @@ HIGHPASS_IIR_100HZ = signal.iirdesign(
     analog=False, ftype='cheby1', output='zpk',
 )
 HIGHPASS_FIR_100HZ = signal.firwin(
-    45, 100, window='hamming', pass_zero=False, width=-100/3., nyq=512.,
+    45, 100, window='hamming', pass_zero=False, width=-100/3., fs=1024.,
 )
 
 BANDPASS_IIR_100HZ_200HZ = signal.iirdesign(
@@ -70,7 +70,7 @@ BANDPASS_IIR_100HZ_200HZ = signal.iirdesign(
     analog=False, ftype='cheby1', output='zpk',
 )
 BANDPASS_FIR_100HZ_200HZ = signal.firwin(
-    45, (100, 200.), window='hamming', pass_zero=False, nyq=512.,
+    45, (100, 200.), window='hamming', pass_zero=False, fs=1024.,
 )
 
 

gwpy/spectrogram/tests/test_spectrogram.py

@@ -175,11 +175,11 @@ class TestSpectrogram(_TestArray2D):
 
         # test simple filter
         a2 = array.filter(*zpk)
-        utils.assert_array_equal(array * fresp, a2)
+        utils.assert_quantity_sub_equal(array * fresp, a2)
 
         # test inplace filtering
         array.filter(lti, inplace=True)
-        utils.assert_array_equal(array, a2)
+        utils.assert_quantity_sub_equal(array, a2)
 
         # test errors
         with pytest.raises(TypeError):

gwpy/table/io/ligolw.py

@@ -58,7 +58,7 @@ try:
 except ImportError:
     pass
 else:
-    NUMPY_TYPE_MAP[LIGOTimeGPS] = numpy.float_
+    NUMPY_TYPE_MAP[LIGOTimeGPS] = numpy.float64
 
 
 # -- utilities ----------------------------------------------------------------

gwpy/table/tests/test_io_ligolw.py

@@ -70,7 +70,7 @@ def test_to_astropy_table_empty():
     tab = EventTable(llwtable, columns=["peak", "ifo"])
     assert set(tab.colnames) == {"peak", "ifo"}
     assert tab['peak'].dtype.type is numpy.object_
-    assert tab['ifo'].dtype.type is numpy.unicode_
+    assert tab['ifo'].dtype.type is numpy.str_
 
 
 @pytest.mark.requires("ligo.lw.lsctables")

gwpy/testing/errors.py

@@ -31,6 +31,7 @@ import requests.exceptions
 NETWORK_ERROR = (
     ConnectionError,
     requests.exceptions.ConnectionError,
+    requests.exceptions.ReadTimeout,
     socket.timeout,
     SSLError,
     URLError,

gwpy/testing/fixtures.py

@@ -52,25 +52,8 @@ def tmpfile():
 
 # -- plotting -----------------------------------------------------------------
 
-def _test_usetex():
-    """Return `True` if we can render figures using `text.usetex`.
-    """
-    from matplotlib import pyplot
-    with rc_context(rc={'text.usetex': True}):
-        fig = pyplot.figure()
-        fig.gca()
-        try:
-            fig.canvas.draw()
-        except RuntimeError:
-            return False
-        else:
-            return True
-        finally:
-            pyplot.close(fig)
-
-
 SKIP_TEX = pytest.mark.skipif(
-    not has_tex() or not _test_usetex(),
+    not has_tex(),
     reason='TeX is not available',
 )
 

gwpy/testing/utils.py

@@ -34,12 +34,14 @@ from numpy.testing import (assert_array_equal, assert_allclose)
 
 from astropy.time import Time
 
-from ..utils.decorators import deprecated_function
+from gwpy.io.cache import file_segment
+from gwpy.utils.decorators import deprecated_function
 
 # -- useful constants ---------------------------------------------------------
 
 TEST_DATA_DIR = os.path.join(os.path.dirname(__file__), 'data')
 TEST_GWF_FILE = os.path.join(TEST_DATA_DIR, 'HLV-HW100916-968654552-1.gwf')
+TEST_GWF_SPAN = file_segment(TEST_GWF_FILE)
 TEST_HDF5_FILE = os.path.join(TEST_DATA_DIR, 'HLV-HW100916-968654552-1.hdf')
 
 
@@ -328,11 +330,18 @@ def TemporaryFilename(*args, **kwargs):  # pylint: disable=invalid-name
             os.remove(name)
 
 
-def test_read_write(data, format,
-                    extension=None, autoidentify=True,
-                    read_args=[], read_kw={},
-                    write_args=[], write_kw={},
-                    assert_equal=assert_array_equal, assert_kw={}):
+def test_read_write(
+    data,
+    format,
+    extension=None,
+    autoidentify=True,
+    read_args=[],
+    read_kw={},
+    write_args=[],
+    write_kw={},
+    assert_equal=assert_quantity_sub_equal,
+    assert_kw={},
+):
     """Test that data can be written to and read from a file in some format
 
     Parameters

gwpy/timeseries/core.py

@@ -793,7 +793,6 @@ class TimeSeriesBase(Series):
                 op_ = ufunc.__name__
             out = out.view(StateTimeSeries)
             out.__metadata_finalize__(orig)
-            out.override_unit('')
             oname = orig.name if isinstance(orig, type(self)) else orig
             vname = value.name if isinstance(value, type(self)) else value
             out.name = '{0!s} {1!s} {2!s}'.format(oname, op_, vname)