How statement and proof provenance work
The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.
- Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
- AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
- AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.
These labels describe origin, not correctness: citations and verification chips remain separate evidence.
Planar Laplace multipole and local translation identities
Statement
For the P2M, M2M, M2L, L2L and L2P formulas of Fixed-order planar Laplace fast multipole procedure, P2M gives the source moments about a leaf center, M2M gives exactly the same moments about a parent center, and the finite M2L coefficients are the degree- through degree- Taylor coefficients of the order- source multipole polynomial about the target center. L2L is exact re-centering of the resulting degree- polynomial, and L2P evaluates it exactly. For every define ; the algorithm stores only . For a far pair and in the target cell, the fixed-branch source potential also has the absolutely convergent outer expansion .
Facts & Assumptions
Given: A far source/target cell pair, its fixed analytic logarithm, finite source moments, and the order- translation formulas.
Finite planar quadtree and separated interaction lists: for a far pair with centers , radii , and , one has and one analytic logarithm is fixed on a disk containing all pair differences.
Fixed-order planar Laplace fast multipole procedure: P2M stores ; M2M uses the binomial shift over nonempty children (empty children have zero moments); M2L uses the displayed finite local coefficients; L2L uses the binomial re-centering formula over nonempty-cell edges; L2P evaluates the local polynomial and its derivative.
Proof
For , the geometric series converges uniformly on each closed disk . Integrating it from to gives the analytic branch identity .
For a nonempty child cell of , write and apply the finite binomial theorem to its th power. Summing over the sources in each nonempty child and then over the nonempty children gives for every ; empty children contribute zero. This is exactly the M2M pass.
If a parent local polynomial is , put and expand each power. The coefficient of is ; hence L2L exactly re-centers the polynomial without adding truncation error.
For a far pair , put with and . Then . Substituting into the branch-normalized series from step 1.1 and summing the finite source set gives . The branch is because both and lie in its disk and the difference of their logarithms is the series that vanishes at . The convergence is absolute because .
For a far pair put , and . By [F1], . The branch-normalized logarithm therefore has expansion .
Expand by the binomial theorem and set . For , the coefficient of is , since . The terms with are for and for ; the terms with are . Summing over sources replaces it by , and retaining yields precisely the M2L formulas in [F2].
Substitution of a target point into the re-centered degree- polynomial is exactly L2P, and termwise differentiation gives its derivative formula. These are finite polynomial identities and require no additional convergence assumption.
Depends on
Used by
Dependency tree · one level
2 results within one dependency step of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.
Sources
- Leslie Greengard and Vladimir Rokhlin, A Fast Algorithm for Particle Simulations, Journal of Computational Physics 73 (1987), 325–348 (standard reference, not scraped)
- Michael O’Neil, Integral Equations and Fast Algorithms, complete course lecture notes (2017) (standard reference, not scraped)
- James Demmel, Fast Hierarchical Methods for the N-body Problem, Part 2, CS267 Lecture 25 (1996) (standard reference, not scraped)