MSGR-based Low Latency Complex Matrix Inversion Architecture

This paper presents a matrix inversion architecture based on the novel Modified Squared Givens Rotations (MSGR) algorithm, which extends the original SGR method to complex valued data, and also corrects erroneous results in the original SGR method when zeros occur on the diagonal of the matrix either initially or during processing. The MSGR algorithm also avoids complex dividers in the matrix inversion, thus minimising the complexity of potential real-time implementations. A systolic array architecture is implemented and FPGA synthesis results indicate a high-throughput low-latency complex matrix inversion solution. © 2008 IEEE.

Keratoconus: an analysis of corneal asymmetry

Background: Keratoconus, a non-inflammatory corneal ectasia, is reported to have bilateral involvement in over 90% of patients. The purpose of this study was to quantify the extent of asymmetry of disease at presentation to a regional corneal clinic.

Electroconvection in a homeotropic bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We characterize the structural transitions in an initially homeotropic bent-rod nematic liquid crystal excited by ac fields of frequency f well above the dielectric inversion point f(i). From the measured principal dielectric constants and electrical conductivities of the compound, the Carr-Helfrich conduction regime is anticipated to extend into the sub-megahertz region. Periodic patterned states occur through secondary bifurcations from the Freedericksz distorted state. An anchoring transition between the bend Freedericksz (1317) and degenerate planar (DP) states is detected. The BF state is metastable well above the Freedericksz threshold and gives way to the DP state, which persists in the field-off condition for several hours. Numerous +1 and -1 umbilics form at the onset of BF distortion, the former being largely of the chiral type. They survive in the DP configuration as linear defects, nonsingular in the core. In the BF regime, not far from fi, periodic Williams-like domains form around the umbilics; they drift along the director easy axis right from their onset. With increasing f, the wave vector of the periodic domains switches from parallel to normal disposition with respect to the c vector. Well above fi, a broadband instability is found.

Drifting periodic structures in a degenerate-planar bent-rod nematic liquid crystal beyond the dielectric inversion frequency

We report on the electric-field-generated effects in the nematic phase of a twin mesogen formed of bent-core and calamitic units, aligned homeotropically in the initial ground state and examined beyond the dielectric inversion point. The bend-Freedericksz (BF) state occurring at the primary bifurcation and containing a network of umbilics is metastable; we focus here on the degenerate planar (DP) configuration that establishes itself at the expense of the BF state in the course of an anchoring transition. In the DP regime, normal rolls, broad domains, and chevrons (both defect-mediated and defect-free types) form at various linear defect-sites, in different regions of the frequency-voltage plane. A significant novel aspect common to all these patterned states is the sustained propagative instability, which does not seem explicable on the basis of known driving mechanisms.

QR decomposition-based matrix inversion for high performance embedded MIMO receivers

Real-time matrix inversion is a key enabling technology in multiple-input multiple-output (MIMO) communications systems, such as 802.11n. To date, however, no matrix inversion implementation has been devised which supports real-time operation for these standards. In this paper, we overcome this barrier by presenting a novel matrix inversion algorithm which is ideally suited to high performance floating-point implementation. We show how the resulting architecture offers fundamentally higher performance than currently published matrix inversion approaches and we use it to create the first reported architecture capable of supporting real-time 802.11n operation. Specifically, we present a matrix inversion approach based on modified squared Givens rotations (MSGR). This is a new QR decomposition algorithm which overcomes critical limitations in other QR algorithms that prohibits their application to MIMO systems. In addition, we present a novel modification that further reduces the complexity of MSGR by almost 20%. This enables real-time implementation with negligible reduction in the accuracy of the inversion operation, or the BER of a MIMO receiver based on this.

Anisotropic seismic inversion using a multigrid Monte Carlo approach

We propose a new approach for the inversion of anisotropic P-wave data based on Monte Carlo methods combined with a multigrid approach. Simulated annealing facilitates objective minimization of the functional characterizing the misfit between observed and predicted traveltimes, as controlled by the Thomsen anisotropy parameters (epsilon, delta). Cycling between finer and coarser grids enhances the computational efficiency of the inversion process, thus accelerating the convergence of the solution while acting as a regularization technique of the inverse problem. Multigrid perturbation samples the probability density function without the requirements for the user to adjust tuning parameters. This increases the probability that the preferred global, rather than a poor local, minimum is attained. Undertaking multigrid refinement and Monte Carlo search in parallel produces more robust convergence than does the initially more intuitive approach of completing them sequentially. We demonstrate the usefulness of the new multigrid Monte Carlo (MGMC) scheme by applying it to (a) synthetic, noise-contaminated data reflecting an isotropic subsurface of constant slowness, horizontally layered geologic media and discrete subsurface anomalies; and (b) a crosshole seismic data set acquired by previous authors at the Reskajeage test site in Cornwall, UK. Inverted distributions of slowness (s) and the Thomson anisotropy parameters (epsilon, delta) compare favourably with those obtained previously using a popular matrix-based method. Reconstruction of the Thomsen epsilon parameter is particularly robust compared to that of slowness and the Thomsen delta parameter, even in the face of complex subsurface anomalies. The Thomsen epsilon and delta parameters have enhanced sensitivities to bulk-fabric and fracture-based anisotropies in the TI medium at Reskajeage. Because reconstruction of slowness (s) is intimately linked to that epsilon and delta in the MGMC scheme, inverted images of phase velocity reflect the integrated effects of these two modes of anisotropy. The new MGMC technique thus promises to facilitate rapid inversion of crosshole P-wave data for seismic slownesses and the Thomsen anisotropy parameters, with minimal user input in the inversion process.

Signatures of delayed detonation, asymmetry, and electron capture in the mid-infrared spectra of supernovae 2003hv and 2005df

We present mid-infrared (5.2-15.2 mu m) spectra of the Type Ia supernovae (SNe Ia) 2003hv and 2005df observed with the Spitzer Space Telescope. These are the first observed mid-infrared spectra of thermonuclear supernovae, and show strong emission from fine-structure lines of Ni, Co, S, and Ar. The detection of Ni emission in SN 2005df 135 days after the explosion provides direct observational evidence of high-density nuclear burning forming a significant amount of stable Ni in a SN Ia. The SN 2005df Ar lines also exhibit a two-pronged emission profile, implying that the Ar emission deviates significantly from spherical symmetry. The spectrum of SN 2003hv also shows signs of asymmetry, exhibiting blueshifted [Co (III)], which matches the blueshift of [Fe (II)] lines in nearly coeval near-infrared spectra. Finally, local thermodynamic equilibrium abundance estimates for the yield of radioactive Ni-56 give M-56Ni approximate to 0.5 M-circle dot, for SN 2003hv, but only M-56Ni approximate to 0.13-0.22 M-circle dot for the apparently subluminous SN 2005df, supporting the notion that the luminosity of SNe Ia is primarily a function of the radioactive 56Ni yield. The observed emission-line profiles in the SN 2005df spectrum indicate a chemically stratified ejecta structure, which matches the predictions of delayed detonation (DD) models, but is entirely incompatible with current three-dimensional deflagration models. Furthermore, the degree that this layering persists to the innermost regions of the supernova is difficult to explain even in a DD scenario, where the innermost ejecta are still the product of deflagration burning. Thus, while these results are roughly consistent with a delayed detonation, it is clear that a key piece of physics is still missing from our understanding of the earliest phases of SN Ia explosions.

Possible evidence of asymmetry in SN 2007rt, a type IIn supernova

An optical photometric and spectroscopic analysis of the slowly-evolving type IIn SN 2007rt is presented, covering a duration of 481 days after discovery. Its earliest spectrum, taken approximately 100 days after the explosion epoch, indicates the presence of a dense circumstellar medium, with which the supernova ejecta is interacting. This is supported by the slowly-evolving light curve. A notable feature in the spectrum of SN 2007rt is the presence of a broad He i 5875 line, not usually detected in type IIn supernovae. This may imply that the progenitor star has a high He/H ratio, having shed a significant portion of its hydrogen shell via mass-loss. An intermediate resolution spectrum reveals a narrow Ha P-Cygni profile, the absorption component of which has a width of 128 km s-1. This slow velocity suggests that the progenitor of SN 2007rt recently underwent mass-loss with wind speeds comparable to the lower limits of those detected in luminous blue variables. Asymmetries in the line profiles of H and He at early phases bears some resemblance to double-peaked features observed in a number of Ib/c spectra. These asymmetries may be indicative of an asymmetric or bipolar outflow or alternatively dust formation in the fast expanding ejecta. In addition, the late time spectrum, at over 240 days post-explosion, shows clear evidence for the presence of newly formed dust.

A high C/O ratio and weak thermal inversion in the atmosphere of exoplanet WASP-12b

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior, as opposed to the silicate-dominated composition found on Earth; the atmosphere can also differ from those in the Solar System. The solar C/O is 0.54 (ref. 3). Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b (ref. 6) that reveals C/O>=1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapour and enhanced in methane, each by more than two orders of magnitude compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T>2,500K) of WASP-12b lacks a prominent thermal inversion (or stratosphere) and has very efficient day-night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.