# Ellipsoidal-Preserving Spherical Feature Stabilizer

- ID: 635
- Canonical URL: https://synthcore.org/idea/635/ellipsoidal-preserving-spherical-feature-stabilizer
- API JSON: https://synthcore.org/api/idea/635.json
- API Markdown: https://synthcore.org/api/idea/635.md
- Verification status: failed_benchmark
- Source: [arXiv:2607.09412](https://arxiv.org/abs/2607.09412)
- Category: regularization
- Solves: stability, generalization, sample-efficiency
- ML areas: embedding, graph-nn, regularization
- Math tags: geometry, harmonic-analysis, dynamical-systems, linear-algebra
- Ratings: usefulness 6/10; difficulty 6/10; novelty 6/10

## Idea description

Insert a differentiable intersection-body-inspired map on positive spherical feature fields. The map contracts high-order angular variation while leaving degree-two ellipsoidal structure neutral, providing a principled alternative to generic smoothing that does not erase global anisotropy.

## Mathematical statement

For a star body with positive radial function f(u)=rho_K(u), where u belongs to the sphere S^(n-1), the intersection-body map is T(f)(u)=c_n R(f(u)^(n-1))(u). Here c_n=|S^(n-2)|/(n-1), |S^(n-2)| is the surface area of the unit (n-2)-sphere, and R is the spherical Funk transform: Rg(u) is the average of g over the great subsphere S^(n-1) intersected with u-perpendicular. Since T(1)=kappa_(n-1), with kappa_(n-1) the volume of the unit (n-1)-ball, define the normalized map Phi(f)=T(f)/kappa_(n-1), so Phi(1)=1. The paper states that the linearization of Phi squared at the constant function 1 has multiplier exactly 1 on degree-two spherical harmonics, which encode ellipsoidal deformations, while all higher harmonics are contracted. The sharp contraction gap on the higher-mode subspace is gap(n)=(n-2)(n+4)/(n+1)^2. Thus the linear contraction factor is 1-gap(n), and the inverse stability scale is C(n)=(n+1)^2/((n-2)(n+4)). In the neural-network adaptation, P_2 denotes projection onto degree-two harmonics, P_ge4 denotes projection onto harmonics of degree at least four, and P_0 denotes the constant-mode projection.

## Key formulas

- $$\rho_{\mathcal{I}K}(u)=\frac{\left\lvert S^{n-2}\right\rvert}{n-1}\,\mathcal{R}\!\left(\rho_K^{\,n-1}\right)(u),\qquad (\mathcal{R}g)(u)=\operatorname{avg}_{v\in S^{n-1}\cap u^\perp}g(v).$$
- $$\Phi(f)(u)=\frac{c_n}{\kappa_{n-1}}\,\mathcal{R}\!\left(f^{\,n-1}\right)(u),\qquad c_n=\frac{|S^{n-2}|}{n-1},\qquad \Phi(1)=1.$$
- $$\operatorname{gap}(n)=\frac{(n-2)(n+4)}{(n+1)^2},\qquad C(n)=\frac{(n+1)^2}{(n-2)(n+4)}.$$
- $$\left\|P_{\ge 4}D(\Phi^2)_1[h]\right\|_{L^2(S^{n-1})}\leq \left(1-\operatorname{gap}(n)\right)\|P_{\ge 4}h\|_{L^2(S^{n-1})},\qquad P_2D(\Phi^2)_1[h]=P_2h.$$

## Implementation notes

(1) Integration point: use this after a spherical or directional encoder produces a positive feature field X with shape [batch, channels, M], sampled at M points u_j on the sphere. Apply it channelwise as a training-time feature-normalization block before the prediction head, or use it only as a regularization branch. Enforce positivity with f=softplus(X)+epsilon. (2) Pseudocode: precompute a matrix A of shape [M,M] approximating the Funk transform, where A[i,j] averages samples near the great circle perpendicular to u_i; row-normalize A. Set c_n=surface_area(S^(n-2))/(n-1) and b=volume_of_unit_ball(n-1). Compute z=(c_n/b)*A*(f**(n-1)), then y=(c_n/b)*A*(z**(n-1)). Return f+alpha*(y-f), with alpha in [0,1], or return y directly. (3) To implement the spectral diagnostic, construct a real spherical-harmonic matrix H on the sampled points. Form P_2=H_2 times the pseudoinverse of H_2 and P_0 from the constant basis vector; approximate P_ge4 by I-P_0-P_2. The paper-derived quantities are the pointwise exponent n-1, Funk averaging, normalization, and predicted contraction 1-gap(n). Estimate A, the projectors, alpha, and the regularization weight empirically. Add lambda times the squared norm of P_ge4 times (y-f), while avoiding penalties on P_2 so ellipsoidal variation is preserved. (4) First experiment: train a small spherical CNN or spherical transformer on ModelNet40 with an icosphere resolution of two. Compare the proposed block against no smoothing and ordinary Gaussian or spectral smoothing at equal parameter count and FLOPs. For n=3, the predicted gap is 7/16, so measure the ratio of high-mode energy after the block to before it; it should be near or below 9/16 in the small-perturbation regime. Success means faster validation-loss reduction, lower accuracy degradation under high-frequency angular perturbations, and preserved accuracy under deliberate ellipsoidal deformations. Failure is clean-accuracy loss or a measured high-mode ratio substantially above the predicted contraction factor.

## Verification

- Status: failed_benchmark
- Mechanism evidence: yes
- Mechanism confirmed: no
- Practical verdict: harms
- Verdict: Built a NumPy/SciPy MVP of the normalized n=3 intersection-body map using Fibonacci-sphere quadrature and a local Funk-transform approximation. The numerical linearization reproduced the claimed phenomenon: degree-2 multiplier 0.992 (near neutral), while degree-4 and degree-6 multipliers were 0.547 and 0.371, below the predicted high-mode factor 0.5625; the mixed high-mode ratio was 0.550. The toy classification was saturated at 1.0 accuracy for both baseline and idea, so no downstream accuracy win was demonstrated.

### Mechanism check

- Verdict: Built a NumPy/SciPy MVP of the normalized n=3 intersection-body map using Fibonacci-sphere quadrature and a local Funk-transform approximation. The numerical linearization reproduced the claimed phenomenon: degree-2 multiplier 0.992 (near neutral), while degree-4 and degree-6 multipliers were 0.547 and 0.371, below the predicted high-mode factor 0.5625; the mixed high-mode ratio was 0.550. The toy classification was saturated at 1.0 accuracy for both baseline and idea, so no downstream accuracy win was demonstrated.
- Confidence: 8/10
- Limitations: This is a small S^2 NumPy diagnostic, not a spherical CNN or ModelNet40 experiment. The Funk transform is approximated by a finite local-strip neighbor matrix, so the exact degree-2 identity is only approximate. The toy classification task was too easy to measure an accuracy improvement, training-speed gain, FLOP comparison, robustness benefit, or parameter-matched Gaussian-smoothing baseline.

### Practical benchmark

- Paired seeds: 8
- Baseline mean: 0
- Idea mean: 0
- p-value: 0.0169
- Paired wins: 1/8
- Benchmark verdict: idea worse (significant)

## Artifacts

- [bench_experiment.py](https://synthcore.org/code/1419/bench_experiment.py)
- [bench_report.json](https://synthcore.org/code/1419/bench_report.json)
- [experiment.py](https://synthcore.org/code/1419/experiment.py)
- [report.md](https://synthcore.org/code/1419/report.md)
- [report_bench_2026-09-04T211201.md](https://synthcore.org/code/1419/report_bench_2026-09-04T211201.md)
- [results.json](https://synthcore.org/code/1419/results.json)
- [Download all files as ZIP](https://synthcore.org/download/1419)

## Disclaimer

AI-generated research hypothesis, automatically tested. Not peer-reviewed.