Fourier, Gauss, Fraunhofer, Porod and the Shape from Moments Problem

We show how the Fourier transform of a shape in any number of dimensions can be simplified using Gauss's law and evaluated explicitly for polygons in two dimensions, polyhedra three dimensions, etc. We also show how this combination of Fourier and Gauss can be related to numerous classical problems in physics and mathematics. Examples include Fraunhofer diffraction patterns, Porods law, Hopfs Umlaufsatz, the isoperimetric inequality and Didos problem. We also use this approach to provide an alternative derivation of Davis's extension of the Motzkin-Schoenberg formula to polygons in the complex plane.Comment: 21 pages, no figure

Large-Scale Assessment of Polygon-Edge Boulder Clustering in the Martian Northern Lowlands

Two features evident in many images of the martian northern low-lands are polygonal fractures (especially northwards of 60N) and meter-scale surface boulders. Since their first observation, several attempts have been made to classify and study these polygons as well as how the forces that form these polygons may modify the surface. Surface boulders have been used as a potential indicator of such modification, though current studies find evidence both for and against their association with the underlying polygons. Both these investigations are limited by the same fundamental challenge: map-ping the location of surface boulders manually is not practical at large scales. Here, we use the Martian Boulder Automatic Recognition System (MBARS) to provide image-wide assessments of boulder location and size, enabling large-scale assessment of boulder populations. To compare these boulder locations with the underlying polygons, we modified the 2-D Fourier analysis described by Orloff in 2013 to analyze boulder locations. When compared with Orloffs observations of polygon scales, this provides an avenue for large-scale comparison of boulder-cluster scale and polygon scale

An Exponentially Convergent Nonpolynomial Finite Element Method for Time-Harmonic Scattering from Polygons

In recent years nonpolynomial finite element methods have received increasing attention for the efficient solution of wave problems. As with their close cousin the method of particular solutions, high efficiency comes from using solutions to the Helmholtz equation as basis functions. We present and analyze such a method for the scattering of two-dimensional scalar waves from a polygonal domain that achieves exponential convergence purely by increasing the number of basis functions in each element. Key ingredients are the use of basis functions that capture the singularities at corners and the representation of the scattered field towards infinity by a combination of fundamental solutions. The solution is obtained by minimizing a least-squares functional, which we discretize in such a way that a matrix least-squares problem is obtained. We give computable exponential bounds on the rate of convergence of the least-squares functional that are in very good agreement with the observed numerical convergence. Challenging numerical examples, including a nonconvex polygon with several corner singularities, and a cavity domain, are solved to around 10 digits of accuracy with a few seconds of CPU time. The examples are implemented concisely with MPSpack, a MATLAB toolbox for wave computations with nonpolynomial basis functions, developed by the authors. A code example is included

Pruned Continuous Haar Transform of 2D Polygonal Patterns with Application to VLSI Layouts

We introduce an algorithm for the efficient computation of the continuous Haar transform of 2D patterns that can be described by polygons. These patterns are ubiquitous in VLSI processes where they are used to describe design and mask layouts. There, speed is of paramount importance due to the magnitude of the problems to be solved and hence very fast algorithms are needed. We show that by techniques borrowed from computational geometry we are not only able to compute the continuous Haar transform directly, but also to do it quickly. This is achieved by massively pruning the transform tree and thus dramatically decreasing the computational load when the number of vertices is small, as is the case for VLSI layouts. We call this new algorithm the pruned continuous Haar transform. We implement this algorithm and show that for patterns found in VLSI layouts the proposed algorithm was in the worst case as fast as its discrete counterpart and up to 12 times faster.Comment: 4 pages, 5 figures, 1 algorith

The Enumeration of Prudent Polygons by Area and its Unusual Asymptotics

Prudent walks are special self-avoiding walks that never take a step towards an already occupied site, and \emph{$k$-sided prudent walks} (with $k=1,2,3,4$) are, in essence, only allowed to grow along $k$ directions. Prudent polygons are prudent walks that return to a point adjacent to their starting point. Prudent walks and polygons have been previously enumerated by length and perimeter (Bousquet-M\'elou, Schwerdtfeger; 2010). We consider the enumeration of \emph{prudent polygons} by \emph{area}. For the 3-sided variety, we find that the generating function is expressed in terms of a $q$-hypergeometric function, with an accumulation of poles towards the dominant singularity. This expression reveals an unusual asymptotic structure of the number of polygons of area $n$, where the critical exponent is the transcendental number $\log_23$ and and the amplitude involves tiny oscillations. Based on numerical data, we also expect similar phenomena to occur for 4-sided polygons. The asymptotic methodology involves an original combination of Mellin transform techniques and singularity analysis, which is of potential interest in a number of other asymptotic enumeration problems.Comment: 27 pages, 6 figure