pith. sign in

arxiv: 2601.07704 · v2 · submitted 2026-01-12 · 🧮 math.NA · cs.NA

TMATDG: applying TDG methods to multiple scattering via T-matrix approximation

Armando Maria Monforte This is my paper

Pith reviewed 2026-05-16 14:53 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords multiple scatteringTrefftz methodsT-matrixDiscontinuous GalerkinHelmholtz equationMATLAB softwarepolygonal scatterers
0
0 comments X

The pith

A MATLAB package couples Trefftz Discontinuous Galerkin methods with T-matrix approximations to solve multiple scattering problems.

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

The paper introduces a MATLAB package that integrates Trefftz Discontinuous Galerkin methods, which solve local Helmholtz scattering problems, with the T-matrix method for handling interactions between multiple scatterers. This combination allows for the treatment of scattering by several polygonal obstacles using numerically computed T-matrices from the TMATROM package. The approach avoids the need to discretize the entire computational domain for each interaction. Readers interested in computational methods for wave propagation would value this as a practical tool for simulating acoustic or electromagnetic scattering in complex geometries.

Core claim

The TMATDG package provides a framework to solve multiple scattering problems by coupling Trefftz Discontinuous Galerkin methods for Helmholtz scattering with the T-matrix method, relying on numerical approximations of T-matrices to handle interactions among polygonal obstacles.

What carries the argument

The coupling of Trefftz Discontinuous Galerkin (TDG) methods with T-matrix approximations, which allows local solutions around each obstacle to be combined globally via scattering matrices.

If this is right

  • Multiple scattering interactions for polygonal obstacles can be computed without full-domain meshing.
  • The method extends TDG applicability from single to multiple scatterers.
  • Reliance on precomputed T-matrices improves efficiency for repeated or parametric studies.

Where Pith is reading between the lines

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

  • This coupling could be adapted for time-domain problems or nonlinear scattering.
  • Similar hybrid approaches might apply to other boundary integral or finite element methods for wave problems.
  • Validation on standard benchmarks like two spheres or cylinders would confirm stability.

Load-bearing premise

The numerical T-matrices from TMATROM can be accurately coupled with TDG without introducing significant errors or instabilities in the multiple scattering solution.

What would settle it

A numerical test showing large discrepancies between the computed far-field patterns and known analytical solutions for a simple two-obstacle configuration would falsify the claim.

Figures

Figures reproduced from arXiv: 2601.07704 by Armando Maria Monforte.

Figure 1
Figure 1. Figure 1: Domain geometry for Dirichlet (left) and transmission (right) problem. [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Real part of total field for the square penetrable scatterer of Section 4.3.1 in 3 [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Real part of total field for the ensemble of scatterers of Section 4.3.2 with two [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Left: the obstacle arrangement of Section 4.3.3. Right: the [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Real part of total field for the scatterer arrangement of Section 4.3.3, with [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
read the original abstract

We present a MATLAB package for the solution of multiple scattering problems, coupling Trefftz Discontinuos Galerkin methods for Helmholtz scattering with the T-matrix method. We rely on the TMATROM package to numerically approximate the T-matrices and deal with multiple scattering problem, providing a framework to handle scattering by polygonal obstacles.

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.

Referee Report

2 major / 2 minor

Summary. The manuscript presents a MATLAB package TMATDG that couples Trefftz Discontinuous Galerkin (TDG) methods for the Helmholtz equation with the T-matrix method for multiple scattering, relying on the external TMATROM package to numerically compute T-matrix approximations for polygonal obstacles.

Significance. If the coupling is shown to be accurate and stable, the package would supply a practical computational framework for multiple-scattering problems involving non-smooth polygonal scatterers, combining local TDG resolution with far-field T-matrix interactions. The work would be of interest to the numerical acoustics and electromagnetics communities provided it includes reproducible validation.

major comments (2)
  1. [Abstract and §2] Abstract and §2: the central claim that the TDG-T-matrix coupling via TMATROM accurately handles multiple scattering by polygons rests on the unverified assumption that numerically computed T-matrices remain well-conditioned and truncation-error controlled when inserted into the TDG formulation; no error bounds, condition-number estimates, or comparison against reference solutions (boundary-integral or Mie-type) are supplied for even two interacting polygons.
  2. [§3] §3 (method description): the coupling procedure is described at a high level without explicit equations showing how the T-matrix far-field operator is inserted into the TDG variational form or how the polygonal corner singularities are accommodated in the T-matrix truncation; this omission prevents verification that the scheme does not introduce spurious artifacts at the interfaces.
minor comments (2)
  1. [Abstract] Abstract: 'Trefftz Discontinuos Galerkin' contains a typographical error and should read 'Trefftz Discontinuous Galerkin'.
  2. [Numerical results] The manuscript should include at least one table or figure reporting L2 or far-field errors versus truncation order for a single polygon and for a two-polygon configuration to substantiate the claims.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and the constructive comments on our manuscript describing the TMATDG MATLAB package. We address each major comment below and indicate the changes we will make in the revised version.

read point-by-point responses

  1. Referee: [Abstract and §2] Abstract and §2: the central claim that the TDG-T-matrix coupling via TMATROM accurately handles multiple scattering by polygons rests on the unverified assumption that numerically computed T-matrices remain well-conditioned and truncation-error controlled when inserted into the TDG formulation; no error bounds, condition-number estimates, or comparison against reference solutions (boundary-integral or Mie-type) are supplied for even two interacting polygons.

    Authors: We acknowledge that the manuscript does not supply explicit error bounds, condition-number estimates, or direct comparisons against reference solutions for the coupled scheme. The current focus is on the software framework and its practical use for polygonal scatterers. In the revised version we will add a dedicated numerical validation subsection that includes (i) comparisons of the TMATDG results against a boundary-integral reference solver for two interacting polygons and (ii) plots of the condition numbers of the T-matrices obtained from TMATROM for increasing truncation orders. These additions will provide concrete evidence of stability and accuracy for the cases considered. revision: yes

  2. Referee: [§3] §3 (method description): the coupling procedure is described at a high level without explicit equations showing how the T-matrix far-field operator is inserted into the TDG variational form or how the polygonal corner singularities are accommodated in the T-matrix truncation; this omission prevents verification that the scheme does not introduce spurious artifacts at the interfaces.

    Authors: We agree that the description in §3 is too high-level. In the revision we will insert the explicit weak-form equations that show how the far-field T-matrix operator enters the TDG variational formulation. We will also add a short paragraph explaining the truncation strategy for the T-matrix and its compatibility with the corner singularities of the polygons, referencing the convergence theory already established in the TMATROM package. These additions should allow readers to verify that the coupling does not introduce spurious interface artifacts. revision: yes

Circularity Check

0 steps flagged

Minor external dependency on TMATROM without internal circular reduction

full rationale

The paper presents a coupling framework that delegates T-matrix computation to the external TMATROM package rather than deriving or fitting those matrices from the TDG formulation itself. No equations or steps within the manuscript reduce a claimed prediction to a fitted input or self-definition by construction, and the central contribution (the MATLAB package for polygonal multiple scattering) retains independent content as a software interface. Any self-citation to TMATROM is not load-bearing for a mathematical derivation but simply identifies the source of the numerical T-matrices, yielding only a low circularity score.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on standard assumptions from scattering theory and numerical PDE methods; no new entities or fitted parameters are mentioned in the abstract.

axioms (1)
  • domain assumption T-matrix approximations can be accurately computed via TMATROM and coupled to TDG for multiple scattering
    Invoked in the description of the package framework for polygonal obstacles.

pith-pipeline@v0.9.0 · 5334 in / 1117 out tokens · 45483 ms · 2026-05-16T14:53:11.045102+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

  1. [1]

    Unified analysis of discontinuous Galerkin methods for elliptic problems

    D. N. Arnold et al. “Unified analysis of discontinuous Galerkin methods for elliptic problems”. In:SIAM J. Numer. Anal.39.5 (2002), pp. 1749–1779

    work page 2002

  2. [2]

    Colton and R

    D. Colton and R. Kress.Inverse acoustic and electromagnetic scattering theory. English. 4th expanded edition. Vol. 93. Appl. Math. Sci. Cham: Springer, 2019

    work page 2019

  3. [3]

    Near- and far-field light scattering by nonspherical particles: Applicability of methods that involve a spherical basis

    V. G. Farafonov, V. B. Il’in, and A. A. Vinokurov. “Near- and far-field light scattering by nonspherical particles: Applicability of methods that involve a spherical basis”. In: Optics and Spectroscopy109.3 (2010), pp. 432–443

    work page 2010

  4. [4]

    A far-field based T-matrix method for three di- mensional acoustic scattering

    M. Ganesh and S. C. Hawkins. “A far-field based T-matrix method for three di- mensional acoustic scattering”. In:Proceedings of the 14th Biennial Computational Techniques and Applications Conference, CTAC-2008. Ed. by G. N. Mercer and A. J. Roberts. Vol. 50. ANZIAM J. Oct. 2008, pp. C121–C136

    work page 2008

  5. [5]

    A far-field based T-matrix method for two dimensional obstacle scattering

    M. Ganesh and S. Hawkins. “A far-field based T-matrix method for two dimensional obstacle scattering”. English. In:ANZIAM Journal51 (2009), pp. C215–C230. 11

    work page 2009

  6. [6]

    Algorithm 975: TMATROM—A T-Matrix Reduced Order Model Software

    M. Ganesh and S. C. Hawkins. “Algorithm 975: TMATROM—A T-Matrix Reduced Order Model Software”. In:ACM Trans. Math. Softw.44.1 (July 2017)

    work page 2017

  7. [7]

    A stochastic pseudospectral and T-matrix algorithm for acoustic scattering by a class of multiple particle configurations

    M. Ganesh and S. Hawkins. “A stochastic pseudospectral and T-matrix algorithm for acoustic scattering by a class of multiple particle configurations”. In:Journal of Quantitative Spectroscopy and Radiative Transfer123 (2013). Peter C. Waterman and his scientific legacy, pp. 41–52

    work page 2013

  8. [8]

    Convergence analysis with parameter esti- mates for a reduced basis acoustic scattering T-matrix method

    M. Ganesh, S. Hawkins, and R. Hiptmair. “Convergence analysis with parameter esti- mates for a reduced basis acoustic scattering T-matrix method”. In:IMA Journal of Numerical Analysis32 (Oct. 2012)

    work page 2012

  9. [9]

    Metamaterial applications of Tmatsolver, an easy-to-use software for simulating multiple wave scattering in two dimensions

    S. C. Hawkins et al. “Metamaterial applications of Tmatsolver, an easy-to-use software for simulating multiple wave scattering in two dimensions”. In:Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences480.2292 (June 2024), p. 20230934

    work page 2024

  10. [10]

    Improving the numerical stability of T-matrix light scattering calculations for extreme particle shapes using the nullfield method with discrete sources

    J. Hellmers, V. Schmidt, and T. Wriedt. “Improving the numerical stability of T-matrix light scattering calculations for extreme particle shapes using the nullfield method with discrete sources”. eng. In:Journal of quantitative spectroscopy and radiative transfer 112.11 (2011), pp. 1679–1686

    work page 2011

  11. [11]

    Homogenized boundary conditions and resonance effects in Faraday cages

    D. P. Hewett and I. J. Hewitt. “Homogenized boundary conditions and resonance effects in Faraday cages”. In:Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences472.2189 (May 2016), p. 20160062

    work page 2016

  12. [12]

    Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of thep-version

    R. Hiptmair, A. Moiola, and I. Perugia. “Plane wave discontinuous Galerkin methods for the 2D Helmholtz equation: analysis of thep-version”. In:SIAM J. Numer. Anal. 49.1 (2011), pp. 264–284

    work page 2011

  13. [13]

    New Generation Finite Element Methods For Forward Seismic Mod- elling

    C. J. Howarth. “New Generation Finite Element Methods For Forward Seismic Mod- elling”. Available onhttps://www.reading.ac.uk/maths-and-stats/publications/ theses- and- dissertations/mathematics- phd- theses. PhD thesis. University of Reading, 2014

    work page 2014

  14. [14]

    A plane wave discontinuous Galerkin method with a Dirichlet- to-Neumann boundary condition for the scattering problem in acoustics

    S. Kapita and P. Monk. “A plane wave discontinuous Galerkin method with a Dirichlet- to-Neumann boundary condition for the scattering problem in acoustics”. In:J. Com- put. Appl. Math.327 (2018), pp. 208–225

    work page 2018

  15. [15]

    Martin.Multiple Scattering

    P. Martin.Multiple Scattering. Interaction of Time-Harmonic Waves with N Obstacles. Vol. 107. Cambridge University Press, 2006

    work page 2006

  16. [16]

    T-matrix computations of light scattering by large spheroidal particles

    M. I. Mishchenko and L. D. Travis. “T-matrix computations of light scattering by large spheroidal particles”. In:Optics Communications109.1 (1994), pp. 16–21

    work page 1994

  17. [17]

    T-matrix computations of light scattering by nonspherical particles: A review

    M. I. Mishchenko, L. D. Travis, and D. W. Mackowski. “T-matrix computations of light scattering by nonspherical particles: A review”. In:Journal of Quantitative Spec- troscopy and Radiative Transfer55.5 (1996). Light Scattering by Non-Spherical Parti- cles, pp. 535–575

    work page 1996

  18. [18]

    Trefftz Discontinuous Galerkin methods for scattering by periodic structures

    A. Moiola and A. M. Monforte. “Trefftz Discontinuous Galerkin methods for scattering by periodic structures”. In:arXiv: 2505.23216(2025)

    work page arXiv 2025

  19. [19]

    New Formulation of Acoustic Scattering

    P. C. Waterman. “New Formulation of Acoustic Scattering”. In:The Journal of the Acoustical Society of America45.6 (June 1969), pp. 1417–1429

    work page 1969

  20. [20]

    Symmetry, Unitarity, and Geometry in Electromagnetic Scattering

    P. C. Waterman. “Symmetry, Unitarity, and Geometry in Electromagnetic Scattering”. In:Phys. Rev. D3 (4 Feb. 1971), pp. 825–839

    work page 1971

  21. [21]

    Matrix formulation of electromagnetic scattering

    P. Waterman. “Matrix formulation of electromagnetic scattering”. In:Proceedings of the IEEE53.8 (1965), pp. 805–812

    work page 1965

  22. [22]

    Improved Mie scattering algorithms

    W. J. Wiscombe. “Improved Mie scattering algorithms”. In:Appl. Opt.19.9 (May 1980), pp. 1505–1509. 12

    work page 1980