Continuous

Continuous

GRG — Continuous Model

Continuous Model: Radial Fields, Well-posedness & Stability

Category: Generative Radial Geometry (GRG)

We describe a continuous model that interprets GRG’s inside-out construction as integral curves of a center-seeded field on a domain Ω ⊂ ℝ². Let the seed be s ∈ Ω. For any point x ≠ s define the unit radial direction er(x) = (x−s)/‖x−s‖ and a tangential direction et(x) orthogonal to er. The guidance field blends radial and tangential components.

Radial Guidance Field

The model uses a field F: Ω\\{s} → ℝ² with radial bias

$$ \mathbf{F}(x) \;=\; a(x)\, \mathbf{e}_r(x) \;+\; b(x)\, \mathbf{e}_t(x), \qquad a(x) > 0 $$

Integral curves γ(τ) solve

$$ \dot{\gamma}(\tau) \;=\; \mathbf{F}\big(\gamma(\tau)\big), \qquad \gamma(0)=x_0,\quad x_0\in\Omega\setminus\{s\} $$

Each turn-angle u ∈ [0,1) selects an initial direction, producing a family of trajectories that reach the boundary. The discrete GRG construction samples these trajectories at uniform steps in turns.

Well-posedness

  • Existence & Uniqueness. If F is locally Lipschitz on Ω\\{s} then trajectories starting away from the seed have local existence and uniqueness while they remain in that domain. Global continuation and boundary arrival require additional hypotheses.
  • Termination. Trajectories stop when they hit ∂Ω. Dirichlet-type constraints prescribe a target boundary, while Neumann-type constraints control flux through ∂Ω.

Stability to Perturbations

If F has Lipschitz constant L, perturbing the seed by ε yields trajectories with deviation bounded by

$$ \|\gamma_\varepsilon(\tau) - \gamma(\tau)\| \;\le\; C\,\varepsilon\, e^{L\tau} $$

In practice we keep L moderate via smoothing of a(x), b(x) and by capping curvature of trajectories.

Regularization & Safety

  • Step lower bound. Minimum angular step Δu_{\min} to avoid local collapse.
  • Curvature damping. Add −λ\,κ(\gamma)\,\mathbf{e}_t to reduce spiraling near thin features.
  • Line-search. Backtrack the last advance if distance to boundary decreases below a threshold.

Multi-Seed and Holes

For non-simply connected domains or shapes with holes we blend fields from auxiliary seeds s_i with weights w_i(x):

$$ \mathbf{F}(x) \;=\; \sum_{i=1}^{m} w_i(x)\,\Big(a_i(x)\,\mathbf{e}_r^{(i)}(x) + b_i(x)\,\mathbf{e}_t^{(i)}(x)\Big), \qquad \sum_i w_i(x)=1,\ \ w_i(x)\ge 0 $$

Hole guards choose weights that vanish across cut-loci to prevent ambiguous paths.

Discrete Connection

Sampling each trajectory at uniform steps in turns yields the GRG primitives used for chord sums and radial strips. The continuous model justifies convergence as resolution increases. With turns-first sampling and dyadic refinement, the discrete GRG primitives converge monotonically as b increases.

Radial field with boundary, hole guard and trajectories Center-seeded field with trajectories to the boundary, auxiliary seed for a hole, and example paths. s (seed) aux seed ∂Ω hole guard
Center-seeded field with trajectories to the boundary, auxiliary seed for a hole, and example paths.

Practical Defaults

  • Field coefficients. a(x)=1 near the seed, tapering to maintain stable speed; b(x) small and smooth to avoid spirals.
  • Adaptive control. Increase samples where the curvature proxy exceeds a threshold; cap Δu.
  • Refinement (dyadic). Use the dyadic level b ∈ ℕ≥0 with N = 2^b samples and step Δu = 2^{-b}; increasing b yields monotone refinement.
  • Reconstruction. Quarter symmetry reconstructs the full figure; multi-seed blending covers non-simply connected domains.
Overview

Overview

GRG Overview

Generative Radial Geometry studies finite radial constructions and their measurement kernels. The current two-dimensional reference uses equal binary sectors, internal chord triangles, and external tangent triangles.

Current Reference

Definitions, kernel formulas, bounds, and exact-arithmetic error proofs.

Research History

Original Projection, Cartesian trapezoidal methods, and earlier Turns formulations retain their own definitions and scope. They are not interchangeable specifications.

Evidence

Geometric accuracy, floating-point stability, computational performance, and originality are separate questions. See Results and the benchmark protocol for their status.

GRG Original Projection

GRG Original Projection

Earlier Formulation

This page preserves a historical GRG method. Its notation and formulas belong to that version. The current binary radial formulation and verification are available here. Historical computational claims require separate validation.

Original Projection — historical method This page presents the original GRG idea: the circle is generated from the inside. The circle area is computed on a quadrant by projecting heights onto the x-axis and using midpoint rectangles (discrete upper); a matching right-endpoint rule gives a discrete lower, forming an under/over corridor that shrinks as $N$ increases (the midpoint rule is second-order). The perimeter is obtained from the sum of equal-angle chords on the quadrant (then $\times 4$). This method remains alongside other GRG blocks and is kept intact for comparison and continued study. See also Core Turns and Δx Trapezoidal.

Original Projection (Quarter → ×4)

Work on the quarter $x\in[0,r]$ with the profile $y(x)=\sqrt{r^{2}-x^{2}}$. Use uniform $\Delta x=r/N$. Midpoint rectangles provide a discrete upper estimate; right-endpoint rectangles provide a discrete lower. Multiply by $4$ to obtain full-circle values.

Setup

$$ \Delta x \;=\; \frac{r}{N}, \qquad x_k \;=\; k\,\Delta x \ (k=0,\dots,N), \qquad x_{k+\tfrac12} \;=\; \Bigl(k+\tfrac12\Bigr)\Delta x, \qquad y(x) \;=\; \sqrt{\,r^{2}-x^{2}\,}. $$

Quarter area: discrete corridor (upper vs lower)

Midpoint (discrete upper, second-order)

$$ A_{\text{quarter}}^{\text{mid}}(N) \;=\; \sum_{k=0}^{N-1} \Delta x \; y\!\bigl(x_{k+\tfrac12}\bigr), \qquad A_{\text{circle}}^{\text{mid}}(N) \;=\; 4\,A_{\text{quarter}}^{\text{mid}}(N). $$

Right-endpoint (discrete lower)

$$ A_{\text{quarter}}^{\text{right}}(N) \;=\; \sum_{k=0}^{N-1} \Delta x \; y\!\bigl(x_{k+1}\bigr), \qquad A_{\text{circle}}^{\text{right}}(N) \;=\; 4\,A_{\text{quarter}}^{\text{right}}(N). $$

Corridor: $A_{\text{circle}}^{\text{right}}(N) \;\le\; \text{Area(circle)} \;\le\; A_{\text{circle}}^{\text{mid}}(N)$, with gap shrinking like $O(N^{-2})$ under smoothness.

Perimeter via equal-angle chords (quadrant → ×4)

Partition the quarter by equal angular steps in turns. Let the quarter have $Q$ equal-angle chords of step $\Delta u=\tfrac{1}{4Q}$. The sum over the quarter is

$$ C_{\text{quarter}}(Q) \;=\; 2\,Q\,r\,s\!\Bigl(\tfrac{\Delta u}{2}\Bigr) \;=\; 2\,Q\,r\,s\!\Bigl(\tfrac{1}{8Q}\Bigr), \qquad C_{\text{circle}} \;\approx\; 4\,C_{\text{quarter}}(Q) \;=\; 8\,Q\,r\,s\!\Bigl(\tfrac{1}{8Q}\Bigr). $$

Equivalently, with $N=4Q$ steps on the full turn: $C_{N}=2\,N\,r\,s\!\bigl(\tfrac{1}{2N}\bigr)$ (see Core Turns).

Figure — Quarter with midpoint and right rectangles

Quarter: midpoint (upper) vs right-endpoint (lower) Quarter arc of a circle (r=100) with Δx strips. Green: midpoint rectangles; orange: right-endpoint rectangles. Both lie inside the quarter; midpoint sits above right. x y midpoint (upper) right (lower)
Midpoint vs right-endpoint corridor on the quarter.

Why use this block

  • Explicit under/over corridor with midpoint/right rules on the quarter.
  • Second-order convergence for the midpoint rule ($O(N^{-2})$) with simple operations.
  • No angle handling; complements chord-based perimeter via equal-angle steps.

How to use

  • Choose $N$ (e.g., $200$–$1000$) for the quarter; multiply results by $4$ for the circle.
  • Step-doubling check: compare $N$ vs $2N$ to estimate the residual gap between upper and lower.
  • Pair with equal-angle chords on the quarter to report both area and perimeter consistently.

See also

Core Turns · Δx Trapezoidal · Applications / Circle

Original Projection

The original GRG area method: sample the quarter along the x-axis with uniform step $\Delta x=r/N$. At each sample, project the height on the circle graph $y(x)=\sqrt{r^{2}-x^{2}}$ and build discrete rectangles. Two rules create a shrinking corridor: midpoint rectangles (upper) and right-endpoint rectangles (lower). Multiply by $4$ for the full circle. No angles are required in this block.

Setup

$$ \Delta x \;=\; \frac{r}{N}, \qquad x_k \;=\; k\,\Delta x \quad (k=0,\dots,N), \qquad x_{k+\tfrac12} \;=\; \Bigl(k+\tfrac12\Bigr)\Delta x, \qquad y(x) \;=\; \sqrt{\,r^{2}-x^{2}\,}. $$

Quarter-area sums (→ ×4 for the circle)

$$ \textbf{Midpoint (upper):}\quad A_{\mathrm{mid}}^{\mathrm{q}}(N) \;=\; \sum_{k=0}^{N-1} y\!\bigl(x_{k+\tfrac12}\bigr)\,\Delta x $$
$$ \textbf{Right-end (lower):}\quad A_{\mathrm{right}}^{\mathrm{q}}(N) \;=\; \sum_{k=0}^{N-1} y\!\bigl(x_{k+1}\bigr)\,\Delta x $$
$$ \textbf{Circle totals:}\quad A_{\mathrm{mid}} \;=\; 4\,A_{\mathrm{mid}}^{\mathrm{q}}, \qquad A_{\mathrm{right}} \;=\; 4\,A_{\mathrm{right}}^{\mathrm{q}} $$
$$ \textbf{Corridor:}\quad A_{\mathrm{right}} \;\le\; A_{\mathrm{true}} \;\le\; A_{\mathrm{mid}} $$
$$ \textbf{Estimate (optional):}\quad A_{\mathrm{avg}} \;=\; \frac{A_{\mathrm{mid}}+A_{\mathrm{right}}}{2} \quad\text{(use half-gap as an error bar).} $$

Notes

  • $y(x)=\sqrt{r^{2}-x^{2}}$ is concave on $[0,r]$; hence midpoint gives an upper estimate and right-endpoint a lower one.
  • Increase $N$ to tighten the corridor; publish $A_{\mathrm{avg}}$ with its half-gap if you need an error bar.

See also: Δx Trapezoidal · Core Turns · Applications / Circle

Circumference & Sphere

Circumference & Sphere

Why “Certainly on the Circumference” (and on the Sphere)

Category: Debate

Most ways to “draw a circle” are approximations: they split the arc into chords, use regular polygons, or compute points from arc lengths (where \(\pi\) eventually appears). In GRG we do the opposite: we do not measure the arc, we enforce the distance. We fix a radius \(R\) and choose a geometric direction; the generated point is the endpoint of a segment of length \(R\) from a center \(O\). By definition, that point is at distance \(R\) from \(O\): therefore it is certainly on the circumference. No approximation, no \(\pi\) needed.

What we prove, in plain words

In 2D, every point built as \((x,y)=(R\cos\varphi,\,R\sin\varphi)\) identically satisfies

$$x^{2}+y^{2}=R^{2}.$$

It does not “tend” to the circle: it is on the circle.

In 3D, take a unit direction \(\mathbf{u}\) and set \(\mathbf{P}=R\,\mathbf{u}\). By construction \(\|\mathbf{P}\|=R\): the point is certainly on the sphere

$$x^{2}+y^{2}+z^{2}=R^{2}.$$

The heart of it is: distance first. If the distance from the center is enforced to be \(R\), every generated point belongs to the locus of points at distance \(R\). That is the very definition of circle/sphere.

“Then what are sine and cosine for?”

Here they are not approximate numbers: they are ratios of sides in similar right triangles. They serve as a language to express the Cartesian projections of the chosen direction: \[ x=R\cos\varphi,\qquad y=R\sin\varphi. \] Even if, in a digital controller, you compute them with finite precision, the geometric law you implement remains exact: the commanded points are those of the ideal object defined by \(x^{2}+y^{2}=R^{2}\) (or \(x^{2}+y^{2}+z^{2}=R^{2}\)).

If you wish to avoid direct computation of \(\sin\) and \(\cos\) altogether, there is an equivalent incremental form: rotate the vector \((x,y)\) by multiplying it with an orthogonal matrix. Orthogonal matrices preserve the norm: every step remains on the circle by construction.

“What about error accumulating at each step?”

Distinguish carefully:

  • Geometric law (GRG): exact. The norm is preserved at each step (rotations), therefore there is no geometric drift.
  • Physical implementation: it can introduce errors (backlash, compliance, friction, encoder resolution, numerical rounding). These are real-world errors, not flaws of the idea.

In digital control, the small loss/gain of norm due to rounding is removed with three simple techniques:

\(\text{(i) adequate precision}\quad \text{(ii) occasional renormalization enforcing }x^{2}+y^{2}=R^{2}\quad \text{(iii) periodic recomputation from }\varphi.\)

In all cases, the reference object remains the perfect circle/sphere.

“But the machine moves in straight lines, not curves!”

What you often see is the visualization (a polyline between samples). It is not the definition of the path. The GRG law commands instantaneous poses that lie on the circle/sphere. The drive’s interpolation between samples may be linear or spline-based: that is a control choice. If you want a “fuller” trajectory, increase the sampling rate or use an arc interpolator; in any case, the nodes remain exactly on the geometric locus.

Why \(\pi\) is not needed

\(\pi\) is needed when you measure the arc (lengths, areas, times along an arc). In GRG we do not measure the arc: we enforce the distance from the center and the direction. Therefore \(\pi\) is irrelevant to guarantee that points lie on a circle/sphere. You may use \(\pi\) later if you care about arc length, but not to construct the point.

Natural extension to the sphere (including the \(Z\) axis)

What holds in \(X\!-\!Y\) holds in \(X\!-\!Y\!-\!Z\). A unit direction in 3D (obtained via Cartesian projections, without “walking on the surface”) is scaled by \(R\). The result has norm \(R\). That is all. Azimuth–elevation formulas are just a convenient way to obtain \((x,y,z)\). The key property remains the preservation of distance.

Practical consequences for robotics, CNC, and 3D printers

If your machine moves on Cartesian axes and you command \((x,y)\) (or \((x,y,z)\)) according to the GRG law, the ideal path is exactly the circumference (or the sphere).

Any deviations you might measure on a real part come from mechanics and numerics: backlash, compliance, calibration, friction, quantization. Improve the machine, the controller, and the metrology, and you approach the ideal without changing the GRG law.

Typical objections, short answers

  • “Without \(\pi\) I cannot make a circle.” False: a circle is the locus of points at distance \(R\) from a center. No \(\pi\) is needed to guarantee that distance.
  • “Sine and cosine are approximate numeric functions.” In geometry they are exact ratios; in control you can use incremental orthogonal rotations that preserve the norm by construction.
  • “Steps accumulate error.” Rotations do not change the norm. In digital control, rounding errors are contained with renormalization or periodic recomputation.
  • “Points on the sphere are not uniformly distributed.” Different topic: sampling uniformity does not affect the fact that each point lies on the sphere.

The key takeaway

GRG does not approximate the circle or the sphere: it imposes them. With a radius \(R\) fixed and a direction chosen, the distance from the center is \(R\) by construction: therefore the point is certainly on the circumference (or on the sphere). In a Cartesian system, the geometric precision of the trajectory does not depend on \(\pi\) nor on arc approximations: it depends only on your machine and numerics.

Certainly on the Circumference — and the Sphere

Category: Rigor of GRG

GRG conventions. Angles in turns \(u\in[0,1)\) with step \(\Delta u=1/N\). Radial computations (perimeter & radial area) on the full turn; Cartesian area on the quarter (×4). Trig appears only as classical cross-check.

Certainly on the Circumference

  • Circle \(R\): \(x^2+y^2=R^2\).

  • Direction (turns): \(u\in[0,1)\) → radial pose (no arc lengths).

  • GRG border: certificate \(x(u)^2+y(u)^2=R^2\) (Pythagoras on quarter or tangent intersections on turn).

Classical (optional): with \(\varphi=2\pi u\): \(x=R\cos\varphi,\ y=R\sin\varphi\).

Theorem

Every GRG point \(P(u)\) satisfies \(x(u)^2+y(u)^2=R^2\) ⇒ it lies on the circle.

Fig. — Radial construction / orthogonal update preserves norm.

Certainly on the Sphere

  • Sphere \(R\): \(x^2+y^2+z^2=R^2\).

  • Unit direction: build \(\mathbf{u}\) constructively; \(\mathbf{P}=R\,\mathbf{u}\) certifies \(\|\mathbf{P}\|=R\).

Classical (optional): \((\varphi,\theta)\) → \(x=R\cos\varphi\,\sin\theta,\ y=R\sin\varphi\,\sin\theta,\ z=R\cos\theta\).

Theorem

\(\|R\,\mathbf{u}\|^2=R^2\|\mathbf{u}\|^2=R^2\) ⇒ point on the sphere.

Fig. — 3D orthogonal rotations preserve the norm.

GRG · Operative summary (turns & quarter)

Params. Pick \(N\). Full turn: \(u_j=j/N\) with \(\Delta u=1/N\). Generate \(P_j=(x_j,y_j)\) with \(x_j^2+y_j^2=R^2\). Quarter: \(\Delta x=R/N\).

Perimeter — full turn (radial)

Inscribed (lower)

\[ L_{\mathrm{in}}(N)=\sum_{j=0}^{N-1}\sqrt{(x_{j+1}-x_j)^2+(y_{j+1}-y_j)^2}. \]

Circumscribed (upper) — tangents at \(P_j\): \(x\,x_j+y\,y_j=R^2\); intersections \(T_j\) with \(\det_j=x_jy_{j+1}-y_jx_{j+1}\), \(T_j^x=\frac{R^2(y_{j+1}-y_j)}{\det_j}\), \(T_j^y=\frac{R^2(x_j-x_{j+1})}{\det_j}\).

\[ L_{\mathrm{out}}(N)=\sum_{j=0}^{N-1}\lVert T_{j+1}-T_j\rVert. \]

Radial area — full turn

\[ A_{\mathrm{in}}(N)=\tfrac12\sum_{j=0}^{N-1}(x_jy_{j+1}-x_{j+1}y_j),\qquad A_{\mathrm{out}}(N)=\tfrac12\sum_{j=0}^{N-1}(T_j^xT_{j+1}^y-T_{j+1}^xT_j^y). \]

Cartesian area — quarter (Δx) → ×4

Abscissas \(x_i=i\Delta x\), \(\Delta x=R/N\); heights \(y_i=\sqrt{R^2-(i\Delta x)^2}\).

\[ A_{\mathrm{right,q}}=\sum_{i=1}^{N}y_i\Delta x,\quad A_{\mathrm{mid,q}}=\sum_{i=1}^{N}y\!\left(x_{i-\tfrac12}\right)\Delta x,\quad A_{\mathrm{trap,q}}=\sum_{i=1}^{N}\tfrac{y_{i-1}+y_i}{2}\Delta x. \]

Bounds: \(A_{\mathrm{right,q}}\le A_{\mathrm{quarter}}\le A_{\mathrm{mid,q}}\); disk via ×4. π remains outside the computation; trig only as cross-check.

The Circle

The Circle

The Circle — Radius-First (Three Methods)

Principle

The circle is generated from the center. The radius is the only primitive; the perimeter is the path traced by the radius, and the area is the surface swept by its rotation. We work in turns (fractions of a full rotation), not in sexagesimal degrees.

Angles in turns & radial node

Let \(u\in[0,1)\) be an angle in turns (period \(1\)). We use period-1 trigonometric functions \(c(u), s(u)\) with \(c(0)=1, s(0)=0\) and \((c(\tfrac14), s(\tfrac14))=(0,1)\). The exact radial node is:

$$ P(u)=\bigl(r\,c(u),\; r\,s(u)\bigr). $$

Bit-friendly turns

Use fixed-point angles \(u = k/2^B\). Picking \(N\) that divides \(2^B\) makes equal-angle steps exact (\(\Delta u=1/N\)).
Optional: LUT + interpolation for \(c(u), s(u)\).

Method scope & scaling (quick key)

  • Radial Sectors (equal-angle, full turn) → no ×4.
  • Δx Trapezoidal (quarter tiling) → ×4 to get the full circle.
  • Original Projection (midpoint/right on a quadrant) → ×4.
Perimeter (equal-angle chords) is naturally full-turn (no ×4) and shared across methods.

Perimeter (shared across methods)

Equal-angle chords (full turn)

$$ \boxed{\, C_N \;=\; 2\,N\,r\,s\!\Bigl(\tfrac{1}{2N}\Bigr) \,} $$
Inscribed regular \(N\)-gon in turns. Monotone ↑ in \(N\). No quadrant scaling.
Perimeter by equal-angle chords (Δu = 1/N turns) Cartesian axes; circle r=100; sample radii and dashed chords (example N=12).
Perimeter by equal-angle chords: \(C_N = 2N r\, s(1/(2N))\) (full turn).

Area — three interchangeable builds (pick one)

Δx Trapezoidal (quarter → ×4)

$$ \Delta x=\frac{r}{N},\quad x_k=k\,\Delta x,\quad y_k=\sqrt{r^2-x_k^2} $$
$$ A_{\text{quarter}}(N)=\frac{\Delta x}{2}\!\left(y_0+2\sum_{k=1}^{N-1}y_k+y_N\right),\qquad A(N)=4\,A_{\text{quarter}}(N). $$
Under-approx on a concave arc; smooth convergence \(O(N^{-2})\).

Original Projection (midpoint/right corridor, quarter → ×4)

$$ \Delta x=\frac{r}{N},\ \ x_{k+\tfrac12}=(k+\tfrac12)\Delta x,\ \ y(x)=\sqrt{r^2-x^2} $$
$$ A^{\text{quad}}_{\text{mid}}=\sum_{k=0}^{N-1} y\!\bigl(x_{k+\tfrac12}\bigr)\,\Delta x,\qquad A^{\text{quad}}_{\text{right}}=\sum_{k=0}^{N-1} y\!\bigl(x_{k+1}\bigr)\,\Delta x, $$
$$ A_{\text{mid/right}}=4\,A^{\text{quad}}_{\text{mid/right}},\qquad A_{\text{right}}\le A_{\text{true}}\le A_{\text{mid}}. $$
Midpoint is second-order; corridor shrinks with \(N\).

Radial Sectors (equal-angle sectors, full turn)

$$ \boxed{\, A_N = \tfrac{N}{2}\,r^2\,s\!\Bigl(\tfrac{1}{N}\Bigr) \,} $$
Sum of \(N\) congruent isosceles triangles with included angle \(1/N\) turns (inscribed lower; monotone ↑ in \(N\)). No ×4.
Quarter area by Δx trapezoids (under the arc) Cartesian x–y axes; quarter arc with Δx verticals and interior trapezoids. Δx trapezoids
Quarter area by Δx trapezoids; multiply by 4 for the full circle.
Radial sector tiling (equal-angle sectors, full turn) Two example sectors on a unit circle (illustrative N=8).
Radial sectors: \(A_N=\tfrac{N}{2}r^2\,s(1/N)\) (full turn, no ×4).

Practical notes

  • Pick \(N\). If angles are \(k/2^B\) turns, choose \(N\mid 2^B\). Use step-doubling (\(N\!\to\!2N\)) to gauge the gap.
  • Monotonicity. \(C_N\) ↑ with \(N\). Δx trapezoids ↑ toward area (under-tiling). Original corridor shrinks. Radial sectors ↑ (inscribed).
  • All in turns. No degrees, no π in the constructive steps.