In this work, we introduce two sets of algorithms inspired by the ideas from modern geometry. One is computational conformal geometry method, including harmonic maps, holomorphic 1-forms and Ricci flow. The other one is optimization method using affine normals.

No abstract uploaded!

This article focuses on the discrete double-curl operator arising in the Maxwell equation that models three dimensional photonic crystals with face centered cubic lattice. The discrete double-curl operator is the degenerate coefficient matrix of the generalized eigenvalue problems (GEVP) due to the Maxwell equation. We derive an eigendecomposition of the degenerate coefficient matrix and explore an explicit form of orthogonal basis for the range and null spaces of this matrix. To solve the GEVP, we apply these theoretical results to project the GEVP to a standard eigenvalue problem (SEVP), which involves only the eigenspace associated with the nonzero eigenvalues of the GEVP and therefore the zero eigenvalues are excluded and will not degrade the computational efficiency. This projected SEVP can be solved efficiently by the inverse Lanczos method. The linear systems within the inverse Lanczos method are well-conditioned and can be solved efficiently by the conjugate gradient method without using a preconditioner. We also demonstrate how two forms of matrix-vector multiplications, which are the most costly part of the inverse Lanczos method, can be computed by fast Fourier transformation due to the eigendecomposition to significantly reduce the computation cost. Integrating all of these findings and techniques, we obtain a fast eigenvalue solver. The solver has been implemented by MATLAB and successfully solves each of a set of $5.184$ million dimension eigenvalue problems within $50$ to $104$ minutes on a workstation with two Intel Quad-C ore Xeon X5687 3.6 GHz CPUs.

Wedge product on deRham complex of a Riemannian manifold M can be pulled back to H^∗(M) via explicit homotopy, constructed using Green’s operator, to give higher product structures. We prove Fukaya’s conjecture which suggests that Witten deformation of these higher product structures have semiclassical limits as operators defined by counting gradient flow trees with respect to Morse functions, which generalizes the remarkable Witten deformation of deRham differential from a statement concerning homology to one concerning real homotopy type of M. Various applications of this conjecture to mirror symmetry are also suggested by Fukaya.

For a sub-Riemannian manifold provided with a smooth volume, we relate the small-time asymptotics of the heat kernel at a point y of the cut locus from x with roughly “how much” y is conjugate to x. This is done under the hypothesis that all minimizers connecting x to y are strongly normal, i.e. all pieces of the trajectory are not abnormal. Our result is a refinement of the one of Leandre 4t log pt(x, y) ! −d2(x, y) for t ! 0, in which only the leading exponential term is detected. Our results are obtained by extending an idea of Molchanov from the Riemannian to the sub-Riemannian case, and some details we get appear to be new even in the Riemannian context. These results permit us to obtain properties of the sub-Riemannian distance starting from those of the heat kernel and vice versa. For the Grushin plane endowed with the Euclidean volume, we get the expansion pt(x, y)  t−5/4 exp(−d2(x, y)/4t) where y is reached from a - Riemannian point x by a minimizing geodesic which is conjugate at y.