Interactive Lab + FAQ

Interactive Lab + FAQ

Generative Radial Geometry (GRG)

A discrete–exact grammar for curves and surfaces using turns, chords, and radial projections. π is not imported as a magic constant; it emerges from dyadic steps and converges in the classical limit.

.

Turns-first definitions

We discretize the unit turn; dyadic refinement keeps operations exact on the mesh and drives accuracy.

$$ b \in \mathbb{N}_{\geq 0}, \qquad N = 2^{b} $$
$$ \Delta u = 2^{-(b+2)} \quad \text{(quarter sampling)} \qquad \Rightarrow \qquad \Delta u = \frac{1}{4N}. $$

For a circle of radius \(R\), the inscribed polygon at \(N\) turns has:

$$ P_{N} = 2 N R \sin\!\left(\frac{\pi}{N}\right), \qquad A_{N} = \frac{1}{2} N R^{2} \sin\!\left(\frac{2\pi}{N}\right). $$
Dyadic polygon and chords Circle with an N=16 inscribed polygon, dyadic rays, and a highlighted chord.
Turns, chords, and dyadic rays. Crisp vectors at any zoom.

Interactive Lab

Switch presets and adjust the dyadic exponent \(b\). For the star case \(r(\theta)=1+a\cos(k\theta)\), the exact area is \(A=\pi\!\left(1+\frac{a^{2}}{2}\right)\).

N (turns)
—
Perimeter (approx)
—
Area (approx)
—
Perimeter true
—
Area true
—
Errors (|P|, |A|)
—

Notes. Circle & atlas use exact finite-step formulas (chords) and converge monotonically. Star-shape uses inscribed polygons; C-shape uses odd–even radial sweep for area (perimeter shown exact for reference).

Critical FAQ (condensed)

Q: Is GRG just numerical integration with a trick?
A: For circles with dyadic bisections steps are exact at finite \(N\) (chords & projections). The limit recovers classical values; nothing is imported.

Q: Where does \(\pi\) come from?
A: From the discrete perimeter/area growth of dyadic polygons. \(\pi\) appears as a rotational-symmetry shadow as \(N\to\infty\).

Q: Non-convex shapes “jump across the outside”.
A: We prevent that with visibility-aware schemes: (i) odd–even radial sweep; (ii) radial atlas (locally star-shaped patches with ± signs for outer components/holes).

Q: Convergence?
A: \(O(N^{-2})\) for \(C^{2}\) boundaries; \(O(N^{-1})\) near corners (mitigated via adaptive dyadics). Inscribed perimeter/area are monotone.

Q: Why is GRG pedagogically useful?
A: You can teach Euclidean geometry from bisections and chords; radians & trig tables emerge rather than being prerequisites.

Technical Appendix — Non-Convex & Complex Parametrizations

  1. Star-Check. If a center sees the whole boundary (star-shaped), sample \(\theta\) in dyadic turns and use \(r(\theta)\). Inscribed area/perimeter are monotone \(\uparrow\) with \(O(N^{-2})\).
  2. Odd–Even radial sweep. For each \(\theta_{k}\), intersect the ray with the boundary, sort radii \(r_{k,j}\), and add inside segments \([r_{k,1},\, r_{k,2}],\ [r_{k,3},\, r_{k,4}], \dots\).
  3. Radial atlas (± signs). Decompose into locally star-shaped patches; outer components \(+\), holes \(−\). Sum areas with signs and perimeters by absolute length. Same \(O(N^{-2})\) patch-wise.
  4. Robustness. \(C^{2} \Rightarrow O(N^{-2})\); corners \(\Rightarrow O(N^{-1})\) unless adaptively refined. Self-intersections: use winding rule instead of parity.

All computations use turns, chords, projections; operations are exact at finite dyadic steps for circular primitives and exact on the mesh for polygons.

Team

Roles are presented by responsibility area. Individual bios can be added later.

Team roles pictogram Minimal icons representing research, algorithms, systems, and editorial review.
Role-based structure. Short scope, clear ownership.

Leadership & Research

  • Sets the research agenda and publishes specifications with reproducible examples.
  • Maintains the “turns-first” glossary and consistency across demos and notes.
  • Approves public statements and abstracts before release.

Methods & Algorithms

  • Designs visibility-aware procedures (odd–even sweep, radial atlas with ± signs).
  • Owns convergence studies and baselines for circles, ellipses, and polygons.
  • Runs complexity/accuracy benchmarks; documents limits and mitigations.

Systems & Frontend

  • Implements interactive canvases; keeps demos accessible and performant.
  • Maintains build tooling and test fixtures for numerical correctness.
  • Coordinates localization and analytics (non-invasive, anonymized).

Editorial Review & QA

  • Checks clarity, tone, and scope; enforces the evidence-led style.
  • Validates formulas and captions (MathJax + SVG semantics and accessibility).
  • Controls versioning of figures; ensures non-repetitive hero imagery.

Collaborations

We welcome collaborations with educators and researchers. A public contact form will be added once finalized.

.
Measuring Shape

Measuring Shape

Generative Radial Geometry (GRG): measuring shape to change perspective

This page presents concise formulas and proofs showing how GRG replaces constants like $\pi$ with shape measures (circles, hemispheres, ellipses, polygons) and with unit-sum quadratures. The goal is twofold: measurement precision and operational advantage (stability, variance/error reduction, traceability), even in asymmetric scenarios.

1) Principle

In GRG, many occurrences of $2\pi$, $4\pi$ or $\pi$ are not “mysteries” but the area/volume of a shape in the natural domain (e.g., directions, frequencies). By changing the measure on that domain (circular units: turn, hemi, sphere) the constant is absorbed and the integral becomes a normalized average.

2) Core theorems

Theorem 1 — Circular measure equivalence

Let $S^1$ be the unit circle of directions with arc-length measure $\lambda$ so that $\lambda(S^1)=2\pi$. Define the normalized measure $\mu=\lambda/(2\pi)$ with $\mu(S^1)=1$. Then for all $g\in L^1(S^1)$:

$$ \int_{S^1} g \, d\lambda \;=\; 2\pi \int_{S^1} g \, d\mu . $$

Proof. By definition of $\mu$. For $S^2$ the same holds with $\lambda(S^2)=4\pi$ and for the hemisphere with $2\pi$.

Theorem 2 — Anisotropic measure (no symmetry required)

Let $\rho(\theta)\ge 0$ with $\int \rho(\theta)\,d\theta = 1$ be an angular density (possibly non-uniform). Define the weighted measure $d\nu(\theta):=\rho(\theta)\,d\theta$. For all $g\in L^1$:

$$ \int g(\theta)\,\rho(\theta)\,d\theta \;=\; \mathbb{E}_{\theta\sim \rho}\!\big[g(\theta)\big]. $$

Interpretation. The “shape” of the domain is deformed according to $\rho$. Precision comes from encoding asymmetry in the measure, not from any pre-existing symmetry.

Theorem 3 — Radial quadratures with unit-sum weights

There exist sets of directions $\{\hat u_k\}_{k=1}^{N}\subset S^2$ and weights $w_k>0$ with $\sum_k w_k = 1$ such that, for every $g$ in a target family (e.g., polynomials in $\mu=\cos\theta$ up to degree $t$, or spherical harmonics up to degree $t$):

$$ \int_{S^2} g(\hat u)\, d\Omega \;=\; \sum_{k=1}^{N} w_k\, g(\hat u_k). $$

Note. Such cubatures (e.g., spherical t-designs or refined symmetric polyhedra) give controllable accuracy with $N=O(t^2)$.

3) Operational propositions: “removing $\pi$”

Proposition A — Fourier pair in turns (Hz), no $2\pi$

With frequency in cycles/s (Hz) and phase in turns, the Fourier pair

$$ X(f)=\int_{-\infty}^{\infty} x(t)\,e^{-i\,2 f t}\,dt , $$
$$ x(t)=\int_{-\infty}^{\infty} X(f)\,e^{+i\,2 f t}\,df . $$

is equivalent to the classical pair with $\omega=2\pi f$. Proof. Substitute $\omega=2\pi f$ and $d\omega=2\pi\,df$.

Proposition B — RC integral: rad/s vs Hz

$$ \int_{0}^{\infty} \frac{d\omega}{1+(\omega\tau)^2} \;=\; \frac{\pi}{2\,\tau} $$
$$ \Longleftrightarrow\quad \int_{0}^{\infty} \frac{df}{1+(2\pi f \tau)^2} \;=\; \frac{1}{4\,\tau} . $$

Proof. Use $\omega=2\pi f$. The constant $\pi$ is a unit-conversion factor.

4) Demonstrative examples (non–“pure geometry”, explicit shapes)

4.1 Radiometry (Lambert): the number is the area of a disk

For a perfectly diffuse surface,

$$ E \;=\; L \int_{\text{hemisphere}} \cos\theta \, d\Omega . $$

In spherical coordinates $d\Omega=\sin\theta\, d\theta\, d\varphi$, $\varphi\in[0,2\pi)$, $\theta\in[0,\tfrac{\pi}{2}]$. Then

$$ \int_{\text{hemisphere}} \cos\theta \, d\Omega \;=\; \int_0^{2\pi}\! d\varphi \int_0^{\pi/2}\! \cos\theta \sin\theta\, d\theta \;=\; 2\pi \cdot \tfrac{1}{2} \;=\; \pi . $$

Hence $E=\pi L$. With the “hemi” measure (hemisphere $=1$) the average is $\langle \cos\theta\rangle_{\text{hemi}} = \tfrac{1}{2}$ and, expressed per-hemi, $E_{\text{hemi}} = L \cdot \tfrac{1}{2}$. With the projected “Lambert” measure $d\Sigma = \tfrac{\cos\theta}{\pi}\, d\Omega$ (unit disk area), $\int d\Sigma=1$ and therefore $E = L$. In all cases the constant is a shape measure.

4.2 Signals/Control (RC): angle as shape

With the substitution $\omega\tau = \tan\theta$, the integral

$$ \int_0^\infty \frac{d\omega}{1+(\omega\tau)^2} $$

becomes a measure of an angle on the unit circle (arc from $0$ to $\tfrac{\pi}{2}$), which explains $\tfrac{\pi}{2}$. In “turn/Hz” units that angle is $\tfrac{1}{4}$, giving $\tfrac{1}{4\tau}$. Same content, different shape and measure.

4.3 Statistics (Cauchy as a projection of the circle)

If $\Theta \sim \operatorname{Unif}(-\tfrac{\pi}{2},\tfrac{\pi}{2})$ and $X=\gamma \tan\Theta$, then

$$ f_X(x) \;=\; \frac{1}{\pi \gamma}\,\frac{1}{1+(x/\gamma)^2} . $$

Interpretation: $X$ is the tangent projection of a point uniformly distributed on a circle. In “turn” units, let $\Phi \sim \operatorname{Unif}(-\tfrac{1}{4},\tfrac{1}{4})$ and $\Theta = 2\pi \Phi$; with respect to the turn-normalized angular measure (total length $=1$) the normalization constant is absorbed. Converting back to radians restores $1/(\pi\gamma)$.

5) Asymmetry: elliptical metrics and weighted measures

For anisotropic phenomena introduce an elliptical metric $M\succ 0$ and define $r_M(x)=\sqrt{x^\top M x}$. The “balls” of $r_M$ are ellipses, and many constants are absorbed by switching to $M$-radial coordinates. Any directional density $\rho(\theta)$ defines $d\nu=\rho\,d\theta$: integration becomes a unit-sum average without invoking symmetry.

6) Precise discrete re-determination (quadratures)

Design sets $\{(\hat u_k,w_k)\}$ with $\sum_k w_k=1$ that integrate exactly a useful family (e.g., $\{1,\cos\theta,\sin\theta,\cos2\theta,\sin2\theta\}$ for low-order anisotropic BRDFs, or spherical harmonics up to degree $t$). This yields a precise discrete scheme: constants like “$4\pi$” or “$2\pi$” become simply “sum of weights = 1”.

7) Practical conversions (GRG ↔ classical form)

Domain Classical form GRG form Observation
Radiometry (Lambert) $E=\pi L$ $E=L$ with unit projected measure $\pi$ = area of the projected disk
Signals (RC) $\displaystyle \int \frac{d\omega}{1+(\omega\tau)^2}=\frac{\pi}{2\tau}$ $\displaystyle \int \frac{df}{1+(2\pi f\tau)^2}=\frac{1}{4\tau}$ $\pi$ = rad↔turn conversion
Statistics (Cauchy) $\frac{1}{\pi\gamma}\frac{1}{1+(x/\gamma)^2}$ $\frac{1}{\gamma}\frac{1}{1+(x/\gamma)^2}$ with angle in turns $\pi$ = length of the angular domain
Fields/Fluxes $4\pi$ factors on $S^2$ Average over sphere with $\sum w_k=1$ $4\pi$ = area of the unit sphere

8) Measurable benefits

  • Numerical stability: fewer fixed large/small factors, better conditioning.
  • Variance reduction: with anisotropic $\rho$, MC error shrinks where it matters.
  • Traceability: results come pre-normalized (averages on 1); conversion back to standards is multiplication by $2\pi$ or $4\pi$.

9) Validation protocol (A/B)

  1. Select 3 cases (Lambert, RC, Cauchy/Maxwell).
  2. Define the GRG measure (turn/hemi/sphere or weighted $\rho$) and the quadratures (degree $t$).
  3. Compare to the classical baseline at equal budget (samples/time).
  4. Metrics: relative error, variance, conservation (energy/flux), stability along pipelines.
  5. Success threshold: ≥50% lower cost at equal error or ≤½ the error at equal cost in ≥2 out of 3 cases.

10) Conclusion

“Removing $\pi$” is not rhetoric: it is the outcome of measuring the domain in the right shape. GRG provides a measurement language (circular units and unit-sum quadratures) that yields cleaner formulas and, crucially, quantifiable benefits. This is a change of perspective that delivers progress: measurable improvement, not stubbornness.

Paradigm Shift

Paradigm Shift

A Constructive View of Radial Geometry

Generative Radial Geometry studies measures through explicitly defined radial constructions. Its current reference is a finite binary sector framework with internal chord and external tangent elements.

Start with Definitions

Choose the radius, a full Turn, and a sector resolution. Define the geometric elements before choosing a kernel or discussing accuracy.

Separate Three Questions

  • Construction: which elements and operations define the object?
  • Measurement: which kernel produces the finite estimate?
  • Verification: under which assumptions do bounds and convergence hold?

Current Results

Read the current formulation and analytical verification. The results concern regular symmetric sectors in exact arithmetic. Computational efficiency and extensions require separate evidence.

Research History

Earlier projection, Cartesian trapezoidal, and Turns approaches remain documented. They should be compared within their own scope. This page replaces the earlier slogan-based presentation with a statement of the constructive method and its verification requirements.

Wedge Visualization (Interactive)

Wedge Visualization (Interactive)

GRG — Wedge Visualization: Interactive Uniform vs Adaptive Discretization

This mini-demo provides concrete, reproducible numbers for wedge-stratified GRG. Adjust the controls to compare a uniform wedge count against an adaptive profile on the outer annulus, and see how certified corridor half-widths contract roughly like 1/\sqrt{W_k} for comparable ring budgets.

How the Demo Works (plain language)

  • Dyadic edges. Rings follow turns-first dyadic refinement with b∈ℕ≥0, N=2^b, and quarter sampling Δu=2^{-b-2}.
  • Corridor half-width. Each sector’s estimate has a certified interval [E^- , E^+]. The half-width is (E^+−E^-)/2; decisions stop when this width ≤ tolerance τ.
  • Field model. We synthesize an anisotropic field with mixed angular orders up to fourth order to stress the wedges; values are computed in-page so the numbers are concrete and repeatable.
$$ F(\rho,\theta)=A(\rho)\Big[1+\epsilon_2\cos 2(\theta-\theta_2)+\epsilon_3\cos 3(\theta-\theta_3)+\epsilon_4\cos 4(\theta-\theta_4)\Big], $$
$$ \textbf{Stop when}\quad \widehat{E}_b^{+}-\widehat{E}_b^{-}\le \tau,\qquad \text{with}\ \ \Delta u=2^{-b-2}. $$

Controls

Results — Per-Sector Estimates on the Outer Annulus

Hover or tap a bar to read exact values (estimate and corridor half-width). The polar plots show the directional pattern and corridor bands.

Uniform W — outer annulus sectors (bars = estimate, whiskers = corridor half-width).
Adaptive outer Wk — tighter corridors where anisotropy is richer.
Polar — Uniform W.
Polar — Adaptive Wk.

Concrete Numbers (auto-computed)

The table aggregates sector values (outer annulus) for the current settings. Replace later with CSVs from your scripts; the structure is the same.

Notes & Practical Guidance

  • When to prefer GRG. As angular richness grows or is localized, adaptive W_k exposes structure with linear O(K·W) cost and certified decisions via τ.
  • When low-order bases suffice. Near-radial scenes with mild 2-fold variation are well captured by a short angular basis; choose GRG when auditability (corridors) matters.
  • Next steps. Publish side-by-side benchmarks (Zernike, B-splines, Fourier–Bessel) with the same corridor-aware stopping rule and public scripts/datasets.
GRG — Wedge Visualization

GRG — Wedge Visualization

Wedge-Stratified GRG — Visualizing Per-Sector Results

This section visualizes the anisotropic, wedge-stratified example discussed earlier. We compare a uniform W=8 discretization against an adaptive scheme on the outer annulus \rho\in[0.7,1), where angular content is richer. The plots show per-sector estimates with certified corridor half-widths.

Setup Recap

Field with mixed angular orders up to fourth order:

$$ F(\rho,\theta)=A(\rho)\,\Big[1+\epsilon_2\cos 2(\theta-\theta_2)+\epsilon_3\cos 3(\theta-\theta_3)+\epsilon_4\cos 4(\theta-\theta_4)\Big], $$

Rings and wedges:

$$ \text{Rings }k=1..K,\ \text{wedges }w=1..W_k,\quad I_{k,w}=\int_{\Theta_w}\!\int_{u_{k-1}}^{u_k} F(\rho,\theta)\,\rho\,d\rho\,d\theta. $$
  • Rings: K=32 (dyadic edges, quarter sampling).
  • Uniform: W=8 for all rings.
  • Adaptive: inner W_k=8, mid W_k=12, outer W_k=16 (used for plots below).

Per-Sector Results on the Outer Annulus

The numbers are the same illustrative values used in the write-up; replace them with your measured outputs when available. Error bars display corridor half-widths.

Left: uniform W=8. Right: adaptive outer ring with W_k=16. Hover/tap a bar to read exact values.

Polar View (Directional Pattern)

A simple polar plot places each sector’s estimate radially. Shaded rings visualize the corridor range.

Polar plot of per-sector estimates. Left: W=8. Right: W_k=16.

Reading the Plots

  • Corridor effect: with higher angular resolution, corridor half-widths shrink (≈ 1/\sqrt{W_k} for comparable ring budgets), revealing higher-order structure.
  • Decision making: choose the smallest dyadic level b such that the corridor width is below the task tolerance \tau.
$$ \textbf{Stop when}\ \ \widehat{E}^{+}_{b}-\widehat{E}^{-}_{b}\le\tau. $$