Optimising Sparse Matrix Vector multiplication for large scale FEM problems on FPGA

Sparse Matrix Vector multiplication (SpMV) is an important kernel in many scientific applications. In this work we propose an architecture and an automated customisation method to detect and optimise the architecture for block diagonal sparse matrices. We evaluate the proposed approach in the context of the spectral/hp Finite Element Method, using the local matrix assembly approach. This problem leads to a large sparse system of linear equations with block diagonal matrix which is typically solved using an iterative method such as the Preconditioned Conjugate Gradient. The efficiency of the proposed architecture combined with the effectiveness of the proposed customisation method reduces BRAM resource utilisation by as much as 10 times, while achieving identical throughput with existing state of the art designs and requiring minimal development effort from the end user. In the context of the Finite Element Method, our approach enables the solution of larger problems than previously possible, enabling the applicability of FPGAs to more interesting HPC problems

Relaminarisation of Re_τ=100 channel flow with globally stabilising linear feedback control

The problems of nonlinearity and high dimension have so far prevented a complete solution of the control of turbulent flow. Addressing the problem of nonlinearity, we propose a flow control strategy which ensures that the energy of any perturbation to the target profile decays monotonically. The controller’s estimate of the flow state is similarly guaranteed to converge to the true value. We present a one-time off-line synthesis procedure, which generalises to accommodate more restrictive actuation and sensing arrangements, with conditions for existence for the controller given in this case. The control is tested in turbulent channel flow (Re_τ = 100) using full-domain sensing and actuation on the wall-normal velocity. Concentrated at the point of maximum inflection in the mean profile, the control directly counters the supply of turbulence energy arising from the interaction of the wall-normal perturbations with the flow shear. It is found that the control is only required for the larger-scale motions, specifically those above the scale of the mean streak spacing. Minimal control effort is required once laminar flow is achieved. The response of the near-wall flow is examined in detail, with particular emphasis on the pressure and wall-normal velocity fields, in the context of Landahl’s theory of sheared turbulence

Entanglement distribution for a practical quantum-dot-based quantum processor architecture

We propose a quantum dot (QD) architecture for enabling universal quantum information processing. Quantum registers, consisting of arrays of vertically stacked self-assembled semiconductor QDs, are connected by chains of in-plane self-assembled dots. We propose an entanglement distributor, a device for producing and distributing maximally entangled qubits on demand, communicated through in-plane dot chains. This enables the transmission of entanglement to spatially separated register stacks, providing a resource for the realization of a sizeable quantum processor built from coupled register stacks of practical size. Our entanglement distributor could be integrated into many of the present proposals for self-assembled QD-based quantum computation (QC). Our device exploits the properties of simple, relatively short, spin-chains and does not require microcavities. Utilizing the properties of self-assembled QDs, after distribution the entanglement can be mapped into relatively long-lived spin qubits and purified, providing a flexible, distributed, off-line resource. © IOP Publishing Ltd and Deutsche Physikalische Gesellschaft

Time-Lapse Acoustic Imaging of Mesoscale and Fine-Scale Variability within the Faroe-Shetland Channel

We describe and analyze the results of a three‐dimensional seismic (i.e. acoustic) reflection survey from the Faroe‐Shetland Channel that is calibrated with near‐coincident hydrographic and satellite observations. 54 vertical seismic transects were acquired over a period of 25 days. On each transect, a 250‐‐400 m band of reflections is observed within the water column. Hydrographic measurements demonstrate that this reflective band is caused by temperature variations within the pycnocline that separates warm, near‐surface waters of Atlantic origin from cold, deep waters which flow southward from the Nordic Seas. Tilting of reflective surfaces records geostrophic shear between these near‐surface and deep waters. Measurements of temporal changes of pycnoclinic depth and of reflection tilt are used to infer the existence of an anticyclonic vortex that advects northeastward. Comparison with satellite measurements of sea‐surface temperature and height suggests that this vortex is caused by meandering of the Continental Slope Current. A model of a Gaussian vortex is used to match seismic and satellite observations. This putative vortex has a core radius of 20—30 km and a maximum azimuthal velocity of 0.3‐‐0.4 m s‐1. It translates at 0.01‐‐0.1 m s‐1. Within the pycnocline, diapycnal diffusivity, K , is estimaed by analyzing the turbulent spectral subrange of tracked reflections. K varies between 10‐5.7 and 10‐5.0 m 2 s‐1 in a pattern that is broadly consistent with translation of the vortex. Our integrated study demonstrates the ability of time‐lapse seismic reflection surveying to dynamically resolve the effects that mesoscale activity has upon deep thermohaline structure on scales from meters to hundreds of kilometers.Natural Environment Research Council (NERC) Engineering and Physical Science Research Council 794 Program Grant EP/K034529/

Breaking a Bloch-wave interferometer: quasiparticle species-specific temperature-dependent nonequilibrium dephasing

Recently, high-order sideband polarimetry has been established as an experimental method that links the polarization of sidebands to an interference of Bloch wavefunctions. However, the robustness of sideband polarizations to increasing dephasing remains to be explored. Here, we investigate the dependence of high-order sideband generation in bulk gallium arsenide on dephasing by tuning temperature. We find that the intensities of the sidebands, but not their polarizations, depend strongly on temperature. Using our polarimetry method, we are able to isolate the contributions of electron-heavy hole (HH) and electron-light hole (LH) pairs to sideband intensities, and separately extract the nonequilibrium dephasing coefficients associated with the longitudinal optical (LO) phonons and acoustic (A) phonons for each species of electron-hole pair. We find that $\Gamma_{\text{HH},\text{A}} = 6.1 \pm 1.6$ $\mu$eV/K, $\Gamma_{\text{LH},\text{A}} < 1.5$ $\mu$eV/K, $\Gamma_{\text{HH},\text{LO}} = 14 \pm 3$ meV, and $\Gamma_{\text{LH},\text{LO}} = 30 \pm 3$ meV.Comment: 10 pages, 4 figure