A multivariate full conformal prediction (FullCP) procedure accepts or rejects each candidate response t through an exact conformal comparison, and the collection of accepted candidates is the prediction region studied here. Under the standard exchangeability framework, with i.i.d. sampling as a canonical special case, conformal prediction supplies its usual finite-sample marginal validity (Shafer & Vovk, 2008; Vovk et al., 2005). This paper addresses a different layer of the problem: the geometry of that exact accepted set and how to represent it faithfully. The validity guarantee is inherited background; the contribution is score-specific exact-set geometry and deterministic reconstruction error control.

The present analysis uses the response-only special case of exchangeable conformal prediction: the observed examples are Y1,…,Ym, with no covariates and no fitted predictive model. This restriction isolates the exact-set geometry induced by the score itself from the additional candidate dependence introduced by predictive-model fitting. Shafer & Vovk (2008) discuss the corresponding special case in which prediction is based on past examples without informative features of the new object. By contrast, several recent multivariate conformal methods formulate prediction for regression pairs (X,Y) and/or incorporate fitted or pre-trained predictors (Diquigiovanni et al., 2022; Gan & Liu, 2025; Johnstone & Ndiaye, 2025). In refitted regression-based FullCP, the trial response can also enter the fitted predictor through the augmented training data (Lei et al., 2018), so the resulting score comparisons can depend on the trial value through model fitting as well as through direct score evaluation. Accordingly, Theorem 1 gives an exact characterization of the response-only score-induced setting studied here; it does not establish common-center geometry for refitted regression-based FullCP.

For example, at a fixed forecast context one may have exchangeable historical low-dimensional multivariate outcomes, with no predictor refit as the trial response varies. The empirical energy-score FullCP region can then be used directly for reporting or downstream candidate-level decisions.

That distinction matters in more than one response dimension. An exact FullCP region can be nonconvex, while a finite collection of directional boundary searches leaves uncertainty between sampled directions. Replacing the resulting boundary by an ellipse or convex hull may be convenient, but if the exact region is nonconvex then convexification can enlarge it and add candidates rejected by the exact conformal comparison. Such an enlargement is not claimed here to destroy marginal coverage validity; rather, it is no longer a faithful representation of the exact accepted set and can change its size, shape, and candidate-level decisions.

For the empirical energy-form pairwise score considered below, the leave-one-out algebra supplies a stronger structure than candidate-score convexity. Each training-versus-candidate comparison becomes exactly a pairwise-dissimilarity sublevel condition. Under the stated diagonal and attainment assumptions, all comparison regions contain a common Fréchet-type minimizer; when those comparison regions are convex, the nontrivial exact FullCP vote aggregate is star-shaped about that same point. Star-shapedness is the useful middle ground: the aggregate need not be convex, but every ray from the common center has a single radial exit, so the whole region can be represented ray by ray.

Power distances ρβ(x,y)=‖x−y‖β provide the main concrete family. Their deterministic common-center geometry holds for β≥1. In one dimension at β=1, the score has the usual energy representation of the continuous ranked probability score (CRPS) (Gneiting & Raftery, 2007); for the fixed-bin empirical construction studied by Toccaceli (2026, ​Sections 8.1 and 8.3), the score geometry proves that every nontrivial exact FullCP set is a nonempty closed interval. This CRPS result is a recognizable specialization of the multivariate mechanism, not a new scoring-rule theorem.

The radial representation becomes certifiable on the regular range 1<β<2. Along each direction, the exact FullCP radius is an order statistic of the component comparison-set radii. Explicit score-derived Lipschitz bounds on those radii then turn finitely many certified ray searches into global lower and upper radial envelopes. The envelope width combines the one-dimensional search tolerance with a directional covering term, and the same-center radial sandwich yields deterministic Hausdorff error bounds. Thus finite-direction reconstruction is accompanied by a global approximation certificate rather than only a plotted collection of rays.

The paper makes three contributions. First, it gives the exact score-difference reduction, common-center comparison geometry, and star-shaped aggregation, including the exact univariate empirical-CRPS interval result. Second, it identifies the radial order-statistic interface and, for 1<β<2, explicit data-checkable radial regularity that makes finite-direction certification possible. Third, it derives a focuswise radial certificate that is deterministically no looser than a generic center-based benchmark and assesses the resulting certificate tightness and nonconvexity diagnostics in a frozen staged two-dimensional study. The analytic example in Section 4 shows explicitly why preserving nonconvexity can matter before the certification machinery is introduced.

The scope is deliberately limited. The geometric statements are specific to the score structure and assumptions developed here; no generic theorem for arbitrary FullCP regions is claimed. The paper does not add a conformal coverage theorem, and response-dependent adaptive partition questions remain outside scope. Directional covering restricts practical certification to low-dimensional multivariate outputs, and no runtime superiority is asserted. The numerical study is descriptive for its frozen designs rather than an estimate of how often nonconvexity occurs in a population of FullCP problems.

The standard variogram score used in multivariate forecast evaluation is not of the empirical observation-pair energy-form V-statistic form assumed in Theorem 1, so Theorem 1 does not directly cover it. This scope statement does not rule out analogous geometry for other variogram constructions.

Root finding and bisection for conformal sets predate this work (Ndiaye & Takeuchi, 2023), as does multi-output reconstruction from finitely many search directions (Johnstone & Ndiaye, 2025). In the latter framework, star-shapedness can be imposed as a structural condition and sampled boundary points can be fit by convex objects; here the common star center is instead derived from the score before directional approximation. Classical radius- R supporting-ball geometry (Balashov & Golubev, 2012) and McShane–Whitney Lipschitz-extension formulas (Arnau et al., 2023) are used later as inherited tools after the score-specific bounds are established. Other certified FullCP approximations use different error metrics and scopes (Razafindrakoto et al., 2026b, 2026a). A companion Lean 4 development checks selected discrete and algebraic components of the present results; the multivariate regularity, adaptive-partition, and coverage questions are outside that formalization.

Section 2 develops the exact score-induced FullCP geometry. Section 3 specializes to power distances, isolates the empirical-CRPS interval result, and gives the exact radial order-statistic representation. Section 4 gives the analytic nonconvex witness before Section 5 develops certified reconstruction for 1<β<2. Section 6 reports the staged numerical assessment, Section 7 records computational and endpoint limitations, and Section 8 synthesizes the scope. Appendices A–D contain the β=1 segment pathology, a three-point robustness witness, the detailed technical regularity calculations, and numerical protocol and reproducibility details.

The setup can be read as a voting rule. For each candidate t, training observation j contributes one comparison vote when its leave-one-out score is at least the candidate score; the exact FullCP region consists of candidates receiving enough such votes. The notation below makes that rule explicit before exploiting the score algebra.

Let m≥1, d≥1, and let y1,…,ym∈ℝd. For a function ρ:ℝd×ℝd→ℝ and an empirical distribution F̂=1m∑r=1mδzr, define the energy-form pairwise score Sρ(F̂,y)=1m∑r=1mρ(zr,y)−12m2∑r,s=1mρ(zr,zs). The second term uses the empirical V-statistic normalization 1/(2m2); no off-diagonal normalization is used.

For a candidate t∈ℝd, augment the observed responses by xi(t)=yi,i=1,…,m,xm+1(t)=t. For r=1,…,m+1, let F̂−r,t=1m∑s=1s≠rm+1δxs(t),αr(t)=Sρ(F̂−r,t,xr(t)), and write α(t)=αm+1(t). For j=1,…,m, define the comparison region Ij={t:αj(t)≥α(t)}. Thus score ties are accepted as comparison votes. Let N(t)=∑j=1m𝟏{t∈Ij},p(t)=1+N(t)m+1, and use throughout the strict conformal convention \varepsilon\}.]]> Γε={t:p(t)>ε}. For 0≤ε<1, set k=⌊(m+1)ε⌋. Thus the exact endpoint regimes are k=0(0≤ε<1m+1)⇒Γε=ℝd, while 1≤k≤m is the nontrivial vote-threshold regime; in that regime the integer-valued vote count gives the exact equivalence Γε={t:N(t)≥k}. At ε=1, the strict convention gives Γ1=⌀. These endpoint, tie, and threshold conventions will not be repeated inside later formal statements.

Candidate-score convexity is not sufficient. To see why the score-specific structure below is needed, temporarily replace the energy-form score by the generic reference-sample-symmetric nonconformity rule A(z;Z)=1|Z|∑w∈Z|z−w2|=1|Z|∑w∈Zℓ(z,w),ℓ(u,v)=|u−v2|. The rule is symmetric in the reference multiset Z, and ℓ(⋅,v) is convex for every fixed v.

Take m=2 and y1=y2=2. For a candidate t, its score is α̃(t)=A(t;{2,2})=|t−4|, which is convex, whereas the two training scores coincide: α̃1(t)=α̃2(t)=A(2;{2,t})=1+12|2−t2|. A direct piecewise calculation gives α̃j(t)≥α̃(t)⇔t≤−4ort≥2. Thus the corresponding full-conformal p-value is p̃(t)={1,t∈(−∞,−4]∪[2,∞),1/3,−4<t<2. At ε=1/2, using the strict convention \varepsilon]]> p̃(t)>ε, Γ̃1/2=(−∞,−4]∪[2,∞), which is disconnected. Equality holds at t=−4 and t=2, so both endpoints are included because comparison votes use ≥.

This example is not asserted to arise from a proper distributional scoring rule. Its role is narrower: candidate-score convexity is not sufficient for connected FullCP regions, even within the empirical-average form with a loss convex in its candidate argument.

The following result isolates the stronger algebraic and geometric structure that is sufficient for the energy-form score.

Before stating the result, define dj=∑i≠jρ(yj,yi),F(t)=∑i=1mρ(t,yi).

Theorem 1 (Exact common-center geometry for empirical energy-form FullCP). Assume ρ is symmetric and has constant diagonal ρ(x,x)=q. Then:

Proof. Fix j and put Bj=∑r,s=1r,s≠jmρ(yr,ys),Tj(t)=∑i≠jρ(t,yi). Leaving t out of the augmented sample gives α(t)=1m{ρ(t,yj)+Tj(t)}−12m2{Bj+2dj+q}, whereas leaving yj out gives αj(t)=1m{ρ(t,yj)+dj}−12m2{Bj+2Tj(t)+q}. Subtracting yields part (i).

For part (ii), global minimality of c gives F(c)≤F(yj)=dj+q, while ρ(c,yj)≥q. Hence ∑i≠jρ(c,yi)=F(c)−ρ(c,yj)≤dj, so c∈Ij for all j.

For part (iii), let t∈Γε. Then t belongs to at least k comparison sets, and the same sets contain c. Convexity places every point of [c,t] in each of those k sets, so every point of the segment receives at least k votes. Thus [c,t]⊆Γε. ◻

The theorem separates three logically distinct ingredients. Symmetry and a constant diagonal suffice for the exact comparison identity. The diagonal lower bound and attainment of F give a common point of all comparison regions. Convexity of the individual comparison regions is needed only for the final segment argument. The preceding counterexample shows that convexity of the candidate score itself is not a substitute for this comparison-set structure. The resulting k-coverage set can be star-shaped without being convex; Section 4 gives an explicit analytic example.

Consider the power-distance dissimilarity ρβ(x,y)=‖x−y‖β. This family gives a concrete deterministic geometry and, on a narrower parameter range, the conventional energy-score interpretation.

Corollary 1 (Common-center star-shapedness for power-distance scores). For β≥1, define Fβ(t)=∑i=1m‖t−yi‖β,Gj(t)=∑i≠j‖t−yi‖β,dj=∑i≠j‖yj−yi‖β. Then Fβ attains a global minimum, each comparison set is Ij={t:Gj(t)≤dj}, and every Ij contains any chosen minimizer c of Fβ. Consequently, for every nontrivial threshold 1≤k≤m, the exact FullCP region is star-shaped about c.

Proof. For β≥1, x↦‖x−y‖β is convex for fixed y, while ρβ is symmetric, has diagonal 0, and is nonnegative. The function Fβ is continuous and coercive, hence attains a global minimum, and each Gj is convex. Theorem 1 applies with q=0. ◻

Finite exits and the m=1 degeneracy. When m=1, G1 is the empty sum and d1=0, so I1=ℝd and finite-radius reconstruction is degenerate. When m≥2, each Gj contains at least one coercive term and Ij={Gj≤dj} is bounded, so every ray from a common center has a finite exit.

For d=1 and β=1, this score is exactly the fixed-bin empirical CRPS score used by Toccaceli (2026, Section 8.1). The univariate case deserves a separate formal statement because it gives the exact interval property for that construction. Write the ordered observations as y(1)≤⋯≤y(m) and define the full empirical median set ℳ=arg⁡mintF1(t)=[y(⌈m/2⌉),y(⌊m/2⌋+1)].

For clarity, the augmented fixed-bin construction sets xi(t)=yi(i≤m),xm+1(t)=t,F̂−r,t=1m∑s=1s≠rm+1δxs(t). Its score is CRPS⁡(F̂−r,t,xr(t))=1m∑s=1s≠rm+1|xs(t)−xr(t)|−12m2∑s,ℓ=1s,ℓ≠rm+1|xs(t)−xℓ(t)|. This is exactly the energy-form score with ρ(x,y)=|x−y| and the same empirical V-statistic normalization.

Corollary 2 (Exact interval property for univariate empirical CRPS FullCP). Let d=1 and β=1. For every nontrivial threshold 1≤k≤m, the fixed-bin empirical CRPS FullCP region is a nonempty closed interval containing ℳ. For m≥2 the interval is compact; for m=1 it is ℝ.

Proof. Theorem 1 applies separately to every c∈ℳ, so every point of ℳ belongs to every comparison region. Each Ij is a closed interval because it is a closed convex sublevel set in ℝ. For 1≤k≤m, Γε=⋃S⊆{1,…,m}|S|=k⋂j∈SIj, so Γε is closed and is nonempty because every point of ℳ receives all m votes. It is star-shaped about every c∈ℳ by Theorem 1; a nonempty star-shaped subset of ℝ is an interval. For m≥2 the comparison sets are bounded and k≥1, hence the interval is bounded and therefore compact. The m=1 case is the degeneracy noted above. ◻

This proves the single-interval property reported for Toccaceli’s fixed-bin empirical-CRPS construction. The interval conclusion follows from the deterministic comparison geometry and does not by itself extend to proper distributional scoring rules in general.

Scoring-rule range. Under the conventional distributional energy-score interpretation, with 𝒫β={P:∫‖x‖βdP(x)<∞}, the power-distance energy score is strictly proper on the appropriate finite-moment class for 0<β<2 (Gneiting & Raftery, 2007). At β=2 it is proper but not strictly proper and its expected score depends only on the forecast mean; for 2]]> β>2 the Euclidean distance-power energy-form score is not proper. Thus strict propriety overlaps the deterministic star-shaped geometry on 1≤β<2. These are classical scoring-rule facts and do not enlarge the deterministic claims above.

Lower-boundary sharpness. The deterministic conclusion should not be extended below β=1. With m=3,β=12,(y1,y2,y3)=(−1,0,2), the unique Fréchet minimizer is c=0, yet at ε=0.8 the exact prediction set contains 0 and 1 while rejecting 0.5. Thus it is not star-shaped about c.

For the remainder of the radial development, assume the comparison sets are star-shaped about a common center c and define, for u∈𝕊d−1, Rj(u)=sup⁡{r≥0:c+ru∈Ij}. For a nontrivial threshold, define likewise Rε(u)=sup⁡{r≥0:c+ru∈Γε}, and let R(1)(u)≤⋯≤R(m)(u) be the ascending order statistics of the component exits.

Proposition 1 (Radial order-statistic representation). For 1≤k=⌊(m+1)ε⌋≤m, Rε(u)=R(m−k+1)(u)for every u∈𝕊d−1.

Proof. Along a fixed ray, star-shapedness makes the membership set of each Ij an initial interval with supremum Rj(u). The supremum radius at which at least k votes remain is therefore the kth-largest component exit, equivalently the (m−k+1)st ascending order statistic. ◻

The endpoint conventions remain those fixed in Section 2: k=0 gives Γε=ℝd and hence Rε≡+∞, while ε=1 gives the empty set. For power distances with m=1, I1=ℝd and R1≡+∞, so finite radial reconstruction begins at m≥2. At β=1, exact common-center star-shaped geometry remains valid but unconditional global radial Lipschitzness can fail; Appendix A gives an exact discontinuous example.

For this section, write Γ=Γε. A two-dimensional analytic example makes the distinction between star-shapedness and convexity explicit before the reconstruction machinery is introduced. Let d=2,m=2,β=32, with y1=(−1,0),y2=(1,0),ε=12. Then k=⌊3/2⌋=1. For j=1, the comparison objective contains only the distance to y2, and d1=‖y1−y2‖3/2=23/2. Therefore I1=B(y2,2),I2=B(y1,2), and the one-vote prediction region is exactly Γε=I1∪I2. The origin is a common center. With u(θ)=(cos⁡θ,sin⁡θ), the two ball exits are cos⁡θ+cos⁡2θ+3and−cos⁡θ+cos⁡2θ+3, so the union has exact radial function R(θ)=|cos⁡θ|+cos⁡2θ+3. This set is star-shaped about 0 but nonconvex: (−1,2),(1,2)∈Γε,(0,2)∉Γε.

Its convex hull is the stadium conv⁡Γε=[−e1,e1]+B(0,2). Direct geometry gives the exact radial convexification discrepancy ‖Rconv⁡Γ−RΓ‖∞=2−3, and the exact Hausdorff discrepancy dH(Γ,conv⁡Γ)=5−2. The latter is attained at the midpoint of the top or bottom flat portion of the stadium boundary, whose distance to either radius- 2 component ball is 5−2. These positive discrepancies are properties of this analytic witness; they do not imply that convex reconstruction is inappropriate in every conformal problem. In particular, because the convex hull is a superset of the exact accepted set, no claim is made that this enlargement itself destroys marginal coverage validity. The point is fidelity: convexification adds candidates rejected by the exact comparison and therefore changes the exact FullCP region.

For a same-scale comparison, the analytic union has area and equivalent-area radius AΓ=16π3+23,req=163+23π≈2.5369255259. Thus its radial gap 2−3 has normalized magnitude 2−3req≈0.1056196525. This analytic example is a stress-test and existence illustration: convexification can be materially nonfaithful, but it makes no representative-magnitude claim. In contrast, conditional on ROBUST_WITNESS, the frozen numerical study in Section 6 had normalized radial-gap median 0.00137, 90th percentile 0.00626, and maximum 0.01770. The frozen DGPs therefore had often-detectable but usually shallow nonconvexity. At these observed magnitudes convexification may be negligible for coarse visualization or coarse summaries, while exact boundary-near candidate membership can nevertheless change.

The nonconvexity is not an artifact of using only two observations. In the three-point example of Appendix B with m=3, β=3/2, and k=2, two symmetric endpoints are accepted while their midpoint is rejected, with nonzero comparison margins. The numerical margins are recorded only in the appendix. The computational implication is that the natural exact target is the star-shaped radial function, rather than a convex surrogate, which motivates the certified reconstruction developed next.

We now restrict to m≥2, 1<β<2, and a nontrivial threshold 1≤k≤m. The goal is to turn the exact radial representation into a global certificate from finitely many ray searches. Throughout, B(a,R)={x∈ℝd:‖x−a‖≤R} denotes the closed Euclidean ball.

The numerical study assesses certificate tightness and the detectability and radial magnitude of nonconvexity under frozen two-dimensional designs; it is not a coverage experiment and its detection frequencies are not estimates of population prevalence. The study was staged. The initial smoke run used 5 independent clouds per (m,data-generating process (DGP)) cell, 40 total, and expansion to 20 independent clouds per cell was allowed only after the seven advance-declared engineering and numerical gates listed in Appendix D passed; soft flags did not automatically block expansion. After those gates passed, the expanded-production β/qvote grid, seed mapping, and no-post-hoc rule were frozen before any expanded-run results were observed. The original five clouds in each cell were retained in the final 20, so the final 160 clouds include the original 40 smoke clouds.

The 20 clouds per (m,DGP) cell were the staged descriptive expansion size, not a target-precision sample size. For each m∈{10,25}, the final dataset contained 20 independent clouds from each of four designs—an isotropic Gaussian, a rotated anisotropic Gaussian, a symmetric Gaussian mixture, and a standardized bivariate t3 distribution.1 The fixed grid was β∈{1.1,1.5,1.9} and qvote∈{0.10,0.20,0.50}, with each cloud reused across the β/qvote settings. The threshold target determines k=max⁡{1,⌊qvotem+1/2⌋}, whose corresponding strict-conformal ε interval is [k/(m+1),(k+1)/(m+1)). Certificate summaries therefore operate at the 480 cloud– β level, whereas witness and radial summaries use repeated β/qvote configurations within cloud; the resulting 1,440 cloud– β– k rows are not independent observations. Values at 16,384 directions occur only for cases that triggered the prespecified fallback rule in Appendix D.

Across the 480 cloud– β analyses, the median ratio L̃/Lold was 0.966, with interquartile range (IQR) approximately [0.941,0.976], equivalent to a median tightening of 3.39%. Descriptively, the median ratio was 0.941 (IQR [0.932,0.955]) for m=10 and 0.975 (IQR [0.972,0.980]) for m=25; no inferential significance is attached to this difference. Consistent with the deterministic ordering in Proposition 2, L̃≤Lold in every one of the 480 analyses. Figure 3 summarizes the certificate comparison.

To fix the numerical scale, on the full validation mesh define Wmax(L,N)=maxu[R¯L,N(u)−R_L,N(u)]. On the 8,192-direction radial mesh, with radial values R̂i, define AΓ=π8192∑i=18192R̂i2,req=AΓ/π. For the operational target Wmax≤0.10req, neither certificate achieved the target at N≤128 in the expanded study. At N=256, the descriptive fraction was 128/1440=8.89% for Lold and 165/1440=11.46% for L̃. These fractions aggregate repeated cloud– β– k configurations and are not independent-observation estimates. The tilde envelope width was smaller in all 7,200 declared analysis-by- N comparisons, as expected from the deterministic certificate ordering.

For the frozen 4,096-direction diagnostic mesh, the radial brackets are R̂i±ei, and we write the finite-grid lower proxy as Lsec−=maxi≠jmax⁡{0,|R̂i−R̂j|−ei−ej}‖ui−uj‖. This is a finite-grid lower proxy, not the true Lipschitz constant.2 The median certificate-to-finite-grid-lower-proxy ratio was 108.8 (IQR [63.5,201.4]) for Lold and 105.9 (IQR [59.9,193.9]) for L̃. Substantial residual slack remains relative to the finite-grid lower proxy. The focuswise tightening therefore removes only part of the gap visible against this numerical lower diagnostic.

Figure 4 reports robust-witness detection rates over the prespecified grid. Operationally, the detector forms witness hulls from the common center and genuinely accepted radial lower endpoints, inspects hull edges joining nonadjacent angular samples, and searches convex combinations after any required inward endpoint backoff; ROBUST_WITNESS requires both endpoints to be robustly accepted and the convex-combination point to be robustly rejected under the stored aggregate-margin intervals, with a component/direct radial cross-check. SUGGESTIVE records a prescribed candidate that remains unresolved after the frozen refinement procedure, whereas NONE_DETECTED means only that this finite diagnostic found no witness; Appendix D gives the search details. Each fixed β/qvote cell contains exactly 160 independent clouds, and the cellwise detection rates range from 40.0% to 81.2%. At the independent-cloud level, 155/160 clouds yielded at least one robust witness somewhere among the nine prespecified β/qvote configurations. These observed detection frequencies do not estimate the population prevalence of nonconvexity.

For the reported magnitude, form the metric hull HM=conv⁡{c+R̂iui} and define the mesh-based radial convexification-gap approximation Δ̂rad=maxu[RHM(u)−R̂(u)]+, reported after normalization by req. Conditional on ROBUST_WITNESS, this normalized approximation had median 0.00137, 90th percentile 0.00626, and maximum 0.01770. These summaries use the final gap convention. At baseline, the fallback is triggered when 10^{-4}]]> max⁡{|Δ̂A−Δ̂8192|,|Δ̂B−Δ̂8192|}/req>10−4 for the two interlaced 4,096-direction submeshes; the nested 16,384-direction refinement is then automatic, and its normalized 8,192-versus-16,384 discrepancy is checked against the same threshold. Exactly three cases triggered this fallback and none remained unresolved; Appendix D gives the frozen rule. In this staged numerical study, nonconvexity was frequently detectable, while its normalized radial magnitude was typically small.

Post-review Monte Carlo uncertainty summaries from the frozen outputs. No simulation was rerun. A deterministic design-stratified cloud bootstrap uses only the frozen certificate and analysis CSVs: within each of the eight (m,DGP) strata, cloud IDs are sorted and 20 whole clouds are resampled with replacement, retaining all their β/qvote rows. With numpy.random.default_rng(20261001) and 20,000 replicates, the same cloud resample computes all requested summaries in each replicate and preserves the balanced 20-per-stratum design. The 2.5th and 97.5th percentiles for the median L̃/Lold are [0.96380,0.96909] overall, [0.93904,0.94371] for m=10, and [0.97457,0.97609] for m=25. Table 1 gives the detection-fraction limits. These quantify finite Monte Carlo resolution of the fixed staged design; they are not population prevalence or coverage confidence intervals and imply no formal sample-size sufficiency or power claim.

Theorem 2 and Corollary 3 reduce global certification to certified values on a finite directional set and a sphere-covering radius. For chord-metric covering of 𝕊d−1, the fixed-dimension small-radius law is Nd(h)=Θd(h−(d−1)),h↓0, consistent with standard spherical covering and packing exponents (Bhattacharya et al., 2022). If a target radial-envelope width τ satisfies \delta]]> τ>δ and 0]]> L̃>0, the deterministic bound δ+2L̃h≤τ is guaranteed whenever h≤τ−δ2L̃. Consequently, for fixed d and fine resolution, Nd=Θd[(2L̃τ−δ)d−1]. If L̃=0, the certified width is already at most δ. This fixed- d notation is not a dimension-growth statement; for fixed small h, log⁡Nd(h)=(d−1)log⁡(1/h)+O(d). The directional burden is therefore a substantive limitation and confines the intended certification regime to low-dimensional multivariate outputs.

The endpoint β=1 separates exact geometry from unconditional certification. Corollary 1 and Theorem 1 still give common-center star-shaped geometry, but the global radial function need not be continuous or Lipschitz; Appendix A gives the exact segment example. Proposition 2 is therefore stated unconditionally only for 1<β<2. At the separate degeneracy m=1, I1=ℝd, so finite power-distance radial reconstruction begins at m≥2.

The upper restriction β<2 also keeps the certified development within the strictly proper range of the conventional energy score, but the deterministic and scoring-rule roles should not be conflated. Convex comparison geometry holds for every β≥1, including β≥2; strict propriety has the separate boundaries summarized in Section 3. Adaptive response-dependent partition issues remain outside scope, and the paper makes no new coverage theorem or runtime-superiority claim.

The numerical assessment is limited to d=2, four prespecified DGPs, and a finite β/qvote grid. NONE_DETECTED does not certify convexity, the observed detection frequencies are not prevalence estimates, and the finite-grid secant statistic is only a lower proxy for the true Lipschitz constant.

The statistical distinction driving the paper is between validity of an exact conformal set and faithful representation of that set. Standard conformal theory supplies the former. For the energy-form score studied here, the paper supplies the latter structure: exact leave-one-out comparisons become pairwise-dissimilarity sublevel sets with a common center, and convexity of those comparison sets makes the nontrivial vote aggregate star-shaped without forcing it to be convex. The univariate empirical-CRPS interval result is one specialization of this mechanism, while the analytic example shows why convexification can change the exact accepted set even when inherited coverage validity is not in question.

For 1<β<2, the star-shaped geometry also supports deterministic reconstruction guarantees. Explicit radial regularity turns finite directional brackets into global radial and Hausdorff certificates, with the formal dominance L̃≤Lold preserved. In the frozen numerical study the median certificate tightening was approximately 3.4%: systematic but modest, with both certificates remaining roughly two orders above the finite-grid lower proxy and hence substantial residual certificate slack. Robust nonconvexity witnesses were often detected but typically shallow in normalized radial magnitude. These observations are descriptive for the frozen staged two-dimensional designs.

The analytic example is a stress-test and existence illustration of materially nonfaithful convexification, not a representative-magnitude claim. At the smaller observed numerical gaps, convexification may be negligible for coarse visualization or coarse summaries, while exact boundary-near candidate membership can still change. The results remain score-specific, low-dimensional, and separate from adaptive-partition coverage questions or runtime claims.

These results concern response-only energy-form FullCP. In regression-based FullCP with candidate-dependent refitting, the trial response can alter the fitted predictor and hence the score comparisons, so the common-center reduction proved here does not automatically extend to that setting. Characterizing analogous geometry under fitted-model constructions is outside the present scope.

Let d≥2, m=3, and y1=−e1,y2=0,y3=e1,β=1. Take k=3, which is equivalent under the strict convention to 34≤ε<1. Since all three comparison votes are required, Γε=I1∩I2∩I3. The middle comparison region is I2={t:‖t+e1‖+‖t−e1‖≤2}. By the triangle inequality, ‖t+e1‖+‖t−e1‖≥‖2e1‖=2, and equality holds precisely on the segment joining the two foci. Hence I2=[−e1,e1]. For t=se1 with −1≤s≤1, ‖t‖+‖t−e1‖≤3,‖t‖+‖t+e1‖≤3, so the whole segment also lies in I1 and I3. Therefore Γε=[−e1,e1].

The Fréchet objective is minimized at c=0. About this common center, Rε(e1)=Rε(−e1)=1, whereas for every direction u not collinear with e1 the ray 0\}]]> {ru:r>0} leaves the segment immediately, so Rε(u)=0. Thus the radial function is discontinuous on 𝕊d−1 when d≥2. Exact common-center star-shaped geometry at β=1 therefore does not imply an unconditional global Lipschitz radial certificate.

Consider d=2,m=3,β=32, y1=(−1,0),y2=(1,0),y3=(0,1.5), and ε=0.6,k=2. A common Fréchet minimizer is c=(0,0.5247333004905787534…). Let A=(0.58,−0.94),B=(−0.58,−0.94),M=(0,−0.94). For Gj(t)=∑i≠j‖t−yi‖3/2 and the corresponding dj, the resulting comparison margins and vote counts are d1−G1d2−G2d3−G3N(⋅)A+0.232492290498−1.215637558202+1.3036116879452B−1.215637558202+0.232492290498+1.3036116879452M−0.170268407563−0.170268407563+1.6254126190661 Because k=2 and equality would count as a vote, these strict numerical margins imply A,B∈Γε,M∉Γε. The witness therefore supplies an accepted-endpoints/rejected-midpoint nonconvexity check with nonzero margins.