Shape Operator Metric for Surface Normal Approximation

This work deals with the problem of practical mesh generation for surface normal approximation. Part of its contribution is in presenting previous work in a unified framework. A new algorithm for surface normal approximation is then introduced which improves upon existing ones in a number of aspects. In particular, it produces better approximations of surfaces both in practice and in the theoretical limit regime. Additionally, it resolves in a simple way some of the problems that previous methods for surface approximation suffered from.Engineering and Applied Science

The geometry of nonlinear least squares with applications to sloppy models and optimization

Parameter estimation by nonlinear least squares minimization is a common problem with an elegant geometric interpretation: the possible parameter values of a model induce a manifold in the space of data predictions. The minimization problem is then to find the point on the manifold closest to the data. We show that the model manifolds of a large class of models, known as sloppy models, have many universal features; they are characterized by a geometric series of widths, extrinsic curvatures, and parameter-effects curvatures. A number of common difficulties in optimizing least squares problems are due to this common structure. First, algorithms tend to run into the boundaries of the model manifold, causing parameters to diverge or become unphysical. We introduce the model graph as an extension of the model manifold to remedy this problem. We argue that appropriate priors can remove the boundaries and improve convergence rates. We show that typical fits will have many evaporated parameters. Second, bare model parameters are usually ill-suited to describing model behavior; cost contours in parameter space tend to form hierarchies of plateaus and canyons. Geometrically, we understand this inconvenient parametrization as an extremely skewed coordinate basis and show that it induces a large parameter-effects curvature on the manifold. Using coordinates based on geodesic motion, these narrow canyons are transformed in many cases into a single quadratic, isotropic basin. We interpret the modified Gauss-Newton and Levenberg-Marquardt fitting algorithms as an Euler approximation to geodesic motion in these natural coordinates on the model manifold and the model graph respectively. By adding a geodesic acceleration adjustment to these algorithms, we alleviate the difficulties from parameter-effects curvature, improving both efficiency and success rates at finding good fits.Comment: 40 pages, 29 Figure

Spectral Numerical Exterior Calculus Methods for Differential Equations on Radial Manifolds

We develop exterior calculus approaches for partial differential equations on radial manifolds. We introduce numerical methods that approximate with spectral accuracy the exterior derivative $\mathbf{d}$, Hodge star $\star$, and their compositions. To achieve discretizations with high precision and symmetry, we develop hyperinterpolation methods based on spherical harmonics and Lebedev quadrature. We perform convergence studies of our numerical exterior derivative operator $\overline{\mathbf{d}}$ and Hodge star operator $\overline{\star}$ showing each converge spectrally to $\mathbf{d}$ and $\star$. We show how the numerical operators can be naturally composed to formulate general numerical approximations for solving differential equations on manifolds. We present results for the Laplace-Beltrami equations demonstrating our approach.Comment: 22 pages, 13 figure

Normal gravity field in relativistic geodesy

Modern geodesy is subject to a dramatic change from the Newtonian paradigm to Einstein's theory of general relativity. This is motivated by the ongoing advance in development of quantum sensors for applications in geodesy including quantum gravimeters and gradientometers, atomic clocks and fiber optics for making ultra-precise measurements of the geoid and multipolar structure of the Earth's gravitational field. At the same time, VLBI, SLR, and GNSS have achieved an unprecedented level of accuracy in measuring coordinates of the reference points of the ITRF and the world height system. The main geodetic reference standard is a normal gravity field represented in the Newtonian gravity by the field of a Maclaurin ellipsoid. The present paper extends the concept of the normal gravity field to the realm of general relativity. We focus our attention on the calculation of the first post-Newtonian approximation of the normal field that is sufficient for applications. We show that in general relativity the level surface of the uniformly rotating fluid is no longer described by the Maclaurin ellipsoid but is an axisymmetric spheroid of the forth order. We parametrize the mass density distribution and derive the post-Newtonian normal gravity field of the rotating spheroid which is given in a closed form by a finite number of the ellipsoidal harmonics. We employ transformation from the ellipsoidal to spherical coordinates to deduce the post-Newtonian multipolar expansion of the metric tensor given in terms of scalar and vector gravitational potentials of the rotating spheroid. We compare these expansions with that of the normal gravity field generated by the Kerr metric and demonstrate that the Kerr metric has a fairly limited application in relativistic geodesy. Finally, we derive the post-Newtonian generalization of the Somigliana formula for the gravity field on the reference ellipsoid.Comment: 39 pages, no figures, accepted to Physical Review

`Thermodynamics' of Minimal Surfaces and Entropic Origin of Gravity

Deformations of minimal surfaces lying in constant time slices in static space-times are studied. An exact and universal formula for a change of the area of a minimal surface under shifts of nearby point-like particles is found. It allows one to introduce a local temperature on the surface and represent variations of its area in a thermodynamical form by assuming that the entropy in the Planck units equals the quarter of the area. These results provide a strong support to a recent hypothesis that gravity has an entropic origin, the minimal surfaces being a sort of holographic screens. The gravitational entropy also acquires a definite physical meaning related to quantum entanglement of fundamental degrees of freedom across the screen.Comment: 12 pages, 1 figur

Geodesics of Random Riemannian Metrics

We analyze the disordered Riemannian geometry resulting from random perturbations of the Euclidean metric. We focus on geodesics, the paths traced out by a particle traveling in this quenched random environment. By taking the point of the view of the particle, we show that the law of its observed environment is absolutely continuous with respect to the law of the random metric, and we provide an explicit form for its Radon-Nikodym derivative. We use this result to prove a "local Markov property" along an unbounded geodesic, demonstrating that it eventually encounters any type of geometric phenomenon. We also develop in this paper some general results on conditional Gaussian measures. Our Main Theorem states that a geodesic chosen with random initial conditions (chosen independently of the metric) is almost surely not minimizing. To demonstrate this, we show that a minimizing geodesic is guaranteed to eventually pass over a certain "bump surface," which locally has constant positive curvature. By using Jacobi fields, we show that this is sufficient to destabilize the minimizing property.Comment: 55 pages. Supplementary material at arXiv:1206.494