Rectangular Full Packed Format for Cholesky's Algorithm: Factorization, Solution and Inversion

We describe a new data format for storing triangular, symmetric, and Hermitian matrices called RFPF (Rectangular Full Packed Format). The standard two dimensional arrays of Fortran and C (also known as full format) that are used to represent triangular and symmetric matrices waste nearly half of the storage space but provide high performance via the use of Level 3 BLAS. Standard packed format arrays fully utilize storage (array space) but provide low performance as there is no Level 3 packed BLAS. We combine the good features of packed and full storage using RFPF to obtain high performance via using Level 3 BLAS as RFPF is a standard full format representation. Also, RFPF requires exactly the same minimal storage as packed format. Each LAPACK full and/or packed triangular, symmetric, and Hermitian routine becomes a single new RFPF routine based on eight possible data layouts of RFPF. This new RFPF routine usually consists of two calls to the corresponding LAPACK full format routine and two calls to Level 3 BLAS routines. This means {\it no} new software is required. As examples, we present LAPACK routines for Cholesky factorization, Cholesky solution and Cholesky inverse computation in RFPF to illustrate this new work and to describe its performance on several commonly used computer platforms. Performance of LAPACK full routines using RFPF versus LAPACK full routines using standard format for both serial and SMP parallel processing is about the same while using half the storage. Performance gains are roughly one to a factor of 43 for serial and one to a factor of 97 for SMP parallel times faster using vendor LAPACK full routines with RFPF than with using vendor and/or reference packed routines

Influence of optical aberrations in an atomic gyroscope

In atom interferometry based on light-induced diffraction, the optical aberrations of the laser beam splitters are a dominant source of noise and systematic effect. In an atomic gyroscope, this effect is dramatically reduced by the use of two atomic sources. But it remains critical while coupled to fluctuations of atomic trajectories, and appears as a main source of noise to the long term stability. Therefore we measure these contributions in our setup, using cold Cesium atoms and stimulated Raman transitions

First electron beam polarization measurements with a Compton polarimeter at Jefferson Laboratory

A Compton polarimeter has been installed in Hall A at Jefferson Laboratory. This letter reports on the first electron beam polarization measurements performed during the HAPPEX experiment at an electron energy of 3.3 GeV and an average current of 40 $\mu$A. The heart of this device is a Fabry-Perot cavity which increased the luminosity for Compton scattering in the interaction region so much that a 1.4% statistical accuracy could be obtained within one hour, with a 3.3% total error

Stability of axial orbits in galactic potentials

We investigate the dynamics in a galactic potential with two reflection symmetries. The phase-space structure of the real system is approximated with a resonant detuned normal form constructed with the method based on the Lie transform. Attention is focused on the stability properties of the axial periodic orbits that play an important role in galactic models. Using energy and ellipticity as parameters, we find analytical expressions of bifurcations and compare them with numerical results available in the literature.Comment: 20 pages, accepted for publication on Celestial Mechanics and Dynamical Astronom

Sensitive gravity-gradiometry with atom interferometry: progress towards an improved determination of the gravitational constant

We here present a high sensitivity gravity-gradiometer based on atom interferometry. In our apparatus, two clouds of laser-cooled rubidium atoms are launched in fountain configuration and interrogated by a Raman interferometry sequence to probe the gradient of gravity field. We recently implemented a high-flux atomic source and a newly designed Raman lasers system in the instrument set-up. We discuss the applications towards a precise determination of the Newtonian gravitational constant G. The long-term stability of the instrument and the signal-to-noise ratio demonstrated here open interesting perspectives for pushing the measurement precision below the 100 ppm level

Bohr-Sommerfeld Quantization of Periodic Orbits

We show, that the canonical invariant part of $\hbar$ corrections to the Gutzwiller trace formula and the Gutzwiller-Voros spectral determinant can be computed by the Bohr-Sommerfeld quantization rules, which usually apply for integrable systems. We argue that the information content of the classical action and stability can be used more effectively than in the usual treatment. We demonstrate the improvement of precision on the example of the three disk scattering system.Comment: revte

A new cell primo-culture method for freshwater benthic diatom communities

A new cell primo-culture method was developed for the benthic diatom community isolated from biofilm sampled in rivers. The approach comprised three steps: (1) scraping biofilm from river pebbles, (2) diatom isolation from biofilm, and (3) diatom community culture. With a view to designing a method able to stimulate the growth of diatoms, to limit the development of other microorganisms, and to maintain in culture a community similar to the original natural one, different factors were tested in step 3: cell culture medium (Chu No 10 vs Freshwater “WC” medium modified), cell culture vessel, and time of culture. The results showed that using Chu No 10 medium in an Erlenmeyer flask for cell culture was the optimal method, producing enough biomass for ecotoxicological tests as well as minimising development of other microorganisms. After 96 h of culture, communities differed from the original communities sampled in the two rivers studied. Species tolerant of eutrophic or saprobic conditions were favoured during culture. This method of diatom community culture affords the opportunity to assess, in vitro, the effects of different chemicals or effluents (water samples andindustrial effluents) on diatom communities, as well as on diatom cells, from a wide range of perspectives

Influence of lasers propagation delay on the sensitivity of atom interferometers

In atom interferometers based on two photon transitions, the delay induced by the difference of the laser beams paths makes the interferometer sensitive to the fluctuations of the frequency of the lasers. We first study, in the general case, how the laser frequency noise affects the performance of the interferometer measurement. Our calculations are compared with the measurements performed on our cold atom gravimeter based on stimulated Raman transitions. We finally extend this study to the case of cold atom gradiometers.Comment: 17 pages, 6 figure

A Storage Ring for Neutral Atoms

We have demonstrated a storage ring for ultra-cold neutral atoms. Atoms with mean velocities of 1 m/s corresponding to kinetic energies of ~100 neV are confined to a 2 cm diameter ring by magnetic forces produced by two current-carrying wires. Up to 10^6 atoms are loaded at a time in the ring, and 7 revolutions are clearly observed. Additionally, we have demonstrated multiple loading of the ring and deterministic manipulation of the longitudinal velocity distribution of the atoms using applied laser pulses. Applications of this ring include large area atom interferometers and cw monochromatic atomic beam generation.Comment: 4 pages, 5 figure

Parallel computation of echelon forms

International audienceWe propose efficient parallel algorithms and implementations on shared memory architectures of LU factorization over a finite field. Compared to the corresponding numerical routines, we have identified three main difficulties specific to linear algebra over finite fields. First, the arithmetic complexity could be dominated by modular reductions. Therefore, it is mandatory to delay as much as possible these reductions while mixing fine-grain parallelizations of tiled iterative and recursive algorithms. Second, fast linear algebra variants, e.g., using Strassen-Winograd algorithm, never suffer from instability and can thus be widely used in cascade with the classical algorithms. There, trade-offs are to be made between size of blocks well suited to those fast variants or to load and communication balancing. Third, many applications over finite fields require the rank profile of the matrix (quite often rank deficient) rather than the solution to a linear system. It is thus important to design parallel algorithms that preserve and compute this rank profile. Moreover, as the rank profile is only discovered during the algorithm, block size has then to be dynamic. We propose and compare several block decomposition: tile iterative with left-looking, right-looking and Crout variants, slab and tile recursive. Experiments demonstrate that the tile recursive variant performs better and matches the performance of reference numerical software when no rank deficiency occur. Furthermore, even in the most heterogeneous case, namely when all pivot blocks are rank deficient, we show that it is possbile to maintain a high efficiency