Jacobi Mapping Approach for a Precise Cosmological Weak Lensing Formalism

Cosmological weak lensing has been a highly successful and rapidly developing research field since the first detection of cosmic shear in 2000. However, it has recently been pointed out in Yoo et al. that the standard weak lensing formalism yields gauge-dependent results and, hence, does not meet the level of accuracy demanded by the next generation of weak lensing surveys. Here, we show that the Jacobi mapping formalism provides a solid alternative to the standard formalism, as it accurately describes all the relativistic effects contributing to the weak lensing observables. We calculate gauge-invariant expressions for the distortion in the luminosity distance, the cosmic shear components and the lensing rotation to linear order including scalar, vector and tensor perturbations. In particular, the Jacobi mapping formalism proves that the rotation is fully vanishing to linear order. Furthermore, the cosmic shear components contain an additional term in tensor modes which is absent in the results obtained with the standard formalism. Our work provides further support and confirmation of the gauge-invariant lensing formalism needed in the era of precision cosmology.Comment: 33 pages, no figures, published in JCA

Sympiler: Transforming Sparse Matrix Codes by Decoupling Symbolic Analysis

Sympiler is a domain-specific code generator that optimizes sparse matrix computations by decoupling the symbolic analysis phase from the numerical manipulation stage in sparse codes. The computation patterns in sparse numerical methods are guided by the input sparsity structure and the sparse algorithm itself. In many real-world simulations, the sparsity pattern changes little or not at all. Sympiler takes advantage of these properties to symbolically analyze sparse codes at compile-time and to apply inspector-guided transformations that enable applying low-level transformations to sparse codes. As a result, the Sympiler-generated code outperforms highly-optimized matrix factorization codes from commonly-used specialized libraries, obtaining average speedups over Eigen and CHOLMOD of 3.8X and 1.5X respectively.Comment: 12 page

The orbit rigidity matrix of a symmetric framework

A number of recent papers have studied when symmetry causes frameworks on a graph to become infinitesimally flexible, or stressed, and when it has no impact. A number of other recent papers have studied special classes of frameworks on generically rigid graphs which are finite mechanisms. Here we introduce a new tool, the orbit matrix, which connects these two areas and provides a matrix representation for fully symmetric infinitesimal flexes, and fully symmetric stresses of symmetric frameworks. The orbit matrix is a true analog of the standard rigidity matrix for general frameworks, and its analysis gives important insights into questions about the flexibility and rigidity of classes of symmetric frameworks, in all dimensions. With this narrower focus on fully symmetric infinitesimal motions, comes the power to predict symmetry-preserving finite mechanisms - giving a simplified analysis which covers a wide range of the known mechanisms, and generalizes the classes of known mechanisms. This initial exploration of the properties of the orbit matrix also opens up a number of new questions and possible extensions of the previous results, including transfer of symmetry based results from Euclidean space to spherical, hyperbolic, and some other metrics with shared symmetry groups and underlying projective geometry.Comment: 41 pages, 12 figure

Structure Preserving Parallel Algorithms for Solving the Bethe-Salpeter Eigenvalue Problem

The Bethe-Salpeter eigenvalue problem is a dense structured eigenvalue problem arising from discretized Bethe-Salpeter equation in the context of computing exciton energies and states. A computational challenge is that at least half of the eigenvalues and the associated eigenvectors are desired in practice. We establish the equivalence between Bethe-Salpeter eigenvalue problems and real Hamiltonian eigenvalue problems. Based on theoretical analysis, structure preserving algorithms for a class of Bethe-Salpeter eigenvalue problems are proposed. We also show that for this class of problems all eigenvalues obtained from the Tamm-Dancoff approximation are overestimated. In order to solve large scale problems of practical interest, we discuss parallel implementations of our algorithms targeting distributed memory systems. Several numerical examples are presented to demonstrate the efficiency and accuracy of our algorithms