skip to main content

Title: On the Reconstruction of Geodesic Subspaces of ℝ^N
We consider the topological and geometric reconstruction of a geodesic subspace of [Formula: see text] both from the Čech and Vietoris-Rips filtrations on a finite, Hausdorff-close, Euclidean sample. Our reconstruction technique leverages the intrinsic length metric induced by the geodesics on the subspace. We consider the distortion and convexity radius as our sampling parameters for the reconstruction problem. For a geodesic subspace with finite distortion and positive convexity radius, we guarantee a correct computation of its homotopy and homology groups from the sample. This technique provides alternative sampling conditions to the existing and commonly used conditions based on weak feature size and [Formula: see text]–reach, and performs better under certain types of perturbations of the geodesic subspace. For geodesic subspaces of [Formula: see text], we also devise an algorithm to output a homotopy equivalent geometric complex that has a very small Hausdorff distance to the unknown underlying space.
; ; ;
Award ID(s):
1854336 1664858 2046730
Publication Date:
Journal Name:
International Journal of Computational Geometry & Applications
Page Range or eLocation-ID:
1 to 27
Sponsoring Org:
National Science Foundation
More Like this
  1. Let [Formula: see text] be a group acting properly and by isometries on a metric space [Formula: see text]; it follows that the quotient or orbit space [Formula: see text] is also a metric space. We study the Vietoris–Rips and Čech complexes of [Formula: see text]. Whereas (co)homology theories for metric spaces let the scale parameter of a Vietoris–Rips or Čech complex go to zero, and whereas geometric group theory requires the scale parameter to be sufficiently large, we instead consider intermediate scale parameters (neither tending to zero nor to infinity). As a particular case, we study the Vietoris–Rips and Čech thickenings of projective spaces at the first scale parameter where the homotopy type changes.
  2. We consider a minimization problem whose objective function is the sum of a fidelity term, not necessarily convex, and a regularization term defined by a positive regularization parameter [Formula: see text] multiple of the [Formula: see text] norm composed with a linear transform. This problem has wide applications in compressed sensing, sparse machine learning and image reconstruction. The goal of this paper is to understand what choices of the regularization parameter can dictate the level of sparsity under the transform for a global minimizer of the resulting regularized objective function. This is a critical issue but it has been left unaddressed. We address it from a geometric viewpoint with which the sparsity partition of the image space of the transform is introduced. Choices of the regularization parameter are specified to ensure that a global minimizer of the corresponding regularized objective function achieves a prescribed level of sparsity under the transform. Results are obtained for the spacial sparsity case in which the transform is the identity map, a case that covers several applications of practical importance, including machine learning, image/signal processing and medical image reconstruction.
  3. Given a closed [Formula: see text]-dimensional submanifold [Formula: see text], encapsulated in a compact domain [Formula: see text], [Formula: see text], we consider the problem of determining the intrinsic geometry of the obstacle [Formula: see text] (such as volume, integral curvature) from the scattering data, produced by the reflections of geodesic trajectories from the boundary of a tubular [Formula: see text]-neighborhood [Formula: see text] of [Formula: see text] in [Formula: see text]. The geodesics that participate in this scattering emanate from the boundary [Formula: see text] and terminate there after a few reflections from the boundary [Formula: see text]. However, the major problem in this setting is that a ray (a billiard trajectory) may get stuck in the vicinity of [Formula: see text] by entering some trap there so that this ray will have infinitely many reflections from [Formula: see text]. To rule out such a possibility, we modify the geometry of a tube [Formula: see text] by building it from spherical bubbles. We need to use [Formula: see text] many bubbling tubes [Formula: see text] for detecting certain global invariants of [Formula: see text], invariants that reflect its intrinsic geometry. Thus, the words “layered scattering” are in the title.more »These invariants were studied by Hermann Weyl in his classical theory of tubes [Formula: see text] and their volumes.

    « less
  4. We consider a stochastic inventory control problem under censored demand, lost sales, and positive lead times. This is a fundamental problem in inventory management, with significant literature establishing near optimality of a simple class of policies called “base-stock policies” as well as the convexity of long-run average cost under those policies. We consider a relatively less studied problem of designing a learning algorithm for this problem when the underlying demand distribution is unknown. The goal is to bound the regret of the algorithm when compared with the best base-stock policy. Our main contribution is a learning algorithm with a regret bound of [Formula: see text] for the inventory control problem. Here, [Formula: see text] is the fixed and known lead time, and D is an unknown parameter of the demand distribution described roughly as the expected number of time steps needed to generate enough demand to deplete one unit of inventory. Notably, our regret bounds depend linearly on L, which significantly improves the previously best-known regret bounds for this problem where the dependence on L was exponential. Our techniques utilize the convexity of the long-run average cost and a newly derived bound on the “bias” of base-stock policies to establishmore »an almost black box connection between the problem of learning in Markov decision processes (MDPs) with these properties and the stochastic convex bandit problem. The techniques presented here may be of independent interest for other settings that involve large structured MDPs but with convex asymptotic average cost functions.« less
  5. The goal of this paper is to study limiting behavior of a self-organized continuous flock evolving according to the 1D hydrodynamic Euler Alignment model. We provide a series of quantitative estimates that show how far the density of the limiting flock is from a uniform distribution. The key quantity that controls density distortion is the entropy [Formula: see text], and the measure of deviation from uniformity is given by a well-known conserved quantity [Formula: see text], where [Formula: see text] is velocity and [Formula: see text] is the communication operator with kernel [Formula: see text]. The cases of Lipschitz, singular geometric, and topological kernels are covered in the study.