Pith. sign in

REVIEW 2 major objections 1 minor 1 cited by

B-spline interpolation on equidistant knots with knot-removal stabilization enables robust high-order matrix compression for the electric field integral equation.

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.3

2026-06-27 05:45 UTC pith:Z7D47MLY

load-bearing objection The knot-removal stabilization lets B-splines replace Lagrange in multilevel EFIE without Runge blowup, but the accuracy preservation after stabilization is asserted more than shown. the 2 major comments →

arxiv 2606.13482 v1 pith:Z7D47MLY submitted 2026-06-11 math.NA cs.NA

A Stabilized Multilevel B-Spline-Based Fast Integral Method for the Solution of the Electric Field Integral Equation

Danijel Juki\'c , Bernd Hofmann , Thomas F. Eibert , Simon B. Adrian This is my paper
classification math.NA cs.NA
keywords B-spline interpolationfast integral methodelectric field integral equationmatrix compressionmultilevel methodRunge phenomenonnull space preservationO(N) complexity
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.

The paper develops a multilevel fast integral method for the electric field integral equation that replaces Lagrange interpolation with B-splines on equidistant knot vectors. This change is intended to eliminate Runge-type instabilities that arise at high interpolation orders in existing FFT-based approaches. The authors introduce a knot-removal stabilization strategy together with exact interlevel transfers via knot insertion to keep the interpolation accurate and well-conditioned. They also supply a factorization that preserves the null space of the scalar potential operator up to machine precision. The resulting scheme is shown to deliver O(N) complexity while remaining compatible with low-frequency preconditioners.

Core claim

Replacing Lagrange interpolation by B-spline interpolation on equidistant knot vectors, stabilized by knot removal and combined with exact interlevel transfers based on knot insertion, yields accurate, well-conditioned multilevel interpolation for the EFIE that preserves the null space of the scalar potential operator up to machine precision and achieves O(N) complexity.

What carries the argument

The knot-removal stabilization strategy for B-spline coefficients on equidistant knot vectors, which prevents ill-conditioning while enabling robust high-order interpolation and exact multilevel transfers.

Load-bearing premise

The knot-removal stabilization strategy combined with exact interlevel transfers based on knot insertion yields accurate, well-conditioned multilevel interpolation without introducing significant new errors or conditioning issues.

What would settle it

Numerical tests on realistic geometries at successively higher interpolation orders that either reproduce the Runge breakdown for Lagrange methods or show loss of accuracy or superlinear scaling once the stabilized B-spline scheme is applied.

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

If this is right

  • Robust high-order interpolation is obtained for both canonical and realistic geometries without the breakdown observed in Lagrange-based methods.
  • The scheme maintains O(N) complexity while preserving the null space of the scalar potential operator up to machine precision.
  • The factorization remains compatible with existing low-frequency preconditioning techniques.
  • Accurate matrix compression is achieved through the stabilized multilevel B-spline interpolation.

Where Pith is reading between the lines

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

  • The same stabilization and transfer technique could be applied to other oscillatory or singular kernels that currently rely on Lagrange interpolation in fast multipole or FFT-accelerated solvers.
  • Higher stable interpolation orders may reduce the number of degrees of freedom required to reach a target accuracy in scattering or antenna problems.
  • The exact knot-insertion transfers could support adaptive refinement strategies without reintroducing conditioning problems at coarse-to-fine interfaces.

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 / 1 minor

Summary. The manuscript presents a multilevel B-spline-based fast integral method for the EFIE that replaces Lagrange interpolation with B-splines on equidistant knot vectors to avoid Runge instabilities. It introduces a knot-removal stabilization strategy together with exact interlevel transfers via knot insertion to produce well-conditioned multilevel interpolation, proposes a factorization that preserves the null space of the scalar potential operator up to machine precision, and reports O(N) complexity together with robust high-order accuracy on both canonical and realistic geometries.

Significance. If the stabilization preserves polynomial reproduction degree and the claimed accuracy, the method would supply a practical high-order alternative to existing FFT-accelerated Lagrange schemes for EFIE solvers, with the added benefit of a null-space-preserving factorization compatible with low-frequency preconditioners. The numerical demonstration on realistic geometries is a concrete strength.

major comments (2)
  1. [stabilization strategy description] The central claim that knot-removal stabilization plus exact knot-insertion transfers yields high-order accurate multilevel interpolation rests on the assertion that the stabilized operator still reproduces polynomials up to the original degree p. No a-priori error bound, Strang-type lemma, or reproduction-property proof is supplied for the stabilized sampling operator; without this, the transition from numerical stability to retained approximation order remains unsecured (see the description of the stabilization strategy and the complexity/accuracy claims).
  2. [numerical-results section] Numerical results are stated to confirm robust high-order behavior and O(N) complexity, yet the manuscript provides neither tabulated error-versus-order data nor a direct comparison of effective polynomial degree before and after knot removal. This omission makes it impossible to verify that the stabilization does not induce order reduction on the EFIE matrix compression (see the numerical-results section).
minor comments (1)
  1. The abstract and introduction would benefit from an explicit statement of the spline degree p employed in the experiments and the precise definition of the sampling matrices whose conditioning is improved by knot removal.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive evaluation of the method's significance and for the constructive major comments. We address each point below and will revise the manuscript to strengthen the presentation where the comments identify gaps.

read point-by-point responses

  1. Referee: The central claim that knot-removal stabilization plus exact knot-insertion transfers yields high-order accurate multilevel interpolation rests on the assertion that the stabilized operator still reproduces polynomials up to the original degree p. No a-priori error bound, Strang-type lemma, or reproduction-property proof is supplied for the stabilized sampling operator; without this, the transition from numerical stability to retained approximation order remains unsecured (see the description of the stabilization strategy and the complexity/accuracy claims).

    Authors: The referee correctly notes the absence of a formal a-priori analysis. The stabilization is constructed so that knot removal is performed only on knots that do not affect the local polynomial reproduction property of the underlying B-spline space, while the exact knot-insertion transfers guarantee that the multilevel operators remain consistent with the original spline space. This design choice is intended to retain the full approximation order p. We will add a concise remark in the revised stabilization-strategy section that explicitly recalls the relevant B-spline reproduction properties and explains why the chosen knot-removal criterion preserves them, thereby making the connection between stability and retained order more transparent. A complete Strang-type lemma or new error bound is beyond the scope of the present work but could be pursued in a follow-up study. revision: partial

  2. Referee: Numerical results are stated to confirm robust high-order behavior and O(N) complexity, yet the manuscript provides neither tabulated error-versus-order data nor a direct comparison of effective polynomial degree before and after knot removal. This omission makes it impossible to verify that the stabilization does not induce order reduction on the EFIE matrix compression (see the numerical-results section).

    Authors: We agree that explicit tabulated verification would strengthen the numerical section. The existing figures already demonstrate that high-order accuracy is maintained on both canonical and realistic geometries without the Runge-type breakdown seen in Lagrange schemes, and that the observed convergence rates are consistent with the nominal spline degree. To make this verification direct, we will add a table in the revised numerical-results section that reports L2 or relative errors versus interpolation order for representative test cases, together with a short comparison of effective reproduction degree before and after stabilization. This addition will confirm that order reduction does not occur. revision: yes

Circularity Check

0 steps flagged

No circularity: new algorithmic construction validated by independent numerical experiments

full rationale

The paper introduces a knot-removal stabilization strategy together with exact interlevel transfers based on knot insertion for B-spline interpolation in a multilevel fast integral method. All central claims (robust high-order accuracy, preservation of null space up to machine precision, O(N) complexity) are derived from the explicit algorithmic steps and are then checked against external numerical benchmarks on canonical and realistic geometries. No step reduces a result to a fitted parameter renamed as prediction, a self-definition, or a load-bearing self-citation; the derivation remains self-contained against the stated numerical evidence.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract contains no explicit free parameters, axioms, or invented entities; the contribution is framed as an algorithmic stabilization technique.

pith-pipeline@v0.9.1-grok · 5761 in / 1048 out tokens · 25188 ms · 2026-06-27T05:45:20.615956+00:00 · methodology

0 comments
read the original abstract

We present a multilevel B-spline-based fast integral method for the solution of the electric field integral equation (EFIE), combining fast Fourier transformation (FFT)-compatible kernel interpolation with robust high-order interpolation. Existing FFT-accelerated global Lagrange-based approaches rely on equidistant interpolation points and can, therefore, suffer from Runge-type instabilities at high interpolation orders, limiting robust high-accuracy compression. In contrast, B-splines on equidistant knot vectors overcome these instabilities and enable robust high-order interpolation for accurate matrix compression. Replacing Lagrange interpolation by B-spline interpolation is, however, non-trivial: B-spline coefficients do not coincide with function values at the interpolation points, and the associated sampling matrices can become ill-conditioned. To address these challenges, we introduce a knot-removal stabilization strategy, combined with exact interlevel transfers based on knot insertion, yielding accurate, well-conditioned multilevel interpolation. Moreover, we propose a factorization strategy that preserves the null space of the scalar potential operator up to machine precision and is compatible with low-frequency preconditioning techniques. Numerical results for both canonical and realistic geometries demonstrate robust high-order interpolation without the breakdown observed for Lagrange-based approaches and confirm $\mathcal{O}(N)$ complexity.

Figures

Figures reproduced from arXiv: 2606.13482 by Bernd Hofmann, Danijel Juki\'c, Simon B. Adrian, Thomas F. Eibert.

Figure 2
Figure 2. Figure 2: Interpolation grid associated with the axis-aligned boxes of the test [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Interpolation points and support of a protruding RWG function. [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: (a) Set of 𝑠 = 5 parent Lagrange polynomials of order 𝑝 = 4 and (b) the two sets of children Lagrange polynomials of order 𝑝 = 4. O (𝑠 6𝑁) to O (𝑠 3 log(𝑠)𝑁) [28]. This improved complexity is maintained by our algorithm. To this end, interlevel transfers are incorporated alongside intralevel translations. The algorithm is formulated on a regular octree hierarchy [22] to exploit the translation invariance o… view at source ↗
Figure 5
Figure 5. Figure 5: (a) Set of 𝑠 = 6 parent B-splines of order 𝑝 = 3 and (b) the two sets of children B-splines. Notably, the number of B-splines remains constant, with 𝑠 = 6 B-splines in both the parent and child clusters for boxes with an edge length smaller than 0.5𝜆. 0 1 (a) 0 0.2 0.4 0.6 0.8 1 0 1 (b) [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: (a) Set of 𝑠 = 9 parent B-splines of order 𝑝 = 3 and (b) the two sets of 𝑠 = 6 children B-splines of order 𝑝 = 3. Notably, the number of B-splines increases in the parent cluster for boxes with an edge length greater than or equal to 0.5𝜆 Combining these requirements, the number of splines and buffer knots on the parent level is given by 𝑠 parent = 𝑠 child , (32) 𝑠 parent b = max  ⌈𝑠 child b /2⌉, 𝑠 parent… view at source ↗
Figure 7
Figure 7. Figure 7: Surface mesh of a sphere with the corresponding axis-aligned boxes [PITH_FULL_IMAGE:figures/full_fig_p007_7.png] view at source ↗
Figure 9
Figure 9. Figure 9: Norm of the inverse B-spline sampling matrix for (a) equidistant [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Geometric conventions for defining an RWG function associated with [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Imaginary part of the surface current density [PITH_FULL_IMAGE:figures/full_fig_p010_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Imaginary part of the surface current density [PITH_FULL_IMAGE:figures/full_fig_p011_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Relative spectral-norm error between the Lagrange-based fast approximation [PITH_FULL_IMAGE:figures/full_fig_p011_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Relative spectral-norm error between the B-spline-based fast approximation [PITH_FULL_IMAGE:figures/full_fig_p012_14.png] view at source ↗
Figure 16
Figure 16. Figure 16: Relative spectral-norm error between the Lagrange-based fast approx 2 [PITH_FULL_IMAGE:figures/full_fig_p012_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Relative spectral-norm error between the B-spline-based fast approx 2 [PITH_FULL_IMAGE:figures/full_fig_p013_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Relative spectral-norm error between the B-spline-based fast ap 2 [PITH_FULL_IMAGE:figures/full_fig_p013_18.png] view at source ↗
Figure 20
Figure 20. Figure 20: Aircraft: performance analysis of the B-splined based compression [PITH_FULL_IMAGE:figures/full_fig_p014_20.png] view at source ↗

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. An Explicit Higher-Order Dual Basis for a Multiplicatively Calder\'on Preconditioned Electric Field Integral Equation

    math.NA 2026-07 conditional novelty 7.0

    An explicit B-spline dual basis generalizes Buffa-Christiansen functions to arbitrary order and keeps Calderón-preconditioned EFIE GMRES iterations constant in unknowns and degree.

Reference graph

Works this paper leans on

47 extracted references · 2 canonical work pages · cited by 1 Pith paper

  1. [1]

    The fast multipole method for the wave equation: A pedestrian prescription,

    R. Coifman, V. Rokhlin, and S. Wandzura, “The fast multipole method for the wave equation: A pedestrian prescription,”IEEE Antennas Propag. Mag., vol. 35, no. 3, pp. 7–12, Jun. 1993

    1993

  2. [2]

    Multilevel fast mul- tipole algorithm for electromagnetic scattering by large complex ob- jects,

    J. Song, Cai-Cheng Lu, and Weng Cho Chew, “Multilevel fast mul- tipole algorithm for electromagnetic scattering by large complex ob- jects,”IEEE Trans. Antennas Propag., vol. 45, no. 10, pp. 1488–1493, Oct. 1997

    1997

  3. [3]

    W. C. Chew, J.-M. Jin, E. Michielssen, and J. Song, Fast and Efficient Algorithms in Computational Electromagnetics(Artech House Antennas and Propagation Library). Boston: Artech House, 2001

    2001

  4. [4]

    A multilevel matrix decomposition algorithm for analyzing scattering from large structures,

    E. Michielssen and A. Boag, “A multilevel matrix decomposition algorithm for analyzing scattering from large structures,”IEEE Trans. Antennas Propag., vol. 44, no. 8, pp. 1086–1093, Aug. 1996

    1996

  5. [5]

    A butterfly-based direct integral-equation solver using hierarchical LU factorization for analyz- ing scattering from electrically large conducting objects,

    H. Guo, Y. Liu, J. Hu, and E. Michielssen, “A butterfly-based direct integral-equation solver using hierarchical LU factorization for analyz- ing scattering from electrically large conducting objects,”IEEE Trans. Antennas Propag., vol. 65, no. 9, pp. 4742–4750, Sep. 2017

    2017

  6. [6]

    AIM: Adaptive integral method for solving largescale electromagnetic scattering and radiation problems,

    E. Bleszynski, M. Bleszynski, and T. Jaroszewicz, “AIM: Adaptive integral method for solving largescale electromagnetic scattering and radiation problems,”Radio Sc., vol. 31, no. 5, pp. 1225–1251, Sep. 1996

    1996

  7. [7]

    A three-dimensional adaptive integral method for scattering from structures embedded in layered media,

    K. Yang and A. E. Yilmaz, “A three-dimensional adaptive integral method for scattering from structures embedded in layered media,” IEEE Trans. Geosci. Remote Sens., vol. 50, no. 4, pp. 1130–1139, Apr. 2012

    2012

  8. [8]

    A more scalable and efficient parallelization of the adaptive integral method—Part I: Algorithm,

    F. Wei and A. E. Yilmaz, “A more scalable and efficient parallelization of the adaptive integral method—Part I: Algorithm,” IEEE Trans. Antennas Propag., vol. 62, no. 2, pp. 714–726, Feb. 2014

    2014

  9. [9]

    High scalability FMM-FFT electromagnetic solver for supercomputer systems,

    J. M. Taboada, L. Landesa, F. Obelleiro, J. L. Rodriguez, J. M. Bertolo, M. G. Araujo, J. C. Mouriño, and A. Gomez, “High scalability FMM-FFT electromagnetic solver for supercomputer systems,”IEEE Antennas Propag. Mag., vol. 51, no. 6, pp. 20–28, Dec. 2009

    2009

  10. [10]

    MLFMA-FFT parallel algorithm for the solution of extremely large problems in electromagnetics,

    J. M. Taboada, M. G. Araujo, F. O. Basteiro, J. L. Rodriguez, and L. Landesa, “MLFMA-FFT parallel algorithm for the solution of extremely large problems in electromagnetics,”Proc. IEEE, vol. 101, no. 2, pp. 350–363, Feb. 2013

    2013

  11. [11]

    The adaptive cross approxi- mation algorithm for accelerated method of moments computations of EMC problems,

    K. Zhao, M. Vouvakis, and J.-F. Lee, “The adaptive cross approxi- mation algorithm for accelerated method of moments computations of EMC problems,”IEEE Trans. Electromagn. Compat., vol. 47, no. 4, pp. 763–773, Nov. 2005

    2005

  12. [12]

    Constructing nested bases approximations from the entries of non-local operators,

    M. Bebendorf and R. Venn, “Constructing nested bases approximations from the entries of non-local operators,”Numer. Math., vol. 121, no. 4, pp. 609–635, Aug. 2012

    2012

  13. [13]

    Fast nested cross approximation algorithm for solving large-scale electromagnetic problems,

    Y. Zhao, D. Jiao, and J. Mao, “Fast nested cross approximation algorithm for solving large-scale electromagnetic problems,”IEEE Trans. Microw. Theory Tech., vol. 67, no. 8, pp. 3271–3283, Aug. 2019

    2019

  14. [14]

    The incomplete adaptive cross ap- proximation for the fast construction ofℋ2-matrices and its application to the electric field integral equation for electrically small problems,

    J. M. Tetzner and S. B. Adrian, “The incomplete adaptive cross ap- proximation for the fast construction ofℋ2-matrices and its application to the electric field integral equation for electrically small problems,” IEEE Trans. Antennas Propag., pp. 1–1, 2026

    2026

  15. [15]

    On the fast matrix multiplication in the boundary element method by panel clustering,

    W. Hackbusch and Z. P. Nowak, “On the fast matrix multiplication in the boundary element method by panel clustering,”Numer. Math., vol. 54, no. 4, pp. 463–491, Jul. 1989

    1989

  16. [16]

    A precorrected-FFT method for electrostatic analysis of complicated 3-D structures,

    J. Phillips and J. White, “A precorrected-FFT method for electrostatic analysis of complicated 3-D structures,”IEEE Trans. Comput.-Aided Des. Integr. Circuits Syst., vol. 16, no. 10, pp. 1059–1072, Oct. 1997. JUKIĆet al.: A STABILIZED MULTILEVEL B-SPLINE-BASED FAST INTEGRAL METHOD FOR THE SOLUTION OF THE ELECTRIC FIELD INTEGRAL EQUATION15

    1997

  17. [17]

    Accelerated Cartesian expansions – A fast method for computing of potentials of the form𝑅−𝜈 for all real 𝜈,

    B. Shanker and H. Huang, “Accelerated Cartesian expansions – A fast method for computing of potentials of the form𝑅−𝜈 for all real 𝜈,” J. Comput. Phys., vol. 226, no. 1, pp. 732–753, Sep. 2007

    2007

  18. [18]

    Low-frequency fast inhomogeneous plane- wave algorithm (LF-FIPWA),

    L. Jiang and W. Chew, “Low-frequency fast inhomogeneous plane- wave algorithm (LF-FIPWA),”Microw. Opt. Technol. Lett., vol. 40, no. 2, pp. 117–122, Jan. 20, 2004

    2004

  19. [19]

    An efficient implementation of the com- bined wideband MLFMA/LF-FIPWA,

    D. Wulf and R. Bunger, “An efficient implementation of the com- bined wideband MLFMA/LF-FIPWA,”IEEE Trans. Antennas Propag., vol. 57, no. 2, pp. 467–474, Feb. 2009

    2009

  20. [20]

    A nondirective plane wave MLFMA stable at low frequencies,

    I. Bogaert, J. Peeters, and F. Olyslager, “A nondirective plane wave MLFMA stable at low frequencies,”IEEE Trans. Antennas Propag., vol. 56, no. 12, pp. 3752–3767, Dec. 2008

    2008

  21. [21]

    A low frequency stable plane wave addition theorem,

    I. Bogaert and F. Olyslager, “A low frequency stable plane wave addition theorem,”J. Comput. Phys., vol. 228, no. 4, pp. 1000–1016, Mar. 2009

    2009

  22. [22]

    ℋ2-matrix approximation of integral operators by interpolation,

    W. Hackbusch and S. Börm, “ℋ2-matrix approximation of integral operators by interpolation,”Appl. Numer. Math., vol. 43, no. 1–2, pp. 129–143, Oct. 2002

    2002

  23. [23]

    Approximation of integral operators by variable-order interpolation,

    S. Börm, M. Löhndorf, and J. M. Melenk, “Approximation of integral operators by variable-order interpolation,”Numer. Math., vol. 99, no. 4, pp. 605–643, Feb. 2005

    2005

  24. [24]

    A new multilevel Green’s function interpolation method for large-scale low-frequency EM simu- lations,

    H. G. Wang, C. H. Chan, and L. Tsang, “A new multilevel Green’s function interpolation method for large-scale low-frequency EM simu- lations,”IEEETrans.Comput.-AidedDes.Integr.CircuitsSyst. ,vol.24, no. 9, pp. 1427–1443, Sep. 2005

    2005

  25. [25]

    The implementation of multilevel Green’s function interpolation method for full-wave electromagnetic problems,

    H. Wang and C. Chan, “The implementation of multilevel Green’s function interpolation method for full-wave electromagnetic problems,” IEEE Trans. Antennas Propag., vol. 55, pp. 1348–1358, Jun. 1, 2007

    2007

  26. [26]

    Comparison of interpolating functions and interpolating points in full-wave multilevel Green’s function in- terpolation method,

    Yan Shi and Chi Hou Chan, “Comparison of interpolating functions and interpolating points in full-wave multilevel Green’s function in- terpolation method,”IEEE Trans. Antennas Propag., vol. 58, no. 8, pp. 2691–2699, Aug. 2010

    2010

  27. [27]

    An improved full-wave multilevel Green’s function interpolation method with RBF-QR technique for fast field evaluation,

    P. Zhao, C. H. Chan, and G. Wang, “An improved full-wave multilevel Green’s function interpolation method with RBF-QR technique for fast field evaluation,”IEEE Access, vol. 5, pp. 10241–10249, 2017

    2017

  28. [28]

    Low-frequency surface integral equation solution by multilevel Green’s function interpolation with fast Fourier transform acceleration,

    D. T. Schobert and T. F. Eibert, “Low-frequency surface integral equation solution by multilevel Green’s function interpolation with fast Fourier transform acceleration,”IEEE Trans. Antennas Propag., vol. 60, no. 3, pp. 1440–1449, Mar. 2012

    2012

  29. [29]

    Fast integral equation solution by multilevel Green’s function interpolation combined with multilevel fast multipole method,

    D. T. Schobert and T. F. Eibert, “Fast integral equation solution by multilevel Green’s function interpolation combined with multilevel fast multipole method,”IEEE Trans. Antennas Propag., vol. 60, no. 9, pp. 4458–4463, Sep. 2012

    2012

  30. [30]

    Investigating a B-spline based acceleration method for the electric field integral equation,

    D. Jukić, B. Hofmann, T. F. Eibert, and S. B. Adrian, “Investigating a B-spline based acceleration method for the electric field integral equation,” inProc. IEEE AP-S/URSI 2025, Ottawa, Canada, Jul. 2025, p. 3148

    2025

  31. [31]

    Investigation of a fast Fourier transform accelerated interpolatoryℋ2-method for the electric field integral equation,

    D. Jukić, B. Hofmann, T. F. Eibert, and S. B. Adrian, “Investigation of a fast Fourier transform accelerated interpolatoryℋ2-method for the electric field integral equation,” inInt. Conf. Electromagn. Adv. Appl. (ICEAA), Palermo, Italy, Sep. 2025, pp. 386–390

    2025

  32. [32]

    A broadband stabilized multilevel fast multipole method leveraging B-splines,

    D. Jukić, B. Hofmann, T. F. Eibert, and S. B. Adrian, “A broadband stabilized multilevel fast multipole method leveraging B-splines,” in Proc. Kleinheubach Conf. 2025, Miltenberg, Germany, Sep. 2025, pp. 1–4

    2025

  33. [33]

    H-field, E-field, and combined- field solutions for conducting bodies of revolution,

    J. R. Mautz and R. F. Harrington, “H-field, E-field, and combined- field solutions for conducting bodies of revolution,”Int. J. Electron. Commun., vol. 32, no. 4, pp. 157–164, 1978

    1978

  34. [34]

    Electromagnetic scattering by surfaces of arbitrary shape,

    S. Rao, D. Wilton, and A. Glisson, “Electromagnetic scattering by surfaces of arbitrary shape,”IEEE Trans. Antennas Propag., vol. 30, no. 3, pp. 409–418, May 1982

    1982

  35. [35]

    Electromagnetic integral equations: Insights in conditioning and pre- conditioning,

    S. B. Adrian, A. Dély, D. Consoli, A. Merlini, and F. P. Andriulli, “Electromagnetic integral equations: Insights in conditioning and pre- conditioning,”IEEE Open J. Antennas Propag., vol. 2, pp. 1143–1174, 2021

    2021

  36. [36]

    Barycentric Lagrange interpolation,

    J.-P. Berrut and L. N. Trefethen, “Barycentric Lagrange interpolation,” SIAM Rev., vol. 46, no. 3, pp. 501–517, Jan. 2004

    2004

  37. [37]

    De Boor, A Practical Guide to Splines: With 32 Figures(Applied Mathematical Sciences), revised ed

    C. De Boor, A Practical Guide to Splines: With 32 Figures(Applied Mathematical Sciences), revised ed. New York: Springer, 2001

    2001

  38. [38]

    L. A. Piegl and W. Tiller, The NURBS Book(Monographs in Visual Communications), 2nd ed. Berlin ; New York: Springer, 1997

    1997

  39. [39]

    R. A. Horn and C. R. Johnson,Topics in Matrix Analysis, Transferred to digital printing. Cambridge: Cambridge Univ. Press, 2010

    2010

  40. [40]

    Schnelle Lösung elektromagnetischer Integralgleichun- gen durch kombinierte Interpolation und Multipol-Faktorisierung Greenscher Funktionen,

    D. Schobert, “Schnelle Lösung elektromagnetischer Integralgleichun- gen durch kombinierte Interpolation und Multipol-Faktorisierung Greenscher Funktionen,” Technische Universität München, Munich, Germany, 2013

    2013

  41. [41]

    L. L. Schumaker, Spline Functions: Basic Theory(Cambridge Mathe- matical Library), 3rd. ed. Cambridge: Cambridge Univ. Press, 2010

    2010

  42. [42]

    Loop-star decomposition of basis functions in the dis- cretization of the EFIE,

    G. Vecchi, “Loop-star decomposition of basis functions in the dis- cretization of the EFIE,”IEEE Trans. Antennas Propag., vol. 47, no. 2, pp. 339–346, Feb. 1999

    1999

  43. [43]

    Loop-star and loop-tree decompositions: Analysis and efficient algorithms,

    F. P. Andriulli, “Loop-star and loop-tree decompositions: Analysis and efficient algorithms,”IEEE Trans. Antennas Propag., vol. 60, no. 5, pp. 2347–2356, May 2012

    2012

  44. [44]

    A multiplicative Calderon preconditioner for the electric field integral equation,

    F. P. Andriulli, K. Cools, H. Bagci, F. Olyslager, A. Buffa, S. Chris- tiansen, and E. Michielssen, “A multiplicative Calderon preconditioner for the electric field integral equation,”IEEE Trans. Antennas Propag., vol. 56, no. 8, pp. 2398–2412, Aug. 2008

    2008

  45. [45]

    Praktische Verfahren der Gleichungsauflösung,

    R. V. Mises and H. Pollaczek-Geiringer, “Praktische Verfahren der Gleichungsauflösung,”Z. Angew. Math. Mech., vol. 9, no. 1, pp. 58–77, Jan. 1929

    1929

  46. [46]

    13, 2026.doi: 10.5281/ZENODO.3865279

    krcools, sbadrian, azuccott, et al., Krcools/BEAST.jl: V2.9.0, ver- sion v2.9.0, Zenodo, Feb. 13, 2026.doi: 10.5281/ZENODO.3865279

    work page doi:10.5281/zenodo.3865279 2026

  47. [47]

    Hofmann, NURBS.jl, version v0.8.0, Zenodo, Nov

    B. Hofmann, NURBS.jl, version v0.8.0, Zenodo, Nov. 4, 2024.doi: 10.5281/zenodo.14031791

    work page doi:10.5281/zenodo.14031791 2024