Optimal Finite Homogeneous Sphere Approximation

Omer Lavi

DOI: https://doi.org/10.1007/s00454-022-00377-w

2022-03-07

Abstract:The two-dimensional sphere can't be approximated by finite homogeneous spaces. We describe the optimal approximation and its distance from the sphere. We compare this distance to the distance achieved by all Platonic and Archimedean solids.

mathematics,computer science, theory & methods

What problem does this paper attempt to address?

A Class of Optimal Transportation Problems on the Sphere

Qirui Li,Xujia Wang

DOI: https://doi.org/10.1360/n012017-00061

2018-01-01

Abstract:In this paper,we study the optimal transportation problem on the unit sphere Sn,with the cost function c(x,y) =F(d(x,y)),where d(x,y) is the spherical distance between x,y ∈ Sn.We allow that F is a bounded function defined only in a subinterval of [0,π],while c takes-∞ when F is not defined.Under some suitable conditions,we show that there is a unique optimal mapping solving the problem.In particular,if F(d) =log(kcosd-1),k ＞ 1,or F(d) =logcosd,then the associated optimal transportation problem is respectively equivalent to the refractor problem in geometric optics,and the Aleksandrov problem in convex geometry.
Variational Sphere Set Approximation for Solid Objects

Rui Wang,Kun Zhou,John Snyder,Xinguo Liu,Hujun Bao,Qunsheng Peng,Baining Guo

DOI: https://doi.org/10.1007/s00371-006-0052-0

2006-01-01

Abstract:We approximate a solid object represented as a triangle mesh by a bounding set of spheres having minimal summed volume outside the object. We show how outside volume for a single sphere can be computed using a simple integration over the object’s triangles. We then minimize the total outside volume over all spheres in the set using a variant of iterative Lloyd clustering that splits the mesh points into sets and bounds each with an outside volume-minimizing sphere. The resulting sphere sets are tighter than those of previous methods. In experiments comparing against a state-of-the-art alternative (adaptive medial axis), our method often requires half as many spheres, or fewer, to obtain the same error, under a variety of error metrics including total outside volume, shadowing fidelity, and proximity measurement.
Computing weak distance between the 2-sphere and its nonsmooth approximations

Kazuki Koga

2024-01-13

Abstract:A novel algorithm is proposed for quantitative comparisons between compact surfaces embedded in the three-dimensional Euclidian space. The key idea is to identify those objects with the associated surface measures and compute a weak distance between them using the Fourier transform on the ambient space. In particular, the inhomogeneous Sobolev norm of negative order for a difference between two surface measures is evaluated via the Plancherel theorem, which amounts to approximating a weighted integral norm of smooth data on the frequency space. This approach allows several advantages including high accuracy due to fast-converging numerical quadrature rules, acceleration by the nonuniform fast Fourier transform, and parallelization on many-core processors. In numerical experiments, the 2-sphere, which is an example whose Fourier transform is explicitly known, is compared with its icosahedral discretization, and it is observed that the piecewise linear approximations converge to the smooth object at the quadratic rate up to small truncation.

Numerical Analysis
Central Limit Theorem for Diophantine approximation on spheres

Zouhair Ouaggag

2024-09-04

Abstract:We prove a Central Limit Theorem and an effective estimate for the counting function of Diophantine approximants on the sphere S$^n$ using homogeneous dynamics on the space of orthogonal lattices.

Number Theory,Dynamical Systems
Optimal polynomial meshes and Caratheodory-Tchakaloff submeshes on the sphere

Paul Leopardi,Alvise Sommariva,Marco Vianello

DOI: https://doi.org/10.48550/arXiv.1612.04952

2016-12-15

Numerical Analysis

Abstract:Using the notion of Dubiner distance, we give an elementary proof of the fact that good covering point configurations on the 2-sphere are optimal polynomial meshes. From these we extract Caratheodory-Tchakaloff (CATCH) submeshes for compressed Least Squares fitting.
Approximation Algorithms for Cutting a Convex Polyhedron out of a Sphere

Xuehou Tan,Gangshan Wu

DOI: https://doi.org/10.1016/j.tcs.2012.03.035

IF: 1.002

2013-01-01

Theoretical Computer Science

Abstract:This paper presents the following approximation algorithms for computing a minimum cost sequence of planes to cut a convex polyhedron P of n vertices out of a sphere Q: an O(n logn) time O(log2 n)-factor approximation, an O(n 1.5 logn) time O(logn)-factor approximation, and an O(1)-factor approximation with exponential running time. Our results significantly improve upon the previous O(n 3) time O(log2 n)-factor approximation solution.
An Optimal Ansatz Space for Moving Least Squares Approximation on Spheres

Ralf Hielscher,Tim Pöschl

DOI: https://doi.org/10.1007/s10444-024-10201-z

2024-10-24

Abstract:We revisit the moving least squares (MLS) approximation scheme on the sphere $\mathbb S^{d-1} \subset \mathbb R^d$, where $d>1$. It is well known that using the spherical harmonics up to degree $L \in \mathbb N$ as ansatz space yields for functions in $\mathcal C^{L+1}(\mathbb S^{d-1})$ the approximation order $\mathcal O \left( h^{L+1} \right)$, where $h$ denotes the fill distance of the sampling nodes. In this paper we show that the dimension of the ansatz space can be almost halved, by including only spherical harmonics of even or odd degree up to $L$, while preserving the same order of approximation. Numerical experiments indicate that using the reduced ansatz space is essential to ensure the numerical stability of the MLS approximation scheme as $h \to 0$. Finally, we compare our approach with an MLS approximation scheme that uses polynomials on the tangent space as ansatz space.

Numerical Analysis
Algorithm for Fast Calculating the Nearest Distance Between Space Point and Arbitrary Polyhedron

方向,鲍虎军,王平安,彭群生

DOI: https://doi.org/10.3321/j.issn:1003-9775.2001.09.005

2001-01-01

Journal of Computer-Aided Design & Computer Graphics

Abstract:The algorithm uses dynamic sphere searching technology which takes the given point as sphere center. A potential facet set can be found quickly and accurately and this set is compact enough to accelerate the distance calculation greatly. Compared with the hierarchical presentation approach, this algorithm avoids frequently calculating the distance between the point and hierarchical structure, thus being suitable to be applied to the cases where a great deal of distance calculations are needed, such as distance fields calculation, collision detection in virtual environment and interference test of robots and NC machining.
Generating Probability Distributions on Intervals and Spheres with Application to Finite Element Method.

Yizhi Sun,Zhilin Sun

DOI: https://doi.org/10.1016/j.camwa.2020.10.017

IF: 3.218

2021-01-01

Computers & Mathematics with Applications

Abstract:This work aims to build a bridge between probability methods and finite element methods. It starts with considering probability distributions supported in an interval [a,b], which incorporate the traditional probability distributions defined on the whole real space as limit cases on the one hand, lead to a type of spherical probability models with wide potential applications on the other hand. This type of probability has scaling and symmetry feature, and sufficient conditions under which a density function can be generated through discrete polynomial spectrum are given in this work followed by concrete examples. The density function ρ obtained in this way has the advantage of being positive definite. Computer based numerical simulation shows that the theoretically verified criteria for probability distribution are almost optimal with respect to our testing examples. After the establishment of an approximation theorem in L1 space, we propose a probabilistic Galerkin scheme that can be either continuous or discontinuous, which is potentially useful to asymptotically solve some PDEs on the sphere locally and globally.
A Convergent Finite Difference Method for Optimal Transport on the Sphere

Brittany Froese Hamfeldt,Axel G. R. Turnquist

DOI: https://doi.org/10.48550/arXiv.2105.03500

2021-05-08

Abstract:We introduce a convergent finite difference method for solving the optimal transportation problem on the sphere. The method applies to both the traditional squared geodesic cost (arising in mesh generation) and a logarithmic cost (arising in the reflector antenna design problem). At each point on the sphere, we replace the surface PDE with a Generated Jacobian equation posed on the local tangent plane using geodesic normal coordinates. The discretization is inspired by recent monotone methods for the Monge-Ampère equation, but requires significant adaptations in order to correctly handle the mix of gradient and Hessian terms appearing inside the nonlinear determinant operator, as well as the singular logarithmic cost function. Numerical results demonstrate the success of this method on a wide range of challenging problems involving both the squared geodesic and the logarithmic cost functions.

Numerical Analysis
Fast Geometric Fit Algorithm for Sphere Using Exact Solution

Sumith YD

DOI: https://doi.org/10.48550/arXiv.1506.02776

2015-06-09

Computer Vision and Pattern Recognition

Abstract:Sphere fitting is a common problem in almost all science and engineering disciplines. Most of methods available are iterative in behavior. This involves fitting of the parameters in a least square sense or in a geometric sense. Here we extend the methods of Thomas Chan and Landau who fitted the 2D data using circle. This work closely resemble their work in redefining the error estimate and solving the sphere fitting problem exactly. The solutions for center and radius of the sphere can be found exactly and the equations can be hard coded for high performance. We have also shown some comparison with other popular methods and how this method behaves.
Polyhedral Approximation of Spectrahedral Shadows via Homogenization

Daniel Dörfler,Andreas Löhne

DOI: https://doi.org/10.1007/s10957-023-02363-5

2024-01-18

Journal of Optimization Theory and Applications

Abstract:This article is concerned with the problem of approximating a not necessarily bounded spectrahedral shadow, a certain convex set, by polyhedra. By identifying the set with its homogenization, the problem is reduced to the approximation of a closed convex cone. We introduce the notion of homogeneous -approximation of a convex set and show that it defines a meaningful concept in the sense that approximations converge to the original set if the approximation error diminishes. Moreover, we show that a homogeneous -approximation of the polar of a convex set is immediately available from an approximation of the set itself under mild conditions. Finally, we present an algorithm for the computation of homogeneous -approximations of spectrahedral shadows and demonstrate it on examples.

mathematics, applied,operations research & management science
Geometric Locations of Points Equally Distance from Two Given Geometric Figures. Part 4: Geometric Locations of Points Equally Remote from Two Spheres

Vladimir Vyshnyepolskiy,E. Zavarihina,D. Peh

DOI: https://doi.org/10.12737/2308-4898-2021-9-3-12-29

2021-12-09

Abstract:The article deals with the geometric locations of points equidistant from two spheres. In all variants of the mutual position of the spheres, the geometric places of the points are two surfaces. When the centers of the spheres coincide with the locus of points equidistant from the spheres, there will be spheres equal to the half-sum and half-difference of the diameters of the original spheres. In three variants of the relative position of the initial spheres, one of the two surfaces of the geometric places of the points is a two-sheet hyperboloid of revolution. It is obtained when: 1) the spheres intersect, 2) the spheres touch, 3) the outer surfaces of the spheres are removed from each other. In the case of equal spheres, a two-sheeted hyperboloid of revolution degenerates into a two-sheeted plane, more precisely, it is a second-order degenerate surface with a second infinitely distant branch. The spheres intersect - the second locus of the points will be the ellipsoid of revolution. Spheres touch - the second locus of points - an ellipsoid of revolution, degenerated into a straight line, more precisely into a zero-quadric of the second order - a cylindrical surface with zero radius. The outer surfaces of the spheres are distant from each other - the second locus of points will be a two-sheet hyperboloid of revolution. The small sphere is located inside the large one - two coaxial confocal ellipsoids of revolution. In all variants of the mutual position of spheres of the same diameters, the common geometrical place of equidistant points is a plane (degenerate surface of the second order) passing through the middle of the segment perpendicular to it, connecting the centers of the original spheres. The second locus of points equidistant from two spheres of the same diameter can be either an ellipsoid of revolution (if the original spheres intersect), or a straight (cylindrical surface with zero radius) connecting the centers of the original spheres when the original spheres touch each other, or a two-sheet hyperboloid of revolution (if continue to increase the distance between the centers of the original spheres).
A hierarchy of eigencomputations for polynomial optimization on the sphere

Nathaniel Johnston,Benjamin Lovitz,Aravindan Vijayaraghavan

2023-10-27

Abstract:We introduce a convergent hierarchy of lower bounds on the minimum value of a real homogeneous polynomial over the sphere. The main practical advantage of our hierarchy over the sum-of-squares (SOS) hierarchy is that the lower bound at each level of our hierarchy is obtained by a minimum eigenvalue computation, as opposed to the full semidefinite program (SDP) required at each level of SOS. In practice, this allows us to go to much higher levels than are computationally feasible for the SOS hierarchy. For both hierarchies, the underlying space at the $k$-th level is the set of homogeneous polynomials of degree $2k$. We prove that our hierarchy converges as $O(1/k)$ in the level $k$, matching the best-known convergence of the SOS hierarchy when the number of variables $n$ is less than the half-degree $d$ (the best-known convergence of SOS when $n \geq d$ is $O(1/k^2)$). More generally, we introduce a convergent hierarchy of minimum eigenvalue computations for minimizing the inner product between a real tensor and an element of the spherical Segre-Veronese variety, with similar convergence guarantees. As examples, we obtain hierarchies for computing the (real) tensor spectral norm, and for minimizing biquadratic forms over the sphere. Hierarchies of eigencomputations for more general constrained polynomial optimization problems are discussed.

Optimization and Control,Data Structures and Algorithms,Algebraic Geometry,Quantum Physics
The Sphere Formula

Oliver Knill

DOI: https://doi.org/10.48550/arXiv.2301.05736

2023-01-14

Abstract:The sphere formula states that in an arbitrary finite abstract simplicial complex, the sum of the Euler characteristic of unit spheres centered at even-dimensional simplices is equal to the sum of the Euler characteristic of unit spheres centered at odd-dimensional simplices. It follows that if a geometry has constant unit sphere Euler characteristic, like a manifold, then all its unit spheres have zero Euler characteristic or the space itself has zero Euler characteristic. Especially, odd-dimensional manifolds have zero Euler characteristic, a fact usually verified either in algebraic topology using Poincaré duality together with Riemann-Hurwitz then deriving it from the existence of a Morse function, using that the Morse indices of the function and its negative add up to zero in odd dimensions. Gauss Bonnet also shows that odd-dimensional Dehn-Sommerville spaces have zero Euler characteristic because they have constant zero curvature. Zero curvature phenomenons can be understood integral geometrically as index expectation or as Dehn-Sommerville relations.

Combinatorics,Discrete Mathematics,General Topology
An Optimal Rigidity Theorem for Complete Submanifolds in a Sphere

Hong-wei Xu,Jiao-feng Zhu

DOI: https://doi.org/10.1007/s11766-008-0212-2

2008-01-01

Abstract:It is proved that if M n is an n-dimensional complete submanifold with parallel mean curvature vector and flat normal bundle in S n+p (1), and if sup M S < α(n,H), where $\alpha (n,{\rm H}) = n + \frac{{n^3 }}{{2(n - 1)}}{\rm H}^2 - \frac{{n(n - 2)}}{{2(n - 1)}}\sqrt {n^2 H^4 + 4(n - 1)H^2 } $, then M n must be the totally umbilical sphere $S^n (\frac{1}{{\sqrt {1 + H^2 } }})$. An example to show that the pinching constant α(n,H) appears optimal is given.
Some novel constructions of optimal Gromov-Hausdorff-optimal correspondences between spheres

Saúl Rodríguez Martín

2024-09-04

Abstract:In this article, as a first contribution, we provide alternative proofs of recent results by Harrison and Jeffs which determine the precise value of the Gromov-Hausdorff (GH) distance between the circle $\mathbb{S}^1$ and the $n$-dimensional sphere $\mathbb{S}^n$ (for any $n\in\mathbb{N}$) when endowed with their respective geodesic metrics. Additionally, we prove that the GH distance between $\mathbb{S}^3$ and $\mathbb{S}^4$ is equal to $\frac{1}{2}\arccos\left(\frac{-1}{4}\right)$, thus settling the case $n=3$ of a conjecture by Lim, Mémoli and Smith.

Metric Geometry
Smooth Distance Approximation

Ahmed Abdelkader,David M. Mount

DOI: https://doi.org/10.48550/arXiv.2308.08791

2023-08-17

Abstract:Traditional problems in computational geometry involve aspects that are both discrete and continuous. One such example is nearest-neighbor searching, where the input is discrete, but the result depends on distances, which vary continuously. In many real-world applications of geometric data structures, it is assumed that query results are continuous, free of jump discontinuities. This is at odds with many modern data structures in computational geometry, which employ approximations to achieve efficiency, but these approximations often suffer from discontinuities. In this paper, we present a general method for transforming an approximate but discontinuous data structure into one that produces a smooth approximation, while matching the asymptotic space efficiencies of the original. We achieve this by adapting an approach called the partition-of-unity method, which smoothly blends multiple local approximations into a single smooth global approximation. We illustrate the use of this technique in a specific application of approximating the distance to the boundary of a convex polytope in $\mathbb{R}^d$ from any point in its interior. We begin by developing a novel data structure that efficiently computes an absolute $\varepsilon$-approximation to this query in time $O(\log (1/\varepsilon))$ using $O(1/\varepsilon^{d/2})$ storage space. Then, we proceed to apply the proposed partition-of-unity blending to guarantee the smoothness of the approximate distance field, establishing optimal asymptotic bounds on the norms of its gradient and Hessian.

Computational Geometry
Sphere Packings in Euclidean Space with Forbidden Distances

Felipe Gonçalves,Guilherme Vedana

2024-02-15

Abstract:In this paper we study the sphere packing problem in Euclidean space where we impose additional constraints on the separations of the center points. We prove that any sphere packing in dimension $48$, with spheres of radii $r$, such that \emph{no} two centers $x_1$ and $x_2$ satisfy $\sqrt{\tfrac{4}{3}} < \frac{1}{2r}|x_1-x_2| <\sqrt{\tfrac{5}{3}}$, has center density less or equal than $(3/2)^{24}$. Equality occurs for lattice packings if and only if the packing is given by a $48$-dimensional even unimodular extremal lattice. This shows that any of the lattices $P_{48p},P_{48q},P_{48m}$ and $P_{48n}$ are optimal for this constrained packing problem, and gives evidence towards the conjecture that extremal lattices are optimal unconstrained sphere packings in $48$ dimensions. We also provide results for packings up to dimension $d\leq 1200$, where we impose constraints on the distance between centers and on the minimal norm of the spectrum, showing that even unimodular extremal lattices are again uniquely optimal. Moreover, in the $1$-dimensional case, where is not at all clear that periodic packings are among those with largest density, we nevertheless give a condition on the set of constraints that allows this to happen, and we develop an algorithm to find these periodic configurations by relating the problem to a question about dominos.

Number Theory,Classical Analysis and ODEs,Metric Geometry
Spherical Designs for Function Approximation and Beyond

Yuchen Xiao,Xiaosheng Zhuang

2023-11-30

Abstract:In this paper, we compare two optimization algorithms using full Hessian and approximation Hessian to obtain numerical spherical designs through their variational characterization. Based on the obtained spherical design point sets, we investigate the approximation of smooth and non-smooth functions by spherical harmonics with spherical designs. Finally, we use spherical framelets for denoising Wendland functions as an application, which shows the great potential of spherical designs in spherical data processing.

Numerical Analysis,Signal Processing

Optimal Finite Homogeneous Sphere Approximation

A Class of Optimal Transportation Problems on the Sphere

Variational Sphere Set Approximation for Solid Objects

Computing weak distance between the 2-sphere and its nonsmooth approximations

Central Limit Theorem for Diophantine approximation on spheres

Optimal polynomial meshes and Caratheodory-Tchakaloff submeshes on the sphere

Approximation Algorithms for Cutting a Convex Polyhedron out of a Sphere

An Optimal Ansatz Space for Moving Least Squares Approximation on Spheres

Algorithm for Fast Calculating the Nearest Distance Between Space Point and Arbitrary Polyhedron

Generating Probability Distributions on Intervals and Spheres with Application to Finite Element Method.

A Convergent Finite Difference Method for Optimal Transport on the Sphere

Fast Geometric Fit Algorithm for Sphere Using Exact Solution

Polyhedral Approximation of Spectrahedral Shadows via Homogenization

Geometric Locations of Points Equally Distance from Two Given Geometric Figures. Part 4: Geometric Locations of Points Equally Remote from Two Spheres

A hierarchy of eigencomputations for polynomial optimization on the sphere

The Sphere Formula

An Optimal Rigidity Theorem for Complete Submanifolds in a Sphere

Some novel constructions of optimal Gromov-Hausdorff-optimal correspondences between spheres

Smooth Distance Approximation

Sphere Packings in Euclidean Space with Forbidden Distances

Spherical Designs for Function Approximation and Beyond