Pith. sign in

REVIEW 2 major objections 5 minor 1 cited by

An open, differentiable WDM transform matches frequency-domain LISA posteriors under stationary noise.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · grok-4.5

2026-07-12 13:15 UTC pith:VP6JG5WP

load-bearing objection Solid methods-and-software paper: explicit discrete WDM reference plus a working NumPy/JAX package that matches FD posteriors under stationary noise; the non-stationary claims are correctly deferred. the 2 major comments →

arxiv 2606.20269 v2 pith:VP6JG5WP submitted 2026-06-18 gr-qc astro-ph.HE

An explicit and differentiable Wilson-Daubechies-Meyer transform for gravitational-wave data analysis

Avi Vajpeyi , Giorgio Mentasti , Quentin Baghi , Ollie Burke , Lorenzo Speri This is my paper
classification gr-qc astro-ph.HE
keywords Wilson-Daubechies-Meyer transformgravitational-wave data analysistime-frequency methodsWDM likelihoodLISA galactic binariesJAX autodiffwavelet packets
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Gravitational-wave analysts often need a time-frequency basis that is orthonormal, critically sampled, and usable inside modern autodiff pipelines. This paper supplies a fully explicit derivation of the Wilson–Daubechies–Meyer (WDM) wavelet-packet transform—covering DC, interior, and Nyquist channels, the packed coefficient layout, and the forward and inverse maps—and ships it as a NumPy/JAX package that reconstructs data to floating-point precision. Under a shared stationary noise model, the resulting WDM-domain Whittle likelihood reproduces frequency-domain Bayesian posteriors for resolved LISA galactic binaries to numerical precision across one hundred independent simulations. The practical payoff is a ready-to-use, GPU-accelerated transform that can be dropped into differentiable inference workflows and later extended to non-stationary noise, gaps, and tiling optimisation for future detectors.

Core claim

The authors show that a carefully specified, packed WDM transform, implemented with exact round-trip reconstruction, yields a WDM-domain likelihood that is numerically equivalent to the ordinary frequency-domain Whittle likelihood for resolved LISA galactic binaries when both analyses share the same stationary noise model. Across 100 independent injections the median Jensen–Shannon divergence between one-dimensional marginals is about 5e-6 bits, confirming that the two representations agree to sampling precision in that controlled setting.

What carries the argument

The packed Wilson–Daubechies–Meyer transform: an orthonormal wavelet-packet map that sends N real samples to an (Nt, Nf+1) coefficient grid (DC, interior, and Nyquist channels), with Meyer-windowed atoms, Wilson phase factors, and exact forward/inverse formulae that make the Gram matrix G†G equal the identity, enabling a diagonal WDM Whittle likelihood under local stationarity.

Load-bearing premise

The claim that noise can be treated as locally stationary so that the WDM covariance is diagonal and off-diagonal pixel correlations can be dropped; if the spectrum changes faster than the tiling scales, that diagonal likelihood fails.

What would settle it

Re-run the same LISA galactic-binary inference with a deliberately non-stationary, time-varying PSD and check whether the WDM and frequency-domain posteriors still agree to the same Jensen–Shannon tolerance; any large systematic offset would falsify the claimed numerical equivalence outside the stationary regime.

Watch this falsifier — get emailed when new claim-graph text bears on it.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

2 major / 5 minor

Summary. The manuscript supplies a self-contained derivation and open-source NumPy/JAX implementation of the discrete Wilson–Daubechies–Meyer (WDM) wavelet-packet transform for gravitational-wave data analysis. It fixes the packed (Nt, Nf+1) coefficient layout, the three-case Wilson basis (DC, interior, Nyquist), Meyer-window conventions, and forward/inverse formulae; proves orthonormality and the stationary diagonal-covariance approximation in the appendices; and validates floating-point-exact round-trip reconstruction. As a controlled end-to-end check, Sec. IV A shows that a WDM-domain Whittle likelihood reproduces frequency-domain posteriors for resolved LISA galactic binaries under a shared stationary noise model (median JSD ~5e-6 bits over 100 seeds).

Significance. The contribution is primarily infrastructural and pedagogical rather than a new transform or a new science result. Its value lies in making an already-used representation fully explicit, reproducible, and differentiable for modern Python/JAX inference pipelines, with GPU-accelerated million-point transforms and machine-precision round-trips. The LISA comparison is a clean numerical equivalence test under the regime where the diagonal approximation is expected to hold. The package, Zenodo archive, and pinned code links are genuine strengths that lower the barrier to adoption and to future work on non-stationary noise, gaps, and tiling optimisation that the paper itself defers.

major comments (2)
  1. The abstract and introduction advertise utility for non-stationary noise, stochastic backgrounds, and data gaps, yet the only end-to-end inference demonstration (Sec. IV A) uses a shared time-averaged stationary PSD in both domains. That choice correctly isolates the basis change and makes the tiny JSD expected, but it does not exercise the local-stationarity regime of Eq. (38)/App. C that is the principal statistical motivation for WDM. The paper already acknowledges this in Sec. V; the abstract and conclusions should be tightened so that the demonstrated claim (stationary numerical equivalence) is not oversold relative to the deferred non-stationary applications.
  2. Sec. III A–B and App. C derive the WDM Whittle likelihood and inner product after dropping DC/Nyquist channels and invoking the diagonal approximation (Eq. 38). The matrix formalism of Sec. II D correctly notes that Σw is rank-deficient when edge channels are retained. The manuscript should state more explicitly, near Eqs. (40)–(42), the quantitative conditions under which the interior-only diagonal form remains adequate for inference (e.g., spectral content near DC/Nyquist, or a bound on residual off-diagonal power after whitening), so that users know when the simple pixel-wise likelihood is safe.
minor comments (5)
  1. Notation for the Meyer taper parameter switches between A (main text) and a (API/package); a single convention would reduce friction for implementers.
  2. Fig. 2 caption and surrounding text are dense; a short sentence stating that panel (a) lives in the packed Nt(Nf+1) space while panel (b) lives in the N-dimensional signal space would help first-time readers.
  3. The evenness requirement on Nt and Nf is stated but not motivated beyond “left to future work”; a one-sentence remark on why the current phase/FFT construction needs it would be useful.
  4. A few typographical inconsistencies remain (e.g., “T erminology”, “F orward”, “ST A TISTICS”, stray /github markers in the source).
  5. Related WDM implementations (Digman/Cornish, Moore, Farr, cWB) are cited; a short comparison table of layout conventions and window defaults would help practitioners choose among them.

Circularity Check

0 steps flagged

No circularity: orthonormality, transforms, and diagonal covariance are derived from the basis definitions; the FD–WDM posterior match is an independent numerical check under shared stationary assumptions.

full rationale

The paper’s load-bearing mathematical claims are self-contained derivations rather than tautologies. Appendix A proves the two orthonormality relations (Eq. 16) directly from the three-case Wilson basis (Eq. 11) and the d=1 Meyer window (Eq. 12), using only the partition-of-unity property of the cosine taper and geometric-series identities; the target result is not assumed. Appendix B derives the packed forward coefficients (Eqs. 18–23) and inverse (Eq. 26) by substituting the same basis into the analysis formula and applying reality/symmetry of the DFT, again without circular substitution. Appendix C obtains the diagonal WDM covariance (Eq. 38) as a first-order approximation under wide-sense stationarity by pulling a slowly varying PSD out of the integral against the localized atoms; the approximation is stated with its failure mode (rapid PSD variation relative to ΔT, ΔF) and is not redefined as exact. The Section IV A LISA galactic-binary comparison evaluates two likelihoods (frequency-domain Whittle vs. WDM Whittle) on identical simulated data, source model, priors, and stationary PSD; the reported median JSD ~5e-6 bits is an empirical numerical equivalence check, not a quantity forced by fitting or by construction. Prior WDM literature is cited as historical background and algorithmic precedent; the present work does not import a uniqueness theorem or ansatz from overlapping authors to forbid alternatives, nor does it rename a known empirical pattern as a new prediction. Round-trip reconstruction to floating-point precision is an independent software validation of the packed layout. Consequently there are no self-definitional loops, fitted-input-as-prediction steps, or load-bearing self-citation chains.

Axiom & Free-Parameter Ledger

3 free parameters · 4 axioms · 0 invented entities

The central claims rest on standard Fourier analysis, the classical Wilson/Daubechies–Meyer construction, and the usual wide-sense-stationary noise model used in GW Whittle likelihoods. Free parameters are user-chosen tiling and window knobs (Nt, A, d), not fitted to data to force the result. No new physical entities are postulated; the package is an implementation of an existing basis.

free parameters (3)
  • Nt (number of time bins) = example Nt=32 (LISA demo); benchmarks use Nt=1024
    User-chosen tiling parameter that sets ΔT and ΔF; must be even. Affects localization trade-off but is not fitted to produce the FD–WDM match.
  • A (Meyer taper width) = default 1/3
    Controls passband/roll-off of the frequency window; default A=1/3, admissible in (0,1/2). Exposed as API parameter a.
  • d (Meyer order) = 1 (fixed in current package)
    Smoothness order of the taper; implementation fixes d=1 (cosine) while retaining d in the API for future work. Cornish used d=4.
axioms (4)
  • standard math Discrete Fourier transform conventions and reality condition for real time series (numpy-style, zero-centered frequencies).
    Section I A; used throughout forward/inverse derivations.
  • domain assumption Wilson basis pairing of positive/negative frequency atoms with phase C_nm yields a real orthogonal packet basis when combined with a Meyer window.
    Taken from Daubechies–Jaffard–Journé / prior WDM GW literature; proved for the discrete packed case in Appendix A.
  • domain assumption Wide-sense stationary (or locally stationary) Gaussian noise with slowly varying PSD relative to the WDM pixel support, so the covariance is approximately diagonal in the WDM domain.
    Appendix C and Section III; underpins the WDM Whittle likelihood used in the LISA example.
  • ad hoc to paper Nt and Nf are even positive integers with N = Nt Nf.
    Stated as a current implementation requirement; relaxing evenness is left to future work.

pith-pipeline@v1.1.0-grok45 · 31550 in / 3041 out tokens · 24160 ms · 2026-07-12T13:15:42.555593+00:00 · methodology

0 comments
read the original abstract

The Wilson-Daubechies-Meyer (WDM) time-frequency transform has been widely used in gravitational-wave astronomy, yet a self-contained, mathematically explicit reference for practitioners remains lacking. This is especially true for those wishing to adopt the transform in modern Python and JAX inference workflows. We present wdm_transform, an open-source Python package implementing the WDM wavelet-packet time-frequency transform, and document its mathematical foundations, statistical properties, and practical implementation for gravitational-wave data analysis. The package supplies NumPy and JAX backends, both transforms (forward and inverse) validated to floating-point precision, with the JAX backend enabling GPU-accelerated transforms of million-point data streams in tens of milliseconds. As a worked example, we verify that the WDM-domain likelihood reproduces frequency-domain posteriors for a resolved LISA galactic binary under a shared stationary noise model, confirming numerical equivalence of the two representations in that controlled setting. This work paves the way for systematic optimisation of WDM tilings, a particularly promising direction for the non-stationary noise, stochastic backgrounds, and data gaps anticipated in future detectors, and for direct comparisons with alternative time-frequency representations needed to meet the challenges of future gravitational-wave data analysis.

Figures

Figures reproduced from arXiv: 2606.20269 by Avi Vajpeyi, Giorgio Mentasti, Lorenzo Speri, Ollie Burke, Quentin Baghi.

Figure 1
Figure 1. Figure 1: Schematic illustration of the WDM forward transform. (a) Time-domain atoms [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Orthogonality structure and empirical decorrelation of the sampled WDM transform. Both panels use the small tiling Nt = Nf = 8, N = NtNf = 64, with packed atoms flattened by α = n(Nf + 1) + m. Panel (a) lives in the Nt(Nf + 1) = 72-dimensional packed coefficient space, while panel (b) lives in the N-dimensional signal space. The asymmetry GG† ∈ C 72×72 versus G†G ∈ C 64×64 is intrinsic to the rectangular e… view at source ↗
Figure 3
Figure 3. Figure 3: Architecture of the wdm transform package. The public API exposes the TimeSeries, FrequencySeries, and WDM datatypes with reversible conversions between them. The forward and inverse WDM transforms, together with the shared window machinery (windows.py), are dispatched to interchangeable NumPy, CuPy, and JAX backends, so the same high-level code runs on CPU, GPU, or TPU. § efficiency improves and, at N = 1… view at source ↗
Figure 4
Figure 4. Figure 4: Forward WDM runtime versus input length. All transforms use a fixed tiling with Nt = 1024 time bins and Nf = N/Nt frequency channels. Top: me￾dian single-transform wall-clock runtime over 7 runs af￾ter one warmup call. The thin solid black line shows the N log2 N reference scaling. Bottom: batched-execution speedup (WDM transform only) tserial/tbatch for a batch size of B = 3, where the serial baseline app… view at source ↗
Figure 5
Figure 5. Figure 5: Injected resolved galactic binary in both fre￾quency and WDM domains. (a): channel-A frequency￾domain data, the instrumental noise PSD, and the injected source, with the analysis band shaded. Bottom: zooms of the source band in the frequency domain (left (b)) and in the whitened WDM time-frequency plane (right (c)), where amplitude is shown in units of the per-pixel noise standard deviation. § the conditio… view at source ↗
Figure 7
Figure 7. Figure 7: Population PP plot. Empirical distribution of the injected truth’s posterior quantile across the 100 seeds, for each parameter, in the WDM domain (solid, lower opac￾ity) and frequency domain (dashed). Grey shading shows the nested 1/2/3σ pointwise confidence bands expected under the uniform null. The two domains track each other closely and every parameter remains within the 3σ band. The small￾est per-para… view at source ↗

Review history (2 revisions) →

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. Modeling non-stationary noise: applications in gravitational wave astronomy

    gr-qc 2026-07 conditional novelty 6.0

    A positive dynamic spectrum S(f,t) generalizes the stationary power spectrum by defining Gramian closed-form noise covariances in Fourier and Wilson-Daubechies wavelet bases for gravitational wave data.

Reference graph

Works this paper leans on

50 extracted references · 20 linked inside Pith · cited by 1 Pith paper

  1. [1]

    Wil- son

    coefficient layout explicitly, spell out the DC and Nyquist channel treatment, and set the normalization, phase, and indexing conventions used by the implemen- tation. These choices are then validated through ex- act round-trip reconstruction tests – applying the for- ward transform followed by its inverse – and empiri- cal checks of coefficient decorrela...

  2. [2]

    (10) with the interior-channel basis element from Eq

    Interior channels(0< m < N f) Starting from Eq. (10) with the interior-channel basis element from Eq. (11), the forward transform can be 3 Our definitions of the Wilson basis functions (11) carefully in- clude them= 0 andm=N f in order for orthogonality to hold. evaluated in the Fourier domain using wnm = N/2−1X l=−N/2 ˜x[l] ˜g∗ nm[l],(17) For a real-valu...

  3. [3]

    DC channel(m= 0) The DC edge channel has a doubled time-shift ex- ponent and noC nm factor in the basis definition (11). 7 Following the same Fourier-domain evaluation yields wn0 = √ 2ℜ   Nt/2−1X l=1 ˜x[l]e4πiln/Nt ˜φ[l]   + 1√ 2 ˜x[0] ˜φ[0]./github (21) Note the factor of 4πin the exponent (versus 2πfor interior channels): the DC basis element oscill...

  4. [4]

    Towards LISA catalogs

    Nyquist channel(m=N f) The Nyquist edge channel is handled analogously to the DC channel but centered at the Nyquist frequency fmax. Its structure mirrors the DC case with appropri- ate frequency shifts. In practice, the Nyquist channel is computed using the same windowed-FFT machinery as the interior channels, with the frequency shift set to mNt/2 =N/2. ...

    2023

  5. [5]

    Now we have to consider three cases

    The second (quasi)normality condition We want to evaluate X l ˜gnm[l]˜g∗ pq[l] (A10) we immediately note that, since the window support in frequency isN t/2 we are forced to haveq=m. Now we have to consider three cases. Casem= 0, for which we have X l ˜gn0[l]˜g∗ p0[l] = 1 2 N/2−1X l=−N/2 e−4πil(n−p)/Nt ˜φ2[l] = 1 2 Nt/2−1X l=−Nt/2 e−4πil(n−p)/Nt ˜φ2[l] = ...

  6. [6]

    (A2) is correctly normalized

    Normalization of the Meyer window We verify that the Meyer window defined in Eq. (A2) is correctly normalized. We compute N/2−1X l=−N/2 ˜φ2[l] = 2 A Nt 2 −1X l=0 2 Nt − 2 Nt + 2 (A+B) Nt 2 −1X l=A Nt 2 2 Nt cos2 " π 2 2l Nt −A B # = 2A− 2 Nt + 4 Nt B Nt 2 −1X k=0 cos2 π 2 2k BNt = = 2A− 2 Nt + 4 Nt B Nt 2 −1X k=0 1 2 + 4 Nt B Nt 2 −1X k=0 1 2 cos 2πk BNt ...

  7. [7]

    Interior channels(0< m < N f) Inserting the interior-channel basis element from Eq. (11) and splitting the sum over positive and negative frequencies gives wnm = 1√ 2 N/2−1X l=0 ˜x[l]e2πiln/Nt C ∗ nm ˜φ l− mNt 2 +C nm ˜φ l+ mNt 2 + 1√ 2 −1X l=−N/2 ˜x[l]e2πiln/Nt C ∗ nm ˜φ l− mNt 2 +C nm ˜φ l+ mNt 2 .(B2) Using the reality condition ˜x[−l] = ˜x[l]∗ togethe...

  8. [8]

    DC channel(m= 0) For the DC edge channel, the basis element has a doubled time-shift exponent and noC nm factor. The same positive/negative frequency splitting yields wn0 = 1√ 2 N/2−1X l=0 ˜x[l]e4πiln/Nt ˜φ[l] +1√ 2 −1X l=−N/2 ˜x[l]e4πiln/Nt ˜φ[l] = √ 2ℜ   Nt/2−1X l=1 ˜x[l]e4πiln/Nt ˜φ[l]   + 1√ 2 ˜x[0] ˜φ[0],(B4) which is Eq. (21) of the main text

  9. [9]

    Nyquist channel(m=N f) wnm = 1√ 2 N/2−1X l=0 ˜x[l]e4πiln/Nt ˜φ l− N 2 + ˜φ l+ N 2 + 1√ 2 −1X l=−N/2 ˜x[l]e4πiln/Nt ˜φ l− N 2 + ˜φ l+ N 2 = 1√ 2 N/2−1X l=0 ˜x[l]e4πiln/Nt ˜φ l− N 2 + ˜φ l+ N 2 + 1√ 2 N/2X l=1 ˜x[−l]e−4πiln/Nt ˜φ −l− N 2 + ˜φ −l+ N 2 .(B5) Thel= 0 term in the first summation vanishes, and thel=N/2 term is treated separately: wnm = √ 2ℜ N/2−...

  10. [10]

    Appendix C: The WDM Covariance Matrix in the stationary case We start from the WDM transform in Eq

    The inverse transform We write the inverse transform of thew nm coefficients to recover the signal ˜x[l] forl≥0 since we know that for l <0 the reality condition holds (˜x[−l] = ˜x∗[l]): ˜x[l] = Nt−1X n=0 NfX m=0 wnmgnm[l] = 1√ 2 X n X 0<m<Nf wnme−2πiln/Nt Cnm ˜φ l− mNt 2 +C ∗ nm ˜φ l+ mNt 2 + 1√ 2 X n wn0e−4πiln/Nt ˜φ[l] +1√ 2 X n wnNf e−4πiln/Nt ˜φ l− N...

  11. [11]

    Kr´ olak and P

    A. Kr´ olak and P. Trzaskoma, Classical and Quantum Gravity13, 813 (1996)

    1996

  12. [12]

    Virtuoso and E

    A. Virtuoso and E. Milotti, Phys. Rev. D109, 102010 (2024)

    2024

  13. [13]

    Robinet, N

    F. Robinet, N. Arnaud, N. Leroy, A. Lundgren, D. Macleod, and J. McIver, SoftwareX12, 100620 (2020)

    2020

  14. [14]

    N. J. Cornish and T. B. Littenberg, Classical and Quantum Gravity32, 135012 (2015)

    2015

  15. [15]

    Chatterji, L

    S. Chatterji, L. Blackburn, G. Martin, and E. Katsavounidis, Class. Quant. Grav.21, S1809 (2004), arXiv:gr-qc/0412119

    Pith/arXiv arXiv 2004

  16. [16]

    Chassande-Mottin and A

    E. Chassande-Mottin and A. Pai, Phys. Rev. D73, 042003 (2006), arXiv:gr-qc/0512137

    Pith/arXiv arXiv 2006

  17. [17]

    Gair and L

    J. Gair and L. Wen, Classical and Quantum Gravity22, S1359–S1371 (2005)

    2005

  18. [18]

    J. R. Gair, I. Mandel, and L. Wen, Class. Quant. Grav.25, 184031 (2008), arXiv:0804.1084 [gr-qc]

    Pith/arXiv arXiv 2008

  19. [19]

    Speri, R

    L. Speri, R. Tenorio, C. Chapman-Bird, and D. Gerosa, Physical Review D113, 10.1103/dh3j-ksfl (2026)

    doi:10.1103/dh3j-ksfl 2026

  20. [20]

    Bandopadhyay, C

    D. Bandopadhyay, C. E. A. Chapman-Bird, and A. Vecchio, Global time-frequency search for stellar-mass binary black holes in LISA (2026), arXiv:2510.19047 [gr-qc]

    Pith/arXiv arXiv 2026

  21. [21]

    Colpiet al., arXiv e-prints , arXiv:2402.07571 (2024), arXiv:2402.07571 [astro-ph.CO]

    M. Colpiet al., arXiv e-prints , arXiv:2402.07571 (2024), arXiv:2402.07571 [astro-ph.CO]

    Pith/arXiv arXiv 2024

  22. [22]

    Luoet al., Class

    J. Luoet al., Class. Quant. Grav.33, 035010 (2016), arXiv:1512.02076 [astro-ph.IM]

    Pith/arXiv arXiv 2016

  23. [23]

    N. J. Cornish, T. B. Littenberg, B. B´ ecsy, K. Chatziioannou, J. A. Clark, S. Ghonge, and M. Millhouse, Phys. Rev. D 103, 044006 (2021), arXiv:2011.09494 [gr-qc]

    Pith/arXiv arXiv 2021

  24. [24]

    Necula, S

    V. Necula, S. Klimenko, and G. Mitselmakher, Journal of Physics: Conference Series363, 012032 (2012)

    2012

  25. [25]

    M. J. Szczepa´ nczyk, F. Salemi, S. Bini, T. Mishra, G. Vedovato, V. Gayathri, I. Bartos, S. Bhaumik, M. Drago, O. Halim, C. Lazzaro, A. Miani, E. Milotti, G. A. Prodi, S. Tiwari, and S. Klimenko, Phys. Rev. D107, 062002 (2023), arXiv:2210.01754 [gr-qc]

    Pith/arXiv arXiv 2023

  26. [26]

    Drago, S

    M. Drago, S. Klimenko, C. Lazzaro, E. Milotti, G. Mitselmakher, V. Necula, B. O’Brian, G. A. Prodi, F. Salemi, M. Szczepanczyk, S. Tiwari, V. Tiwari, V. Gayathri, G. Vedovato, and I. Yakushin, SoftwareX14, 100678 (2021), arXiv:2006.12604 [gr-qc]. 23

    Pith/arXiv arXiv 2021

  27. [27]

    W. G. Anderson, P. R. Brady, J. D. E. Creighton, and E. E. Flanagan, Phys. Rev. D63, 042003 (2001), arXiv:gr- qc/0008066

    arXiv 2001

  28. [28]

    Klimenko and G

    S. Klimenko and G. Mitselmakher, Class. Quant. Grav.21, S1819 (2004)

    2004

  29. [29]

    M. C. Digman and N. J. Cornish, WDMWaveletTransforms: Fast forward and inverse WDM wavelet transforms, Astro- physics Source Code Library, record ascl:2307.037 (2023), ascl:2307.037

    2023

  30. [30]

    M. C. Digman and N. J. Cornish, Phys. Rev. D108, 023022 (2023), arXiv:2212.04600 [gr-qc]

    Pith/arXiv arXiv 2023

  31. [31]

    M. C. Digman and N. J. Cornish, Astrophys. J.940, 10 (2022), arXiv:2206.14813 [astro-ph.IM]

    Pith/arXiv arXiv 2022

  32. [32]

    Pearson and N

    N. Pearson and N. J. Cornish, Phys. Rev. D113, 064033 (2026), arXiv:2509.05479 [gr-qc]

    Pith/arXiv arXiv 2026

  33. [33]

    N. J. Cornish, Phys. Rev. D102, 124038 (2020), arXiv:2009.00043 [gr-qc]

    Pith/arXiv arXiv 2020

  34. [34]

    N. J. Cornish, arXiv e-prints , arXiv:2511.10632 (2025), arXiv:2511.10632 [gr-qc]

    Pith/arXiv arXiv 2025

  35. [35]

    N. J. Cornish, WDM Transform: Codes to compute the WDM wavelet transform,https://github.com/ eXtremeGravityInstitute/WDM_Transform(2020), eXtreme Gravity Institute, Montana State University. Implements fast TaylorT, TaylorF, and sparse WDM transforms for binary chirp signals in C. Accessed: May 2026

    2020

  36. [36]

    M. C. Digman, WDMWaveletTransforms: Fast forward and inverse WDM wavelet transforms in Python,https:// github.com/XGI-MSU/WDMWaveletTransforms(2022), eXtreme Gravity Institute, Montana State University. Produced under NASA LISA Preparatory Science Grant 80NSSC19K0320. Released v0.0.1, April 2025; GPLv2+ licence. Accessed: May 2026

    2022

  37. [37]

    C. J. Moore, WDM GW wavelets: A fast, JAX-based Python implementation of the Wilson–Daubechies–Meyer wavelet transform for gravitational wave data,https://cjm96.github.io/WDM_GW_wavelets/(2025), gitHub repository:https: //github.com/cjm96/WDM_GW_wavelets. Includes time-delay filter routines for space-based detector responses. Accessed: May 2026

    2025

  38. [38]

    GitLab:https://gitlab.com/gwburst/public/library

    cWB Team, coherent WaveBurst (cWB), version cWB-6.4.6.0,https://gwburst.gitlab.io/documentation/latest/ html/(2024), public release used for the LVK O4b analysis. GitLab:https://gitlab.com/gwburst/public/library. GPLv3 licence. Accessed: May 2026

    2024

  39. [39]

    W. M. Farr, WDMWavelets.jl: WDM wavelets in Julia,https://github.com/farr/WDMWavelets.jl(2025), gitHub repository

    2025

  40. [40]

    Krishnan, A

    B. Krishnan, A. M. Sintes, M. A. Papa, B. F. Schutz, S. Frasca, and C. Palomba, Phys. Rev. D70, 082001 (2004), arXiv:gr-qc/0407001

    Pith/arXiv arXiv 2004

  41. [41]

    Palomba, P

    C. Palomba, P. Astone, and S. Frasca, Classical and Quantum Gravity22, S1255 (2005)

    2005

  42. [42]

    Savalle, J

    E. Savalle, J. Gair, L. Speri, and S. Babak, Phys. Rev. D106, 022003 (2022)

    2022

  43. [43]

    Daubechies, S

    I. Daubechies, S. Jaffard, and J. L. Journe, SIAM J. Math. Anal.22, 554 (1991)

    1991

  44. [44]

    Burke, S

    O. Burke, S. Marsat, J. R. Gair, and M. L. Katz, Phys. Rev. D111, 124053 (2025), arXiv:2502.17426 [gr-qc]

    Pith/arXiv arXiv 2025

  45. [45]

    Bartoloet al., J

    N. Bartoloet al., J. Cosmol. Astropart. Phys.2022, 009 (2022), arXiv:2201.08782 [astro-ph.CO]

    Pith/arXiv arXiv 2022

  46. [46]

    Bayle, M

    J.-B. Bayle, M. Le Jeune, and J. Menu, jaxgb: Fast LISA response for Galactic binaries using JAX,https://pypi.org/ project/jaxgb/(2025), version 0.2.1, BSD-3-Clause

    2025

  47. [47]

    D. Phan, N. Pradhan, and M. Jankowiak, arXiv e-prints , arXiv:1912.11554 (2019), arXiv:1912.11554 [stat.ML]

    Pith/arXiv arXiv 1912

  48. [48]

    Klimenko, S

    S. Klimenko, S. Mohanty, M. Rakhmanov, and G. Mitselmakher, Phys. Rev. D72, 122002 (2005), arXiv:gr-qc/0508068

    Pith/arXiv arXiv 2005

  49. [49]

    Klimenkoet al., Phys

    S. Klimenkoet al., Phys. Rev. D93, 042004 (2016), arXiv:1511.05999 [gr-qc]

    Pith/arXiv arXiv 2016

  50. [50]

    Vajpeyi, G

    A. Vajpeyi, G. Mentasti, Q. Baghi, O. Burke, and L. Speri, pywavelet/wdm transform: v0.05 (2026)

    2026