arxiv: 2605.13738 · v1 · submitted 2026-05-13 · 🌌 astro-ph.IM · astro-ph.CO· gr-qc· hep-ex

Recognition: no theorem link

Inpainting over the cracks: challenges of applying pre-merger searches for massive black hole binaries to realistic LISA datasets

Gareth Cabourn Davies , Ian Harry

Authors on Pith no claims yet

Pith reviewed 2026-05-14 17:38 UTC · model grok-4.3

classification 🌌 astro-ph.IM astro-ph.COgr-qchep-ex

keywords LISAmassive black hole binariespre-merger detectioninpaintinggravitational wave data analysisdata gapssignal recovery

0 comments

The pith

Inpainting recovers 14 pre-merger massive black hole binary signals in realistic LISA data even with gaps.

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

This paper tests two approaches for spotting massive black hole binary mergers in LISA data before they happen. The goal is to give electromagnetic telescopes enough advance warning to observe the event at the same time. They apply a standard zero-latency filter and a new inpainting method to the Sangria-HM challenge dataset. Both recover the 14 loudest signals expected to be visible at least half a day early. The inpainting approach succeeds even when three-day gaps are inserted into the data stream.

Core claim

The authors show that an inpainting technique can identify pre-merger signals from massive black hole binaries in realistic LISA datasets even when data gaps are present. They recover all 14 expected signals from the Sangria-HM dataset at least half a day before merger. The method also allows identification of quieter signals in regions with overlapping mergers by subtracting previously detected ones.

What carries the argument

The inpainting technique, which fills missing data segments to enable continuous searches for gravitational-wave signals.

If this is right

Pre-merger alerts for electromagnetic follow-up become feasible even with realistic data interruptions.
Overlapping signals can be separated by sequential detection and subtraction of louder events.
The zero-latency filter provides a baseline that inpainting extends to gapped data.
Early sky localization estimates can be updated as more data arrives before merger.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Real LISA data pipelines will require robust gap-filling to achieve reliable early-warning science.
The approach could extend to other space-based gravitational-wave observatories facing similar data interruptions.
Further tests on varied overlap configurations would clarify how many signals can be peeled away sequentially.

Load-bearing premise

The Sangria-HM simulated dataset captures the essential noise properties and signal overlaps of real LISA data without introducing artificial biases.

What would settle it

Applying the inpainting search to actual LISA flight data and failing to recover a known massive black hole binary at the predicted pre-merger time would show the method does not generalize.

Figures

Figures reproduced from arXiv: 2605.13738 by Gareth Cabourn Davies, Ian Harry.

**Figure 1.** Figure 1: FIG. 1: Comparison of power spectral densities used in this work. The power spectral density is estimated from the [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗

**Figure 2.** Figure 2: FIG. 2: The method used to smooth over the dip in the [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3: Template bank sizes assuming different [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4: Optimal signal-to-noise ratio (SNR) build-up [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5: Signal 10 results, showing results with and [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6: Search results for the zero-latency filter search [PITH_FULL_IMAGE:figures/full_fig_p006_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7: Signal-to-noise ratio (SNR) as a function of the forecast merger time and the time before merger. Produced [PITH_FULL_IMAGE:figures/full_fig_p008_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8: Timeline of results for Signal Zero, [PITH_FULL_IMAGE:figures/full_fig_p009_8.png] view at source ↗

**Figure 9.** Figure 9: FIG. 9: Signal-to-Noise Ratio (SNR) vs forecast end time results for the inpainting search in the congested region of [PITH_FULL_IMAGE:figures/full_fig_p010_9.png] view at source ↗

**Figure 10.** Figure 10: FIG. 10: Timeline of results in the congested period of days 110-140 of the Sangria-HM dataset. These are [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗

read the original abstract

A key science target of the Large Interferometer Space Antenna (LISA) is to carry out multi-messenger observations of massive black hole binaries, observing the merger simultaneously in gravitational waves and with electromagnetic observatories. Identifying that a merger is happening and providing an updating estimate of the sky location in the hours, days and weeks before the merger is critical to enable electromagnetic observations of the merger event. In this work we demonstrate and compare two methods for premerger identification of massive black hole binaries; a zero-latency filter approach and, for the first time, an approach using an ``inpainting'' technique. We apply these methods to the LISA Data Challenge dataset 2a--Sangria-HM--and demonstrate the successful recovery of the 14 signals in the dataset that we expected to be identifiable at least half a day before merger. We show that the inpainting method can identify premerger signals even when gaps are present in the data, demonstrating the recovery of a signal even when 3 day-long data gaps are added to the 14 days preceding merger. Finally, we explore the challenge of overlapping signals, using a region of overlapping signals in the Sangria-HM dataset, all of which merge within a 10-day window, and show how removing signals that have been confidently identified from the data allows us to identify quieter signals in the same period.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

Inpainting recovers the targeted pre-merger signals on Sangria-HM even after added gaps, but the whole result sits on how faithfully that one simulation matches real LISA noise and gap statistics.

read the letter

The main takeaway is that inpainting lets them pull out the 14 expected massive black hole binary signals from the Sangria-HM dataset at least half a day before merger, and the same signals stay recoverable after they insert 3-day gaps into the final 14 days. They also show that subtracting the louder overlapping signals in a 10-day window lets the quieter ones come through. This is the first time inpainting has been tried for pre-merger LISA searches, and the comparison to the zero-latency filter is straightforward and useful on the public data challenge set. The demonstration is concrete: they report which signals they recover and what happens when gaps are added. That gives people working on early-alert pipelines something they can actually test against the same dataset. The soft spot is that everything rests on Sangria-HM being a good stand-in for real LISA data. The noise model, gap timing, and signal density are fixed in this one injection campaign, so any mismatch in non-stationarity or overlapping-signal statistics would change the recovery rates. No alternative simulations or injected real-data proxies are shown, which keeps the evidence at the level of a solid proof-of-concept rather than a broad validation. This paper is for LISA data analysts who need practical ways to handle gaps in pre-merger searches. A reader building or testing alert pipelines would get direct value from the inpainting example and the overlapping-signal test. It deserves peer review because the method is new in this setting, the dataset is public, and the results are reported clearly enough for referees to check the assumptions and ask for more cross-checks if needed. I would send it to referees.

Referee Report

3 major / 2 minor

Summary. The manuscript presents and compares two methods—a zero-latency filter and an inpainting technique—for identifying massive black hole binary mergers in LISA data prior to merger. Using the Sangria-HM dataset from the LISA Data Challenge, the authors demonstrate recovery of 14 pre-selected signals at least half a day before merger, show that inpainting allows recovery even with 3-day data gaps inserted in the 14 days before merger, and illustrate how removing confidently identified signals aids in detecting overlapping quieter signals within a 10-day window.

Significance. If the results hold, this work provides practical tools for early warning of LISA mergers, enabling multi-messenger observations. The application to a public data challenge dataset and the handling of realistic features like gaps and overlaps are strengths, offering a foundation for operational LISA data analysis pipelines.

major comments (3)

[Results on Sangria-HM dataset] The central empirical support consists of recovery demonstrations on a single controlled injection campaign (Sangria-HM) with fixed noise model and gap statistics. No cross-checks against alternative LISA simulations or injected real-data proxies are described, which is load-bearing for the claim that the methods address 'realistic LISA challenges'.
[Overlapping signals analysis] The section exploring overlapping signals shows qualitative improvement by removing identified signals, but lacks quantitative metrics such as detection thresholds, false alarm rates, or SNR improvements to substantiate the claim of identifying quieter signals.
[Methods and results] Detailed error analysis, validation metrics (e.g., recovery time distributions, false positive rates), and comparison baselines to existing search methods are not visible in the reported results, weakening the assessment of the methods' performance.

minor comments (2)

[Abstract] The abstract states 'successful recovery' without specifying the criteria or metrics used to define success.
[Notation and terminology] Ensure consistent use of terms like 'inpainting' and 'zero-latency filter' throughout, with clear definitions in the methods section.

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for their thoughtful review and constructive feedback on our manuscript. We address each of the major comments below, providing clarifications and outlining revisions to improve the paper.

read point-by-point responses

Referee: [Results on Sangria-HM dataset] The central empirical support consists of recovery demonstrations on a single controlled injection campaign (Sangria-HM) with fixed noise model and gap statistics. No cross-checks against alternative LISA simulations or injected real-data proxies are described, which is load-bearing for the claim that the methods address 'realistic LISA challenges'.

Authors: The Sangria-HM dataset is specifically constructed to simulate realistic LISA data challenges, including instrumental noise and data gaps, as part of the LISA Data Challenge. While we agree that additional validation on other datasets would be valuable, this work focuses on demonstrating the methods on this public benchmark. In the revised manuscript, we will expand the discussion section to explicitly address the scope of our claims and the need for future cross-validation with other simulations. revision: partial
Referee: [Overlapping signals analysis] The section exploring overlapping signals shows qualitative improvement by removing identified signals, but lacks quantitative metrics such as detection thresholds, false alarm rates, or SNR improvements to substantiate the claim of identifying quieter signals.

Authors: We acknowledge that the analysis of overlapping signals is presented qualitatively. To address this, we will add quantitative metrics in the revised version, including estimated SNR improvements and approximate false alarm rate reductions for the quieter signals after subtracting the confidently identified ones. revision: yes
Referee: [Methods and results] Detailed error analysis, validation metrics (e.g., recovery time distributions, false positive rates), and comparison baselines to existing search methods are not visible in the reported results, weakening the assessment of the methods' performance.

Authors: Some validation is included, such as the recovery of all 14 expected signals at least half a day before merger. However, we agree that more detailed metrics would strengthen the presentation. We will revise the results section to include recovery time distributions, estimates of false positive rates derived from the dataset, and comparisons to standard search techniques like matched filtering. revision: yes

Circularity Check

0 steps flagged

No significant circularity; empirical demonstration on external public dataset

full rationale

The paper applies zero-latency filtering and inpainting to recover 14 pre-selected signals from the external LISA Data Challenge Sangria-HM dataset (including a test with added 3-day gaps). Recovery is measured directly against the known injected signals in that fixed public benchmark. No derivation chain exists that reduces by construction to fitted parameters, self-definitions, or self-citations; the reported success rates are falsifiable against the external injections. The methods are therefore self-contained against an independent data challenge.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The claims depend on the fidelity of the simulated dataset and standard assumptions in gravitational wave data analysis.

free parameters (1)

Detection thresholds
Likely used to identify signals but not specified in abstract.

axioms (1)

domain assumption The LISA Data Challenge dataset 2a Sangria-HM accurately represents realistic LISA observations including noise and signal characteristics.
All results are based on performance on this simulated dataset.

pith-pipeline@v0.9.0 · 5561 in / 1163 out tokens · 60113 ms · 2026-05-14T17:38:20.480333+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

40 extracted references · 33 canonical work pages · 17 internal anchors

[1]

Laser Interferometer Space Antenna

P. Amaro-Seoaneet al.(LISA), Laser Interferometer Space Antenna (2017), arXiv:1702.00786 [astro-ph.IM]

work page internal anchor Pith review Pith/arXiv arXiv 2017
[2]

LISA Definition Study Report

M. Colpiet al., LISA Definition Study Report (2024), arXiv:2402.07571 [astro-ph.CO]

work page internal anchor Pith review Pith/arXiv arXiv 2024
[3]

P. A. Seoaneet al.(LISA), Astrophysics with the Laser Interferometer Space Antenna, Living Rev. Rel.26, 2 (2023), arXiv:2203.06016 [gr-qc]

work page arXiv 2023
[4]

K. A. Arnaudet al., The Mock LISA Data Challenges: An overview, AIP Conf. Proc.873, 619 (2006), arXiv:gr- qc/0609105

work page arXiv 2006
[5]

K. A. Arnaudet al.(Mock LISA Data Challenge Task Force), A How-To for the Mock LISA Data Challenges, 12 AIP Conf. Proc.873, 625 (2006), arXiv:gr-qc/0609106

work page internal anchor Pith review Pith/arXiv arXiv 2006
[6]

K. A. Arnaudet al., Report on the first round of the Mock LISA data challenges, Class. Quant. Grav.24, S529 (2007), arXiv:gr-qc/0701139

work page internal anchor Pith review Pith/arXiv arXiv 2007
[7]

Report on the second Mock LISA Data Challenge

S. Babaket al.(Mock LISA Data Challenge Task Force), Report on the second Mock LISA Data Challenge, Class. Quant. Grav.25, 114037 (2008), arXiv:0711.2667 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[8]

The Mock LISA Data Challenges: from Challenge 1B to Challenge 3

S. Babaket al., The Mock LISA Data Challenges: From Challenge 1B to Challenge 3, Class. Quant. Grav.25, 184026 (2008), arXiv:0806.2110 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[9]

K. A. Arnaudet al., An Overview of the second round of the Mock LISA Data Challenges, Class. Quant. Grav. 24, S551 (2007), arXiv:gr-qc/0701170

work page internal anchor Pith review Pith/arXiv arXiv 2007
[10]

The Mock LISA Data Challenges: from Challenge 3 to Challenge 4

S. Babaket al.(Mock LISA Data Challenge Task Force), The Mock LISA Data Challenges: From Challenge 3 to Challenge 4, Class. Quant. Grav.27, 084009 (2010), arXiv:0912.0548 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2010
[11]

Baghi (LDC Working Group), The LISA Data Chal- lenges, in56th Rencontres de Moriond on Gravitation (2022) arXiv:2204.12142 [gr-qc]

Q. Baghi (LDC Working Group), The LISA Data Chal- lenges, in56th Rencontres de Moriond on Gravitation (2022) arXiv:2204.12142 [gr-qc]

work page arXiv 2022
[12]

lal.in2p3.fr/()

The New LISA Data Challenges,https://lisa-ldc. lal.in2p3.fr/()
[13]

Cabourn Davieset al., Premerger observation and characterization of massive black hole binaries, Phys

G. Cabourn Davieset al., Premerger observation and characterization of massive black hole binaries, Phys. Rev. D111, 043045 (2025), arXiv:2411.07020 [hep-ex]

work page arXiv 2025
[14]

Houba, S

N. Houba, S. H. Strub, L. Ferraioli, and D. Giardini, Detection and prediction of future massive black hole mergers with machine learning and truncated waveforms, Phys. Rev. D110, 062003 (2024), arXiv:2405.11340 [astro-ph.IM]

work page arXiv 2024
[15]

Ruan and Z.-K

W.-H. Ruan and Z.-K. Guo, Premerger detection of mas- sive black hole binaries using deep learning, Phys. Rev. D109, 123031 (2024), arXiv:2402.16282 [astro-ph.IM]

work page arXiv 2024
[16]

Pre-Merger Localization of Gravitational-Wave Standard Sirens With LISA: Triggered Search for an Electromagnetic Counterpart

B. Kocsis, Z. Haiman, and K. Menou, Pre-Merger Lo- calization of Gravitational-Wave Standard Sirens With LISA: Triggered Search for an Electromagnetic Coun- terpart, Astrophys. J.684, 870 (2008), arXiv:0712.1144 [astro-ph]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[17]

S. T. McWilliams, R. N. Lang, J. G. Baker, and J. I. Thorpe, Sky localization of complete inspiral-merger- ringdown signals for nonspinning massive black hole bi- naries, Phys. Rev. D84, 064003 (2011), arXiv:1104.5650 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2011
[18]

Saini, S

P. Saini, S. A. Bhat, and K. G. Arun, Premerger localization of intermediate mass binary black holes with LISA and prospects of joint observations with Athena and LSST, Phys. Rev. D106, 104015 (2022), arXiv:2208.03004 [gr-qc]

work page arXiv 2022
[19]

Dal Canton, A

T. Dal Canton, A. Mangiagli, S. C. Noble, J. Schnittman, A. Ptak, A. Klein, A. Sesana, and J. Camp, De- tectability of modulated X-rays from LISA’s supermas- sive black hole mergers, Astrophys. J.886, 146 (2019), arXiv:1902.01538 [astro-ph.HE]

work page arXiv 2019
[20]

Castelli, Q

E. Castelli, Q. Baghi, J. G. Baker, J. Slutsky, J. Bobin, N. Karnesis, A. Petiteau, O. Sauter, P. Wass, and W. J. Weber, Extracting gravitational wave signals from LISA data in the presence of artifacts, Class. Quant. Grav.42, 065018 (2025), arXiv:2411.13402 [gr-qc]

work page arXiv 2025
[21]

Zackay, T

B. Zackay, T. Venumadhav, J. Roulet, L. Dai, and M. Zaldarriaga, Detecting gravitational waves in data with non-stationary and non-Gaussian noise, Phys. Rev. D104, 063034 (2021), arXiv:1908.05644 [astro-ph.IM]

work page arXiv 2021
[22]

in2p3.fr/challenge2a()

The New LISA Data Challenges,https://lisa-ldc. in2p3.fr/challenge2a()
[23]

García-Quirós, M

C. Garc´ ıa-Quir´ os, M. Colleoni, S. Husa, H. Estell´ es, G. Pratten, A. Ramos-Buades, M. Mateu-Lucena, and R. Jaume, Multimode frequency-domain model for the gravitational wave signal from nonprecessing black-hole binaries, Phys. Rev. D102, 064002 (2020), arXiv:2001.10914 [gr-qc]

work page arXiv 2020
[24]

Annual modulation of the Galactic binary confusion noise bakground and LISA data analysis

N. Seto, Annual modulation of the galactic binary con- fusion noise background and LISA data analysis, Phys. Rev. D69, 123005 (2004), arXiv:gr-qc/0403014

work page internal anchor Pith review Pith/arXiv arXiv 2004
[25]

A Solution to the Galactic Foreground Problem for LISA

J. Crowder and N. Cornish, A Solution to the Galactic Foreground Problem for LISA, Phys. Rev. D75, 043008 (2007), arXiv:astro-ph/0611546

work page internal anchor Pith review Pith/arXiv arXiv 2007
[26]

S. H. Strub, L. Ferraioli, C. Schmelzbach, S. C. St¨ ahler, and D. Giardini, Global analysis of LISA data with Galactic binaries and massive black hole binaries, Phys. Rev. D110, 024005 (2024), arXiv:2403.15318 [gr-qc]

work page arXiv 2024
[27]

C. J. Moore, R. H. Cole, and C. P. L. Berry, Gravitational-wave sensitivity curves, Class. Quant. Grav.32, 015014 (2015), arXiv:1408.0740 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2015
[28]

Marsat, J

S. Marsat, J. G. Baker, and T. Dal Canton, Explor- ing the Bayesian parameter estimation of binary black holes with LISA, Phys. Rev. D103, 083011 (2021), arXiv:2003.00357 [gr-qc]

work page arXiv 2021
[29]

Building a stochastic template bank for detecting massive black hole binaries

S. Babak, Building a stochastic template bank for detect- ing massive black hole binaries, Class. Quant. Grav.25, 195011 (2008), arXiv:0801.4070 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2008
[30]

I. W. Harry, B. Allen, and B. S. Sathyaprakash, A Stochastic template placement algorithm for gravita- tional wave data analysis, Phys. Rev. D80, 104014 (2009), arXiv:0908.2090 [gr-qc]

work page internal anchor Pith review Pith/arXiv arXiv 2009
[31]

A chi-squared time-frequency discriminator for gravitational wave detection

B. Allen,χ 2 time-frequency discriminator for gravita- tional wave detection, Phys. Rev. D71, 062001 (2005), arXiv:gr-qc/0405045

work page internal anchor Pith review Pith/arXiv arXiv 2005
[32]

Application of a Zero-latency Whitening Filter to Compact Binary Coalescence Gravitational-wave Searches

L. Tsukada, K. Cannon, C. Hanna, D. Keppel, D. Meacher, and C. Messick, Application of a Zero- latency Whitening Filter to Compact Binary Coalescence Gravitational-wave Searches, Phys. Rev. D97, 103009 (2018), arXiv:1708.04125 [astro-ph.IM]

work page internal anchor Pith review Pith/arXiv arXiv 2018
[33]

C. D. Capano, M. Cabero, J. Westerweck, J. Abedi, S. Kastha, A. H. Nitz, Y.-F. Wang, A. B. Nielsen, and B. Krishnan, Multimode Quasinormal Spectrum from a Perturbed Black Hole, Phys. Rev. Lett.131, 221402 (2023), arXiv:2105.05238 [gr-qc]

work page arXiv 2023
[34]

Allen, W

B. Allen, W. G. Anderson, P. R. Brady, D. A. Brown, and J. D. E. Creighton, FINDCHIRP: An Algorithm for detection of gravitational waves from inspiraling com- pact binaries, Phys. Rev. D85, 122006 (2012), arXiv:gr- qc/0509116

work page arXiv 2012
[35]

Katz and J

M. Katz and J. Roberts, mikekatz04/bbhx: New release! (2023)

2023
[36]

Cannon, S

K. Cannon, S. Caudill, C. Chan, B. Cousins, J. D. Creighton, B. Ewing, H. Fong, P. Godwin, C. Hanna, S. Hooper, R. Huxford, R. Magee, D. Meacher, C. Mes- sick, S. Morisaki, D. Mukherjee, H. Ohta, A. Pace, S. Privitera, I. de Ruiter, S. Sachdev, L. Singer, D. Singh, R. Tapia, L. Tsukada, D. Tsuna, T. Tsutsui, K. Ueno, A. Viets, L. Wade, and M. Wade, Gstlal...

work page arXiv 2021
[37]

J. D. Hunter, Matplotlib: A 2d graphics environment, Computing in Science & Engineering9, 90 (2007)

2007
[38]

C. R. Harris, K. J. Millman, S. J. van der Walt, R. Gom- mers, P. Virtanen, D. Cournapeau, E. Wieser, J. Tay- lor, S. Berg, N. J. Smith, R. Kern, M. Picus, S. Hoyer, 13 M. H. van Kerkwijk, M. Brett, A. Haldane, J. F. del R´ ıo, M. Wiebe, P. Peterson, P. G´ erard-Marchant, K. Shep- pard, T. Reddy, W. Weckesser, H. Abbasi, C. Gohlke, and T. E. Oliphant, Arr...

2020
[39]

A. Nitz, I. Harry, D. Brown, C. M. Biwer, J. Willis, T. D. Canton, C. Capano, T. Dent, L. Pekowsky, G. S. C. Davies, S. De, M. Cabero, S. Wu, A. R. Williamson, B. Machenschalk, D. Macleod, F. Pannarale, P. Kumar, S. Reyes, dfinstad, S. Kumar, M. T´ apai, L. Singer, P. Ku- mar, veronica villa, maxtrevor, B. U. V. Gadre, S. Khan, S. Fairhurst, and A. Tolley...

2024
[40]

Virtanen, R

P. Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P. Pe- terson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wilson, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙I. Po- lat, Y. Feng, E. W. Moore, J. VanderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. Henr...

2020