Pith. sign in

REVIEW 2 major objections 5 minor 49 references

An explicit high-order dual basis generalizes Buffa-Christiansen functions and keeps Calderón-preconditioned EFIE iterations low and constant in mesh size and polynomial degree.

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 →

T0 review · grok-4.5

2026-07-13 06:16 UTC pith:E5VYK6PF

load-bearing objection Explicit arbitrary-order dual B-spline basis that finally makes multiplicative Calderón preconditioning practical for higher-order isogeometric EFIE, with iteration counts that stay flat in both N and p. the 2 major comments →

arxiv 2607.08848 v1 pith:E5VYK6PF submitted 2026-07-09 math.NA cs.NAphysics.comp-ph

An Explicit Higher-Order Dual Basis for a Multiplicatively Calder\'on Preconditioned Electric Field Integral Equation

classification math.NA cs.NAphysics.comp-ph MSC 65N3865R2078M15
keywords B-splinesCalderón preconditioningEFIEdual basisBuffa-Christiansen functionsisogeometric analysishigher-order basis functionsGMRES
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 electric field integral equation is a standard way to compute electromagnetic scattering from perfect conductors, but when it is discretized with higher-order basis functions its matrix becomes badly conditioned as the mesh is refined or the polynomial degree rises. Multiplicative Calderón preconditioning cures that pathology for the lowest-order Rao-Wilton-Glisson basis by pairing it with a carefully constructed dual basis (Buffa-Christiansen functions). This paper supplies the first explicit dual basis that works for arbitrary polynomial degree when the primal space is built from divergence-conforming B-splines. The dual functions are obtained by linear combination of a refined B-spline set whose knots are shifted so that the mixed Gram matrix stays well-conditioned. On spheres, cubes, plates, a spaceplane and a car the resulting preconditioned system converges in a few dozen GMRES iterations that do not grow with either the number of unknowns or the polynomial degree, thereby unlocking the accuracy advantages of higher-order bases without the usual iterative-cost penalty.

Core claim

The authors construct an explicit dual basis for any polynomial degree that is dual (in the Buffa-Christiansen sense) to a divergence-conforming B-spline primal basis on curvilinear quadrilateral patches; when this dual is used inside a multiplicative Calderón preconditioner the discrete operator remains spectrally well-behaved and the GMRES iteration count stays low and essentially independent of both the number of unknowns and the polynomial degree.

What carries the argument

The dual basis is formed by taking linear combinations, with coefficients in {1, 1/2}, of a refined B-spline space whose knot vectors are non-uniformly shifted so that Greville sites of the dual align with those of the primal; the refinement doubles the number of functions per parametric direction and exchanges the roles of loops and stars on the Greville mesh.

Load-bearing premise

The particular non-uniform knot shift that aligns Greville sites is assumed to keep the Gram-matrix condition number bounded for all practical mesh sizes and polynomial degrees; the paper only shows this numerically up to a few hundred unknowns per patch and degree four.

What would settle it

Compute the condition number of the mixed Gram matrix (or the GMRES iteration count of the preconditioned EFIE) for the same family of dual bases on a sequence of meshes with N_uv well beyond 100 and polynomial degrees p_uv greater than 4; unbounded growth would falsify the claim that the refinement rule is sufficient.

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 paper constructs an explicit higher-order dual basis for divergence-conforming B-spline (isogeometric) discretizations of the EFIE and uses it to form a multiplicative Calderón preconditioner. The dual functions are local linear combinations of refined Curry–Schoenberg B-splines on a Greville mesh that generalizes the barycentric refinement of Buffa–Christiansen functions; a non-uniform knot shift is introduced so that the mixed Gram matrix remains well-conditioned. The three classical dual-basis properties (divergence conformity, bounded Gram condition number, and vanishing of the dual-primal scalar-potential product via swapped loop-star association) are verified analytically. Numerical experiments on spheres, cubes, plates, a spaceplane and a car show that GMRES iteration counts stay low and essentially independent of both the number of unknowns and the polynomial degree.

Significance. An explicit dual basis of arbitrary polynomial degree removes a long-standing obstacle to higher-order multiplicative Calderón preconditioning of the EFIE. Because the construction recovers the classical BC functions at lowest order and inherits the loop-star structure already available for B-splines, it supplies a reusable building block for other dual-basis techniques (quasi-Helmholtz projectors, combined-field formulations, etc.). The numerical evidence that iteration counts remain constant with both h- and p-refinement is practically important: it finally lets higher-order bases deliver their theoretical accuracy advantage without a prohibitive growth in solver cost. The work is therefore a solid and timely contribution to computational electromagnetics and isogeometric boundary-element methods.

major comments (2)
  1. §III-C, Eqs. (26)–(33) and Figs. 7–8: the claim that the mixed Gram matrix remains well-conditioned for arbitrary N_uv and p_uv rests entirely on a numerically tuned non-uniform knot shift whose clipping parameter a is chosen ad hoc. While the reported range (N_uv up to ~100, p up to 4) is already useful, a short analytic argument or a sharper asymptotic bound would strengthen the central claim that the preconditioner is robust for all practical refinements.
  2. §IV-D and Fig. 18: the high-frequency complexification of the dual wavenumber employs a free parameter H that is selected by inspecting a histogram of mean curvature. The paper should either provide a reproducible, geometry-independent rule for choosing H or demonstrate that a simple default (e.g., the median or a fixed fraction of the maximum) already yields iteration counts comparable to the hand-tuned value.
minor comments (5)
  1. Abstract and throughout: “unknonws” is a recurring typographical error for “unknowns”.
  2. §III-B, Fig. 3 and surrounding text: the precise assignment of the six refined functions and the weights {1,1/2} for a general dual edge would be clearer if a short algorithmic description or a small table of local indices were added.
  3. §III-E: the remark that the primal and refined Bézier meshes are in general non-aligned is important for implementers; a brief note on how the sparse Gram matrix is still assembled efficiently (intersection mesh) would help reproducibility.
  4. References [26],[27] are conference abstracts of the same work; they could be moved to a “preliminary results” footnote to avoid self-citation inflation.
  5. Fig. 9 caption: the RWG triangulation is said to use 8100 functions; confirming that the same number is used for every p_uv would make the accuracy comparison more transparent.

Circularity Check

1 steps flagged

No significant circularity: dual-basis construction is explicit and independent of the numerical confirmation of constant GMRES iterations.

specific steps
  1. self citation load bearing [§III-A (primal basis) and §III-D (properties iii)]
    "we choose the same divergence conforming basis functions as introduced in [4], since we can leverage the loop-star decomposition of [4] and its specific graph properties. ... the relations TΦ,k Λ=0, ΛT TΦ,k=0, ... holding for the primal basis carry over to eTΦ,˜k Σ=0, ΣT eTΦ,˜k=0 ... for the dual basis"

    The dual-basis construction and the proof that property (iii) holds rely on the loop-star matrices and graph properties previously defined by the same authors. Those prior results are used as established tools rather than as circular definitions of the dual basis itself, so the circularity is minor and not load-bearing for the central claim of constant GMRES iterations.

full rationale

The paper constructs an explicit higher-order dual basis by linear combination of refined B-spline functions on a non-uniformly refined Greville mesh (Eqs. 18–33, Figs. 2–5), then verifies the three required properties (divergence conformity, well-conditioned Gram matrix, vanishing dual-primal scalar-potential product) from that construction and the interchange of loop/star roles (§III-B,D). The subsequent numerical results (§IV) are independent confirmation that the resulting Calderón-preconditioned system yields low, essentially constant GMRES iteration counts. Self-citations to the authors’ prior B-spline EFIE and loop-star work ([4] and related) supply the primal space and mapping matrices as established tools; they do not define the dual basis or force the iteration-count claim by construction. The only mild self-referential element is the non-uniform knot-shift rule (Eqs. 26–33), whose sufficiency for bounded Gram conditioning is shown only numerically (Figs. 7–8); that is a weakest-assumption issue, not a circular reduction of a claimed prediction to its inputs. Score 1 reflects ordinary, non-load-bearing self-citation of prior tools.

Axiom & Free-Parameter Ledger

2 free parameters · 3 axioms · 1 invented entities

The central claim rests on standard Calderón theory, the known B-spline primal space of the authors’ earlier work, and a small set of geometric and refinement choices introduced in this paper. No free parameters are fitted to the scattering data that demonstrate constant iteration counts; the only tunable quantities appear in the knot-refinement rule and in the optional high-frequency complexification of the dual wavenumber.

free parameters (2)
  • clipping parameter a in μ_a = a = p_u/(2N_u-1-p_u)
    Chosen as a = p_u/(2N_u-1-p_u) to keep knots well-separated near the boundary; not derived from a uniqueness theorem.
  • mean-curvature value H used in complexified dual wavenumber = 0.5 m^{-1}
    Selected as the mode of the curvature histogram (H=0.5 m^{-1}) after inspecting the space-plane geometry; other values give higher iteration counts.
axioms (3)
  • domain assumption Calderón identity T_k^{2} = -I/4 + K^{2} holds for the continuous EFIE operator on Lipschitz surfaces.
    Invoked in §II-B as the analytic foundation of multiplicative preconditioning; standard in the literature.
  • domain assumption The B-spline primal basis of Hofmann et al. (2024) admits a discrete loop-star decomposition associated with the Greville mesh.
    Used throughout §III to swap the roles of loops and stars for the dual basis; taken from the authors’ prior work.
  • ad hoc to paper Linear combinations of refined Curry-Schoenberg B-splines with coefficients in {1,1/2} (plus corner weights) remain divergence-conforming and produce a well-conditioned mixed Gram matrix under the proposed non-uniform refinement.
    The explicit dual construction of §III-B–C; verified numerically but not proved for arbitrary p and N.
invented entities (1)
  • explicit higher-order dual B-spline basis (generalized Buffa-Christiansen functions) no independent evidence
    purpose: Provide a dual space that satisfies the three algebraic properties required for a multiplicative Calderón preconditioner at arbitrary polynomial degree.
    Defined by local linear combinations of a refined B-spline basis on a non-uniformly refined knot vector; reduces to classical BC functions when p=1.

pith-pipeline@v1.1.0-grok45 · 23874 in / 2756 out tokens · 29320 ms · 2026-07-13T06:16:44.326070+00:00 · methodology

0 comments
read the original abstract

One of the most effective means to precondition the electric field integral equation (EFIE) discretized with Rao-Wilton-Glisson (RWG) functions is the multiplicative Calder\'on preconditioner employing Buffa-Christiansen (BC) functions as a basis dual to the RWG basis. It results in a formulation that is free from the dense-discretization and the low-frequency breakdown. To generalize the multiplicative Calder\'on preconditioner from the low-order BC and RWG basis to higher orders, we utilize B-spline-based basis functions and establish the first explicit high-order dual basis. It can be regarded as a generalization of the BC functions to arbitrary polynomial degrees and constitutes a fundamental building block for other approaches that rely on a dual basis. Numerical results for the obtained preconditioner demonstrate a low and constant number of generalized minimum residual (GMRES) iterations independent of the number of unknonws and the polynomial degree for canonical and realistic perfectly electrically conducting (PEC) scatterers; a key to enable the full potential of higher-order bases.

Figures

Figures reproduced from arXiv: 2607.08848 by Bernd Hofmann, Francesco P. Andriulli, Simon B. Adrian, Thomas F. Eibert.

Figure 2
Figure 2. Figure 2: Refined basis functions ˆ𝒇𝑖 𝑗 𝑑 and example of a basis function dual to the primal basis function shown in [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Forming dual basis functions from the refined basis: six refined [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: Functions 𝜇(𝑢) and 𝜇𝑎 (𝑢) used for the definition of the refined knot vectors 𝑈ˇ and 𝑈˜ , respectively. nature when combining the dual basis functions around the faces of the Greville mesh (which is essential for property iii) to hold also for open geometries): due to the Curry-Schoenberg normalization, the weights can be interpreted as the amount of charge flowing in and out of each Greville mesh face. He… view at source ↗
Figure 5
Figure 5. Figure 5: Forming dual basis functions on a patch touching the boundaries of [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: Conditioning of the Gram matrix for a closed surface for different [PITH_FULL_IMAGE:figures/full_fig_p006_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Conditioning of the Gram matrix for an open surface for different [PITH_FULL_IMAGE:figures/full_fig_p006_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Scattering of a plane wave from a sphere of [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 11
Figure 11. Figure 11: Condition number for a sphere of 𝑟s = 1 m discretized with 1452 unknowns for different frequencies. is computed for the electric field. In all cases, we have also computed the errors for the magnetic field and the far field (FF) and verified that they match the electric field error. As can be seen from [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 10
Figure 10. Figure 10: Scattering of a plane wave from a sphere of [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 13
Figure 13. Figure 13: Scattering of a plane wave from a square plate of edge length [PITH_FULL_IMAGE:figures/full_fig_p009_13.png] view at source ↗
Figure 12
Figure 12. Figure 12: Scattering of a plane wave from a cube of edge length [PITH_FULL_IMAGE:figures/full_fig_p009_12.png] view at source ↗
Figure 14
Figure 14. Figure 14: Model of a spaceplane with 81 NURBS patches, induced surface current densities and GMRES residual for the illumination with a plane wave, where the shuttle has a maximum extent of 0.08𝜆 and a discretization with 23 328 basis functions. (a) 𝑝𝑢𝑣 = 1 (b) 𝑝𝑢𝑣 = 2 (c) 𝑝𝑢𝑣 = 3 −20 −10 0 10 20 30 40 50 in m−1 (d) mean curvature 𝐻 0 5,000 10,000 15,000 20,000 25,000 30,000 35,000 40,000 100 10−1 10−2 10−3 10−4 10… view at source ↗
Figure 15
Figure 15. Figure 15: Scattering from the spaceplane for the illumination with a plane wave, where the shuttle has a maximum extent of [PITH_FULL_IMAGE:figures/full_fig_p010_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Model of a car with 1032 NURBS patches, induced surface current densities and GMRES residual for the illumination with a plane wave, where the car has a maximum extent of 8.3𝜆 and a discretization with 132 072 basis functions. such that the shuttle has a maximum extent of 0.08𝜆. The fast convergence of GMRES within 45 iterations, independent of the polynomial degree, compared to several thousand iteration… view at source ↗
Figure 17
Figure 17. Figure 17: Extract of the histogram for mean curvature [PITH_FULL_IMAGE:figures/full_fig_p011_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Number of iterations for the 21𝜆 spaceplane for different values of 𝐻 in 𝑘˜ = 𝑘 − j0.4𝑘 1/3𝐻2/3 for the Calderón preconditioner for 𝑝𝑢𝑣 = 2. Table II Solution times for the scattering from the car. 𝑝𝑢𝑣 = 1 𝑝𝑢𝑣 = 2 𝑝𝑢𝑣 = 3 𝑝𝑢𝑣 = 4 No preconditioner 2.3 h 2.75 h 11.3 h 37.0 h This work 17.8 min 17.8 min 17.8 min 17.8 min but explicitly constructed the matrices and performed the full matrix-vector products u… view at source ↗

discussion (0)

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

Reference graph

Works this paper leans on

49 extracted references · 1 linked inside Pith

  1. [1]

    Zur Formulierung eines allgemeinen Beugungsproblems durch eine Integralgleichung,

    A.-W. Maue, “Zur Formulierung eines allgemeinen Beugungsproblems durch eine Integralgleichung,”Zeitschrift Phys., vol. 126, pp. 601–618, Jul. 1949

  2. [2]

    Galerkin boundary element methods for electromagnetic scattering,

    A. Buffa and R. Hiptmair, “Galerkin boundary element methods for electromagnetic scattering,” inTopics in Computational Wave Propa- gation, M. Ainsworth, P. Davies, D. Duncan, B. Rynne, and P. Martin, Eds., ed. by T. J. Barth, M. Griebel, D. E. Keyes, R. M. Nieminen, D. Roose, and T. Schlick, vol. 31, Berlin, Heidelberg: Springer Berlin Heidelberg, 2003, ...

  3. [3]

    On a well-conditioned electric field integral operator for multiply connected geometries,

    F. P. Andriulli, K. Cools, I. Bogaert, and E. Michielssen, “On a well-conditioned electric field integral operator for multiply connected geometries,”IEEE Trans. Antennas Propag., vol. 61, no. 4, pp. 2077– 2087, Apr. 2013

  4. [4]

    Low-frequency stabilization for the B-spline-based isogeometric discretization of the electric field integral equation,

    B. Hofmann, M. Mirmohammadsadeghi, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “Low-frequency stabilization for the B-spline-based isogeometric discretization of the electric field integral equation,”IEEE Trans. Antennas Propag., vol. 72, no. 4, pp. 3558–3571, Apr. 2024

  5. [5]

    An excitation-aware and self-adaptive frequency normalization for low- frequencystabilizedelectricfieldintegralequationformulations,

    B. Hofmann, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “An excitation-aware and self-adaptive frequency normalization for low- frequencystabilizedelectricfieldintegralequationformulations,”IEEE Trans. Antennas Propag., vol. 71, no. 5, pp. 4301–4314, May 2023

  6. [6]

    A low-frequency stable, excitation agnostic discretization of the right- hand side for the electric field integral equation on multiply-connected geometries,

    B. Hofmann, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “A low-frequency stable, excitation agnostic discretization of the right- hand side for the electric field integral equation on multiply-connected geometries,”IEEE Trans. Antennas Propag., vol. 71, no. 12, pp. 9277– 9288, Dec. 2023

  7. [7]

    The construction of some efficient preconditioners in the boundary element method,

    O. Steinbach and W. Wendland, “The construction of some efficient preconditioners in the boundary element method,”Advances in Com- putational Mathematics, vol. 9, no. 1, pp. 191–216, Sep. 1, 1998

  8. [8]

    Operator preconditioning,

    R. Hiptmair, “Operator preconditioning,”Comp. Math. Appl., vol. 52, no. 5, pp. 699–706, Sep. 1, 2006

  9. [9]

    A dual finite element complex on the barycentric refinement,

    A. Buffa and S. Christiansen, “A dual finite element complex on the barycentric refinement,”Math. Comp., vol. 76, no. 260, pp. 1743–1769, Oct. 2007

  10. [10]

    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

  11. [11]

    An introduction to operator precondition- ing for the fast iterative integral equation solution of time-harmonic scattering problems,

    X. Antoine and M. Darbas, “An introduction to operator precondition- ing for the fast iterative integral equation solution of time-harmonic scattering problems,”Multiscale Sci. Eng., vol. 3, no. 1, pp. 1–35, Mar. 1, 2021

  12. [12]

    Quadrilateral barycentric basis functions for surface integral equations,

    R. Chang and V. Lomakin, “Quadrilateral barycentric basis functions for surface integral equations,”IEEE Trans. Antennas Propag., vol. 61, no. 12, pp. 6039–6050, Dec. 2013

  13. [13]

    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

  14. [14]

    Parametric mapping of vector basis functions for surface integral equation formulations,

    A. F. Peterson and K. R. Aberegg, “Parametric mapping of vector basis functions for surface integral equation formulations,”Appl. Comput. Electromagn. Soc. J. ACES, pp. 107–115, 1995

  15. [15]

    Higher order interpolatory vector bases for computational electromagnetics,

    R. Graglia, D. Wilton, and A. Peterson, “Higher order interpolatory vector bases for computational electromagnetics,”IEEE Trans. Anten- nas Propag., vol. 45, no. 3, pp. 329–342, Mar. 1997

  16. [16]

    A. F. Peterson,Mapped Vector Basis Functions for Electromagnetic Integral Equations. Morgan & Claypool Publishers, Dec. 2006

  17. [17]

    High- order div- and quasi curl-conforming basis functions for Calderón multiplicative preconditioning of the EFIE,

    F. Valdés, F. P. Andriulli, K. Cools, and E. Michielssen, “High- order div- and quasi curl-conforming basis functions for Calderón multiplicative preconditioning of the EFIE,”IEEE Trans. Antennas Propag., vol. 59, no. 4, pp. 1321–1337, Apr. 2011

  18. [18]

    On a novel Calderón preconditioning strategy based on high-order quasi- Helmholtz projectors,

    J. Bourhis, D. Franzo, A. Merlini, and F. P. Andriulli, “On a novel Calderón preconditioning strategy based on high-order quasi- Helmholtz projectors,” inProc. Int. Conf. Electromagn. Adv. Appl. (ICEAA), Sep. 2024, pp. 700–703

  19. [19]

    Preconditioning strategies and conformal discretizations empowered by high-order projectors,

    J. Bourhis, A. Merlini, and F. P. Andriulli, “Preconditioning strategies and conformal discretizations empowered by high-order projectors,” in URSI Int. Symp. Electrom. Theory (EMTS), Jun. 2025, pp. 1–4

  20. [20]

    Subdivision based isogeometric analysis technique for electric field integral equations for simply connected structures,

    J. Li, D. Dault, B. Liu, Y. Tong, and B. Shanker, “Subdivision based isogeometric analysis technique for electric field integral equations for simply connected structures,”J. Comput. Phys., vol. 319, pp. 145–162, Aug. 15, 2016

  21. [21]

    Mathematical analysis of variational isogeometric methods,

    L. B. Veiga, A. Buffa, G. Sangalli, and R. Vázquez, “Mathematical analysis of variational isogeometric methods,”Acta Numer., vol. 23, pp. 157–287, May 2014

  22. [22]

    An isogeometric boundary element method for electromagnetic scattering with compat- ible B-spline discretizations,

    R. N. Simpson, Z. Liu, R. Vázquez, and J. A. Evans, “An isogeometric boundary element method for electromagnetic scattering with compat- ible B-spline discretizations,”J. Comp. Physics, vol. 362, pp. 264–289, Jun. 1, 2018

  23. [23]

    J. Dölz, S. Kurz, S. Schöps, and F. Wolf, “Isogeometric boundary elements in electromagnetism: Rigorous analysis, fast methods, and HOFMANNet al.: AN EXPLICIT HIGHER-ORDER DUAL BASIS FOR A MULTIPLICATIVELY CALDERÓN PRECONDITIONED EFIE 13 examples,”SIAM J. Sci. Comput., vol. 41, no. 5, B983–B1010, Jan. 2019

  24. [24]

    A mixed finite element method for 2-nd order elliptic problems,

    P. A. Raviart and J. M. Thomas, “A mixed finite element method for 2-nd order elliptic problems,” inMathematical Aspects of Finite Element Methods, I. Galligani and E. Magenes, Eds., ser. Lecture Notes in Mathematics, Berlin, Heidelberg: Springer, 1977, pp. 292–315

  25. [25]

    Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement,

    T. J. R. Hughes, J. A. Cottrell, and Y. Bazilevs, “Isogeometric analysis: CAD, finite elements, NURBS, exact geometry and mesh refinement,” Comput. Methods Appl. Mech. Eng., vol. 194, no. 39, pp. 4135–4195, Oct. 1, 2005

  26. [26]

    A multiplicative Calderón preconditioner for the B-spline-based electric field integral equation,

    B. Hofmann, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “A multiplicative Calderón preconditioner for the B-spline-based electric field integral equation,” inProc. Int. Conf. Electromagn. Adv. Appl. (ICEAA), Lisbon, Portugal, Sep. 2024

  27. [27]

    B- spline-based higher-order discretization of the magnetic field integral equation,

    B. Hofmann, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “B- spline-based higher-order discretization of the magnetic field integral equation,”inProc.KleinheubacherTagung,Miltenberg,Germany,Sep. 2024

  28. [28]

    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, pp. 157–164, Apr. 1978

  29. [29]

    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., pp. 1143–1174, Oct. 2021

  30. [30]

    Magnetic and combined field integral equations based on the quasi- Helmholtz projectors,

    A. Merlini, Y. Beghein, K. Cools, E. Michielssen, and F. P. Andriulli, “Magnetic and combined field integral equations based on the quasi- Helmholtz projectors,”IEEE Trans. Antennas Propag., vol. 68, no. 5, pp. 3834–3846, May 2020

  31. [31]

    Préconditionneurs analytiques de type Calderón pour les formulations intégrales des problèmes de diffraction d’ondes,

    M. Darbas, “Préconditionneurs analytiques de type Calderón pour les formulations intégrales des problèmes de diffraction d’ondes,” Département de Génie Mathématique, Nat. Inst. Appl. Sci. Toulouse, Toulouse, France, 2004

  32. [32]

    Well-conditioned boundary integral equa- tion formulations for the solution of high-frequency electromagnetic scattering problems,

    Y. Boubendir and C. Turc, “Well-conditioned boundary integral equa- tion formulations for the solution of high-frequency electromagnetic scattering problems,”Computers & Mathematics with Applications, vol. 67, no. 10, pp. 1772–1805, Jun. 1, 2014

  33. [33]

    Time do- main Calderón identities and their application to the integral equation analysis of scattering by PEC objects Part I: Preconditioning,

    K. Cools, F. P. Andriulli, F. Olyslager, and E. Michielssen, “Time do- main Calderón identities and their application to the integral equation analysis of scattering by PEC objects Part I: Preconditioning,”IEEE Trans. Antennas Propag., vol. 57, no. 8, pp. 2352–2364, Aug. 2009

  34. [34]

    de Boor,A Practical Guide to Splines(Applied Mathematical Sciences 27), 1st ed

    C. de Boor,A Practical Guide to Splines(Applied Mathematical Sciences 27), 1st ed. New York, NY: Springer, Nov. 2001

  35. [35]

    Isogeometric methods for computational electromagnetics: B-spline and T-spline discretizations,

    A. Buffa, G. Sangalli, and R. Vázquez, “Isogeometric methods for computational electromagnetics: B-spline and T-spline discretizations,” J. Comput. Phys., Physics-compatible numerical methods, vol. 257, pp. 1291–1320, Jan. 15, 2014

  36. [36]

    Piegl,The NURBS Book

    L. Piegl,The NURBS Book. Berlin, Heidelberg: Springer Berlin Hei- delberg, 1995

  37. [37]

    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

  38. [38]

    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, 931 pp

  39. [39]

    A stabilized multilevel B-spline-based fast integral method for the solution of the electric field integral equation,

    D. Jukić, B. Hofmann, T. F. Eibert, and S. B. Adrian, “A stabilized multilevel B-spline-based fast integral method for the solution of the electric field integral equation,” Jun. 11, 2026. arXiv: 2606.13482, pre-published

  40. [40]

    Cools, S

    K. Cools, S. Adrian, and collaborators,BEAST.jl, version v2.9.0, Genève, Switzerland, Feb. 13, 2026. [Online]. Available: https://doi. org/10.5281/zenodo.3865279

  41. [41]

    Hofmann,NURBS.jl, version v0.8.0, Genève, Switzerland, Nov

    B. Hofmann,NURBS.jl, version v0.8.0, Genève, Switzerland, Nov. 4,

  42. [42]

    Available: https://doi.org/10.5281/zenodo.14031791

    [Online]. Available: https://doi.org/10.5281/zenodo.14031791

  43. [43]

    Quadrature forℎ 𝑝-Galerkin BEM in R3,

    S. A. Sauter and C. Schwab, “Quadrature forℎ 𝑝-Galerkin BEM in R3,”Numer. Math., vol. 78, no. 2, pp. 211–258, Dec. 1997

  44. [44]

    S. A. Sauter and C. Schwab,Boundary Element Methods. Berlin, Heidelberg: Springer Berlin Heidelberg, 2011

  45. [45]

    Jin,Theory and Computation of Electromagnetic Fields, Second edition

    J.-M. Jin,Theory and Computation of Electromagnetic Fields, Second edition. Hoboken, New Jersey: John Wiley & Sons, Inc, 2015

  46. [46]

    SphericalScattering: A Julia package for electromagnetic scattering from spherical objects,

    B. Hofmann, P. Respondek, and S. B. Adrian, “SphericalScattering: A Julia package for electromagnetic scattering from spherical objects,” J. Open Source Softw., vol. 8, no. 91, pp. 1–4, Nov. 2023

  47. [47]

    Low- frequency-stabilized electric field integral equation on topologically non-trivial geometries for arbitrary excitations,

    B. Hofmann, T. F. Eibert, F. P. Andriulli, and S. B. Adrian, “Low- frequency-stabilized electric field integral equation on topologically non-trivial geometries for arbitrary excitations,” inProc. IEEE Anten- nas Propag. Soc. Int. Symp. URSI Nat. Radio Sci. Meeting, Denver, CO, Jul. 2022

  48. [48]

    A numerical comparison of an isogeometric and a parametric higher order Raviart–Thomas approach to the electric field integral equation,

    J. Dölz, S. Kurz, S. Schöps, and F. Wolf, “A numerical comparison of an isogeometric and a parametric higher order Raviart–Thomas approach to the electric field integral equation,”IEEE Trans. Antennas Propag., vol. 68, no. 1, pp. 593–597, Jan. 2020

  49. [49]

    Preconditioning the EFIE on screens,

    R. Hiptmair and C. Urzúa-Torres, “Preconditioning the EFIE on screens,”Math. Models Methods Appl. Sci., vol. 30, no. 9, Sep. 2020