From ca6aeabbe6c3628d898a8db97904c26c72620d25 Mon Sep 17 00:00:00 2001 From: Dylan Roy <249851642+d6g8k5htny-coder@users.noreply.github.com> Date: Mon, 5 Oct 2026 16:17:15 -0700 Subject: [PATCH] Add source-bound IBA2-012 finite-radius mixed-soft cluster candidate OpenAI / GPT-6 Astra Pro author continuation iba2-012-publication-audit-20261005-r3, acting for Dylan Roy. Four-file additive packet only; P and SC exact inputs retained. Includes full norm-tail argument and unequal-soft-scale corollary. Author-side candidate, nonauthor review required; scientific effect NONE. Pickup main#259/6005046599. --- .../iba2_012_finite_radius_20261005/PROOF.md | 253 ++++++++++++++++++ .../iba2_012_finite_radius_20261005/README.md | 29 ++ .../SOURCES.json | 57 ++++ .../test_finite_radius_cluster.py | 179 +++++++++++++ 4 files changed, 518 insertions(+) create mode 100644 reviews/iba2_012_finite_radius_20261005/PROOF.md create mode 100644 reviews/iba2_012_finite_radius_20261005/README.md create mode 100644 reviews/iba2_012_finite_radius_20261005/SOURCES.json create mode 100644 reviews/iba2_012_finite_radius_20261005/test_finite_radius_cluster.py diff --git a/reviews/iba2_012_finite_radius_20261005/PROOF.md b/reviews/iba2_012_finite_radius_20261005/PROOF.md new file mode 100644 index 000000000..5ed684c9b --- /dev/null +++ b/reviews/iba2_012_finite_radius_20261005/PROOF.md @@ -0,0 +1,253 @@ +# IBA2-012: a finite-radius additional-soft cluster estimate + +**Object:** IBA2-012-FINITE-R-CLUSTER-20261005-R3-v1.1. +**Author:** OpenAI / GPT-6 Astra Pro, continuation session `iba2-012-publication-audit-20261005-r3`, acting for Dylan Roy. +**Disposition:** NEW AUTHOR-SIDE PROOF CANDIDATE, conditional on the exact imported source interfaces below; nonauthor mathematical review required. +**Scientific effect:** NONE. No existing status, theorem, proof, register, review, flag or prize is changed. This is an additive publication of the preceding local proposal, not an accepted theorem or a review of its own mathematics. The author is exposed to the predecessor audit and other agents' reports; independence credit is zero. + +## 1. What this adds and what it does not + +SC §4 already establishes qualitative exclusion of further soft transverse eigenvalues from the fixed-radius, rescaled nonempty cluster measure. Codex `/root`'s main#259 comments 6003533425 and 6003547392 isolate and quantify the corresponding **limiting** measure: q additional hard eigenvalues in `(0,eta]` cost `eta^[q(q+7)/2]`. That work is credited, not reclaimed. + +A limiting estimate does not by itself allow `eta = r`. The present proposal proves an original-law **finite-r** bound, including norm tails, so this substitution becomes valid at fixed positive compact marks and fixed observation radius. It is separate from the preceding author's small-k cusp estimate (main#259/6004452810): that estimate concerns a different normalization and range of marks. The present argument does not consume that new, not-yet-reviewed lemma. + +This is not an exhaustive IBA2-012 closure, a rate of convergence of the full cluster distribution, a sharp lower bound, an all-mark result, or an elder matching theorem. No uniformity is claimed as the observation radius, dimension, volume, or marks escape their specified sets. + +## 2. Exact bindings and consumed interfaces + +Repository: `d6g8k5htny-coder/Math-`. Source cut: + +`04f47774bace871e3ad7dedc409e04f0280dfdd7`. + +| Tag | Path | Complete Git blob | Consumed material | +|---|---|---|---| +| P | `imports/lifetime_parent_20260925/UNIFORM_MATRIX_CAP_AND_LIFETIME.md` | `dfed3b8d318a3ab1950957f393307733a4bef3f2` | §3 (3.5); §4 full-field regression and conditional moments; §5 endpoint bounds and full normalizer (5.5); §6 typed determinant inequality (6.2); §7 Lebesgue spectral Jacobian | +| SC | `frontiers/spectral_cluster_closure_20260929/PROOF.md` | `16c56821b52fd76b0be791622b9c3809eafde75a` | §3 deterministic vector exclusion (11); §4 is the qualitative predecessor, not a quantitative premise | + +The consumed interfaces were re-read at this source cut and have the same complete Git blobs as the preceding e31b6ad8 and c08e8359 cuts. The core finite-radius proof is retained; v1.1 adds the unequal-scale corollary below and publication/reproduction details. The new proof does not need C6's all-order factorial moment theorem, remote kernels, the cubic classifier, or a full-variation convergence rate. A source-interface read is not represented as a fresh independent review of all ancestors. + +Immutable sources: +- https://github.com/d6g8k5htny-coder/Math-/blob/04f47774bace871e3ad7dedc409e04f0280dfdd7/imports/lifetime_parent_20260925/UNIFORM_MATRIX_CAP_AND_LIFETIME.md +- https://github.com/d6g8k5htny-coder/Math-/blob/04f47774bace871e3ad7dedc409e04f0280dfdd7/frontiers/spectral_cluster_closure_20260929/PROOF.md + +## 3. Precise statement + +Fix `d >= 3`, `m=d-1`, `L>0`, a compact birth set B, a compact positive mark interval `[k_-,k_+]` with `k_->0`, and a fixed `R>=1`. Use P's exact variance-one periodized field, original Gaussian pin law Q, typed maximum/saddle weight W, and **full** normalizer Z. All frames are allowed, with constants uniform on the compact frame space. + +Write `A_o = D_y^2 f(0)` and order the eigenvalues of `-A_o` as + + h_1 <= h_2 <= ... <= h_m. + +The h_i may be negative at finite r. Let N_R count additional critical points in the physical ball of radius Rr about the midpoint, with heights strictly between the pin heights, excluding the pins. Choose `K >= 1` as a fixed dimension-dependent multiple of `1+||f||_{C^4}` large enough to bound the operator derivatives in P and SC. + +For every fixed real `a>=0`, and every integer `1<=q<=m-1`, set + + beta_q = 4q + q(q-1)/2 = q(q+7)/2. + +There are `C0` such that, for `00, h_{q+1}<=eta] + <= C r^3 (eta+r)^beta_q. (FR) + +The event does not require h_1,...,h_{q+1} to be positive. Hence it also covers the narrower SC strip with positive hard eigenvalues, or `|h_{q+1}|<=eta`. + +The q=0 version is the usual fixed-R `O(r^3)` occurrence bound. For d=2 there is no additional transverse direction, so the q>=1 assertion is empty. + +## 4. Proof + +**Charts and measurability.** Reduce r_* so r<=1 and the closed radius-Rr ball and its surrounding cylinder embed in one torus chart. For fixed pins and marks, the event N_R>0 is Borel: exhaust the punctured closed ball by compact sets at distance at least 1/n from the pins and require the critical value to lie in [b-kr^3+1/n,b-1/n]. Existence of a zero gradient satisfying these closed conditions is closed in C^1, since its location lies in a fixed compact set; the countable union is precisely the nonempty event. Eigenvalue and norm conditions are Borel as well. No extra zero-count moment theorem is needed. + +### 4.1 Original-law moment and density facts + +On typed support, set `B_M=-D_y^2f(M)>0`, with ordered eigenvalues + + 0=z_*r^2` and `W/r^2<=CK^(2d)`. + +### 4.2 Convert near occurrence into a first-eigenvalue width + +Let + + E_large = {rK > k_-/(4R)}. + +Off E_large, SC (11) applies. If there is any additional critical point in the ball, there is one in the surrounding SC cylinder. The contrapositive of its exclusion criterion gives + + h_1 <= (8RK+56RK^2/k)r <= C_R rK^2. + +The Hessian Lipschitz bound and the min-max eigenvalue inequality give, for every i, + + |lambda_i-h_i| <= CrK. + +Therefore, on `{N_R>0,h_{q+1}<=eta}` outside E_large, + + lambda_1 <= C_R rK^2, + lambda_{q+1} <= eta+CrK. (2) + +This step is the additional localization not supplied by a bare determinant-layer statement. It uses no hard-eigenvalue inverse. + +### 4.3 Retain all endpoint determinant factors + +P (6.2), with `M3<=K` and a larger fixed constant, yields + + W/r^2 <= C K^2 prod_{i=1}^m lambda_i(lambda_i+c rK). + (3) + +In particular the saddle factor is retained even where a transverse eigenvalue is smaller than r. Replacing `lambda_i+c rK` by `lambda_i` would be invalid in that range. + +Split `K` into bands `[t,2t)`, `t=2^j`, j>=0. Let + + delta=c_1 rt, + A_t=c_2 rt^2, + B_t=eta+c_3 rt, + +where constants are enlarged so that `delta<=B_t` and (2) is contained in `lambda_1<=A_t`, `lambda_{q+1}<=B_t` on that band. + +By (1), for an arbitrarily large fixed s, + + P_Q(K>=t | B_M) <= C_s t^(-s)(1+||B_M||)^s. (4) + +After bounding `K^(a+2)` by `(2t)^(a+2)`, conditioning on B_M, and multiplying its density, the remaining polynomial from (4) must be absorbed **before** dropping any small-eigenvalue Gaussian factor: + + (1+||B_M||)^s exp(-c||B_M||^2) + <= C_s exp(-(c/2)||B_M||^2). (5) + +The constant may depend on the chosen s. Crucially, no power of t depending on s remains. This ordering prevents a circular choice of the tail moment. + +### 4.4 Spectral integration with two different widths + +Put `p=q+1` and `v_p=p(p-1)/2`. The Lebesgue spectral Jacobian is the Vandermonde times finite angular measure, independent of whether the Gaussian density is rotationally invariant. + +On the actual ordered integration domain, all of the first p eigenvalues are at most B_t. Their internal Vandermonde is therefore at most `B_t^v_p`. Every cross difference between a soft and a remaining hard eigenvalue is bounded by the hard eigenvalue; absorb these and all hard-hard differences into a fixed polynomial in the remaining eigenvalues. The hard determinant factors obey + + lambda_i(lambda_i+delta) <= C t(1+lambda_i)^2, + +because r<=1 and t>=1. The hard Gaussian integrals are thus at most `C t^(m-p)`. + +Only after those bounds are imposed on the actual domain, extend the nonnegative remaining soft integrals to independent intervals: lambda_1 in `[0,A_t]` and each of the other q soft variables in `[0,B_t]`. This extension remains valid even if A_t>B_t: the Vandermonde estimate was made before extension, not asserted on the enlarged domain. + +The first integral is + + int_0^{A_t} x(x+delta) dx + = A_t^3/3+delta A_t^2/2 <= C r^3 t^6. (6) + +Each other integral is at most + + int_0^{B_t} x(x+delta) dx <= (5/6)B_t^3, (7) + +using delta<=B_t. The numerical 5/6 may be absorbed into C. No dimension-dependent singular integral remains. + +Since + + 3q+v_p = 3q+q(q+1)/2 = beta_q, + B_t <= C t(eta+r), + +the contribution of the jth band, after dividing by the full Z and using (3)-(7), is at most + + C_s r^3(eta+r)^beta_q + t^[a+8+(m-p)+beta_q-s]. (8) + +Choose an integer `s>a+9+(m-p)+beta_q`. The dyadic sum converges. Its constant depends on the fixed parameters, q and a, but not on eta or r. This proves (FR) off E_large. + +### 4.5 The large-norm branch is paid for, not discarded + +For any M>0, the full normalizer and original Q moments give + + E_{Q^W}[K^a ; E_large] + <= C E_Q[K^(a+2d); K>k_-/(4Rr)] + <= C_M r^M. (9) + +Choose an integer M>=3+beta_q. Since eta+r>=r and r<=1, + + r^M <= r^(3+beta_q) <= r^3(eta+r)^beta_q. + +Adding (9) proves (FR), including eta=0. This finishes the analytic argument. + +## 5. Quantitative consequences + +For eta>=r, (FR) is `O(r^3 eta^beta_q)`. For eta<=r, the estimate saturates at `O(r^[3+beta_q])`. There is no claim of an extra eta power below the perturbation scale r. + +| Additional soft directions q | Total transverse corank addressed | beta_q | Bound at eta=r | +|---:|---:|---:|---:| +| 1 | 2 | 4 | O(r^7) | +| 2 | 3 | 9 | O(r^12) | +| 3 | 4 | 15 | O(r^18) | +| 4 | 5 | 22 | O(r^25) | + +Rows require `d>=q+2`. The exponent equals the compact p-soft count `3p+p(p-1)/2` at p=q+1, but now the count has an original-law proof with a nonempty spatial event and all derivative tails retained. + +More generally, for eta=r^alpha with alpha>0, + + E_{Q^W}[K^a ; N_R>0, h_{q+1}<=r^alpha] + = O(r^[3+beta_q min(alpha,1)]). (10) + +This is a joint event/derivative-weight estimate, not a rate for arbitrary test functions of the entire cluster process. Exact-count events `N_R=j` are subevents and inherit the same upper bound. Inserting a factor `N_R^j` would require a separate count-moment argument; it is not silently included. + +The estimate does not control an effective quartic-zero sector with a hard transverse Hessian, repeated nonzero hard eigenvalues, remote collisions, the limit R->infinity, or k->0. Those need their own source-bound statements. The preceding small-k cusp lemma and this fixed-positive-k cluster lemma cannot be combined into an all-mark theorem by relabeling constants. + +## 6. Why limiting mass alone cannot justify the diagonal + +For 00 is an event only. No unexplained determinant or factorial moment is added. +8. **Scope:** fixed d,L,R and positive compact marks throughout; all status changes withheld. + +The smallest useful nonauthor review is the chain SC(11) -> (2), the matrix-conditioned dyadic integration (3)-(8), and exception absorption (9), against the two exact blobs. No full re-audit of C6/remote ancestors is required for this proposition because they are not consumed. + +The accompanying standard-library exact checker tests the exponent ledger and soft polynomial integrals, including degree-changing negative controls. These are algebra controls, not formal verification of Gaussian regression or the continuum exclusion argument. + +## 8. Unequal soft scales: the same finite-radius argument + +**Corollary FR-mixed.** Keep every hypothesis of (FR). Fix `q` with `1<=q<=m-1`, put `p=q+1`, and choose deterministic thresholds `0<=eta_2<=...<=eta_p<=1`. Then, with C independent of all these thresholds and r, + + E_{Q^W}[K^a; N_R>0, h_j<=eta_j for every 2<=j<=p] + <= C r^3 prod_{j=2}^p(eta_j+r)^(j+2). (FR-mixed) + +The remaining eigenvalues have no hard lower bound. The all-soft case `p=m` is included. This is still a nonempty event bound, not a count-moment or convergence-rate theorem. + +*Proof.* In each norm band `[t,2t)`, keep `lambda_1<=A_t=c_2 r t^2` from (2) and set `B_{j,t}=eta_j+c_3 rt`, increasing c_3 so `delta<=B_{j,t}`. The actual ordered domain satisfies `lambda_j<=B_{j,t}`. Bound each internal Vandermonde factor by its larger eigenvalue, so + + prod_{1<=i=3+beta_q`, is absorbed because `prod(eta_j+r)^(j+2)>=r^beta_q`. This proves the result, including zero thresholds. + +For positive exponents with `eta_j=r^alpha_j` (ordered consistently), + + contribution = O(r^[3+sum_{j=2}^p(j+2)min(alpha_j,1)]). + +For example `q=2`, `eta_2=r^(1/2)` and `eta_3=r^(1/4)` give `O(r^(25/4))`; equal thresholds recover (FR). The factors have different powers: exchanging two exponents is not innocuous. This corollary adds hierarchical softening to the equal-width slice, not a general classification of arbitrary effective-jet or collision sectors. + +## 9. External reconnaissance boundary + +Consensus Primary again returned its monthly quota error in this continuation; no Consensus paper was retrieved or used. Public primary metadata for Armentano–Azais–Leon, arXiv:2304.07424, confirms that the current landing page is v4, while the project cites v3. No version substitution is made here. Azais–Delmas, arXiv:1911.02300v3, studies isotropic critical-point correlations; it is not a replacement for this exact-periodic determinant-weighted pin law. This proposition imports only P and SC as listed above and makes no literature novelty claim. + +## 10. Publication and reproduction + +Run `python -B -S test_finite_radius_cluster.py` and `python -B -O -S test_finite_radius_cluster.py` from this packet. The twelve test methods include the retained nine equal-scale controls and three mixed-scale checks. They are exact finite algebra tests, not a simulation, Lean replay or proof of the imported Gaussian interfaces. The newly added helper tests were run before implementing the helpers and failed for their two missing definitions; their successor run is recorded separately. No shared workflow is changed and the existing repository CI does not automatically run this standalone packet. Nonauthor analytic review and applicable current-candidate hosted checks remain required before integration. diff --git a/reviews/iba2_012_finite_radius_20261005/README.md b/reviews/iba2_012_finite_radius_20261005/README.md new file mode 100644 index 000000000..13d039726 --- /dev/null +++ b/reviews/iba2_012_finite_radius_20261005/README.md @@ -0,0 +1,29 @@ +# Finite-radius additional-soft cluster bounds + +Author-side mathematical candidate for IBA2-012. Scientific effect **NONE**. +Actual performer: OpenAI / GPT-6 Astra Pro, continuation `iba2-012-publication-audit-20261005-r3`, acting for Dylan Roy. This is the same author continuation as the preceding local proof, not an independent review. + +For fixed dimension, torus side, observation radius and compact positive marks, PROOF.md proves the proposed bound + + E_{Q^W}[K^a; N_R>0, h_{q+1}<=eta] <= C r^3 (eta+r)^{q(q+7)/2}. + +It also proves an unequal-scale refinement, with factors (eta_j+r)^{j+2}. The result is a finite-radius event/derivative-weight estimate. It does not establish a full cluster convergence rate, count-weighted factorial bound, all-mark statement, or general IBA2-012 closure. The separate small-k cusp proposal is not consumed. + +## Reproduce + +From this directory: + +```sh +python -B -S test_finite_radius_cluster.py +python -B -O -S test_finite_radius_cluster.py +``` + +Each runs 12 test methods and prints an exact algebra ledger. Normal and optimized stdout must agree. The tests check finite polynomial integration, exponent accounting, domain-extension bounds, large-norm absorption algebra, the limiting-measure counterexample, and unequal-scale accounting. Degree-changing variants are tested within the suite. No Lean execution, random-field simulation or analytic proof acceptance follows from these tests. + +SOURCES.json records the three other packet files by bytes, SHA-256 and Git blob, and the exact P/SC interfaces. Its own blob is recorded in the PR. Only these four files are introduced. No shared workflow, formal manifest, parent proof, status, scientific flag, review or prize is changed. Existing repository CI does not execute this standalone suite automatically. + +## Review and integration + +Pickup: main#259 comment6005046599. Credit for the limiting additional-hard slice remains with comments6003533425/6003547392. The preceding local archive has SHA-256 f2bf763ea204fc293d64fbae2b9eb68e7e37102e66ce438cc859d2465296d562. + +Nonauthor review should cover SC(11) to endpoint localization, P's conditional moments/full normalizer, Vandermonde bounds before domain extension, norm-polynomial absorption before discarding Gaussian decay, the retained saddle r term, and the unequal-scale refinement. Preserve fixed d,L,R and k_->0 throughout. The author will not self-merge. Applicable exact-candidate checks, ownership reconciliation and a guarded integration remain required; a green check is not a theorem disposition. diff --git a/reviews/iba2_012_finite_radius_20261005/SOURCES.json b/reviews/iba2_012_finite_radius_20261005/SOURCES.json new file mode 100644 index 000000000..634011ece --- /dev/null +++ b/reviews/iba2_012_finite_radius_20261005/SOURCES.json @@ -0,0 +1,57 @@ +{ + "schema": 1, + "scientific_effect": "NONE", + "repository": "d6g8k5htny-coder/Math-", + "source_cut": "04f47774bace871e3ad7dedc409e04f0280dfdd7", + "files": [ + { + "path": "PROOF.md", + "bytes": 17362, + "git_blob": "5ed684c9bb446efc9bef5f3980717c6dc039abfd", + "sha256": "cad7fa81a70dc218286d7ec62f7a86daf337677d65398ff906c455157067f93f" + }, + { + "path": "test_finite_radius_cluster.py", + "bytes": 6977, + "git_blob": "41c8f24d21658ad152ed64e19cbae130f910c643", + "sha256": "9aca39263f02d4092a604fbf52e41c670c325aa0854871cd8b3f8249ec7ee4a5" + }, + { + "path": "README.md", + "bytes": 2556, + "git_blob": "13d0397269d8efa62f8e1c91bf785d7ec4415a9b", + "sha256": "c5465da6e0048de5f92f66a57e2b0fd94ad3368e449b613c3fb9ea931005271f" + } + ], + "sources": [ + { + "id": "P", + "path": "imports/lifetime_parent_20260925/UNIFORM_MATRIX_CAP_AND_LIFETIME.md", + "git_blob": "dfed3b8d318a3ab1950957f393307733a4bef3f2", + "consumed": "(3.5), (4.1)-(4.3), (5.1)-(5.5), (6.2); Lebesgue spectral change of variables in section7. No numerical or acceptance upgrade.", + "fresh_read_lines": [ + [ + 86, + 124 + ], + [ + 120, + 282 + ] + ] + }, + { + "id": "SC", + "path": "frontiers/spectral_cluster_closure_20260929/PROOF.md", + "git_blob": "16c56821b52fd76b0be791622b9c3809eafde75a", + "consumed": "Section3 deterministic exclusion (11) and its proof only; section4 credited as predecessor, not quantitative premise.", + "fresh_read_lines": [ + [ + 110, + 235 + ] + ] + } + ], + "evidence_limits": "Source hashes establish identity, not proof acceptance. Tests are exact finite algebra; no hosted or Lean result is asserted in this inventory. The preceding small-k author proposal is not a premise." +} diff --git a/reviews/iba2_012_finite_radius_20261005/test_finite_radius_cluster.py b/reviews/iba2_012_finite_radius_20261005/test_finite_radius_cluster.py new file mode 100644 index 000000000..41c8f24d2 --- /dev/null +++ b/reviews/iba2_012_finite_radius_20261005/test_finite_radius_cluster.py @@ -0,0 +1,179 @@ +"""Exact algebra checks for FR; not a Gaussian-field simulation or Lean proof.""" +from fractions import Fraction as F +from math import ceil +import json +from pathlib import Path +import unittest + + +def mul(a, b): + out = {} + for ka, va in a.items(): + for kb, vb in b.items(): + key = tuple(x+y for x, y in zip(ka, kb)) + out[key] = out.get(key, F(0)) + va*vb + return {k: v for k, v in out.items() if v} + + +def mono(n, i): + key = [0]*n + key[i] = 1 + return {tuple(key): F(1)} + + +def add(a, b, scale=1): + out = dict(a) + for key, value in b.items(): + out[key] = out.get(key, F(0)) + scale*value + return {k: v for k, v in out.items() if v} + + +def soft_integral(p, mutant=None): + """Integrate prod x_i(x_i+r) V(x) over 0 tuple[int, ...]: + """Powers for eigenvalues numbered 2,...,q+1, not a probability bound.""" + if type(q) is not int or q < 0: + raise ValueError('q must be a nonnegative integer') + return tuple(j+2 for j in range(2,q+2)) + + +def mixed_radial_exponent(alphas: tuple[F, ...]) -> F: + """Exact homogeneity of the proved upper envelope at eta_j=r**alpha_j.""" + if any(a <= 0 for a in alphas): + raise ValueError('all exponents must be positive') + return F(3)+sum((w*min(F(a),F(1)) for w,a in zip(mixed_power_vector(len(alphas)),alphas)),F(0)) + + +class FiniteRadiusAlgebra(unittest.TestCase): + def test_first_band_integral(self): + for r in (F(1,2),F(1,16),F(1,100)): + for t in (F(1),F(2),F(8)): + A=r*t*t; delta=r*t + actual=A**3/3+delta*A*A/2 + self.assertEqual(actual,r**3*(t**6/3+t**5/2)) + self.assertLessEqual(actual,F(5,6)*r**3*t**6) + + def test_other_band_integral(self): + for b in (F(1,100),F(1),F(9)): + for ratio in (F(0),F(1,4),F(1)): + delta=b*ratio + self.assertLessEqual(b**3/3+delta*b*b/2,F(5,6)*b**3) + + def test_exact_full_soft_homogeneity(self): + for p in range(1,6): + got=soft_integral(p) + self.assertEqual(set(got), {3*p+p*(p-1)//2}) + self.assertGreater(next(iter(got.values())),0) + self.assertEqual(soft_integral(1),{3:F(5,6)}) + + def test_diagonal_and_codimension_ledgers(self): + for q in range(1,13): + p=q+1; beta=q*(q+7)//2 + self.assertEqual(beta,3*q+p*(p-1)//2) + self.assertEqual(3+beta,3*p+p*(p-1)//2) + self.assertEqual(p*(p+1)//2,p+p*(p-1)//2) + self.assertEqual([3+q*(q+7)//2 for q in range(1,5)],[7,12,18,25]) + + def test_dyadic_exponent_can_be_made_summable(self): + for m in range(2,10): + for q in range(1,m): + for a in (F(0),F(1,2),F(7)): + beta=q*(q+7)//2; p=q+1 + moment=ceil(a+10+m-p+beta) + exponent=a+8+m-p+beta-moment + self.assertLessEqual(exponent,-2) + + def test_large_norm_tail_absorption(self): + for q in range(1,5): + beta=q*(q+7)//2; M=3+beta + for r in (F(1,2),F(1,4),F(1,16)): + for eta in (F(0),r*r,r,F(1)): + self.assertLessEqual(r**M,r**3*(eta+r)**beta) + + def test_limiting_measure_does_not_supply_diagonal_rate(self): + ratios=[] + for r in (F(1,2),F(1,4),F(1,8),F(1,16)): + mass=r**4/4+r**2 + self.assertEqual(mass/r**4,F(1,4)+1/(r*r)) + ratios.append(mass/r**4) + self.assertTrue(all(x