Streaming Elements for FPGA Signal and Image Processing Accelerators

Field programmable gate array devices boast abundant resources with which custom accelerator components for signal, image and data processing may be realised; however, realising high performance, low cost accelerators currently demands manual register transfer level design. Software-programmable ’soft’ processors have been proposed as a way to reduce this design burden but they are unable to support performance and cost comparable to custom circuits. This paper proposes a new soft processing approach for FPGA which promises to overcome this barrier. A high performance, fine-grained streaming processor, known as a Streaming Accelerator Element, is proposed which realises accelerators as large scale custom multicore networks. By adopting a streaming execution approach with advanced program control and memory addressing capabilities, typical program inefficiencies can be almost completely eliminated to enable performance and cost which are unprecedented amongst software-programmable solutions. When used to realise accelerators for fast fourier transform, motion estimation, matrix multiplication and sobel edge detection it is shown how the proposed architecture enables real-time performance and with performance and cost comparable with hand-crafted custom circuit accelerators and up to two orders of magnitude beyond existing soft processors.

An R-matrix approach for plasma modelling and the interpretation of astrophysical observation

Over the last decade an Auburn-Rollins-Strathclyde consortium has developed several suites of parallel R-matrix codes [1, 2, 3] that can meet the fundamental data needs required for the interpretation of astrophysical observation and/or plasma experiments. Traditionally our collisional work on light fusion-related atoms has been focused towards spectroscopy and impurity transport for magnetically confined fusion devices. Our approach has been to provide a comprehensive data set for the excitation/ionization for every ion stage of a particular element. As we progress towards a burning fusion plasma, there is a demand for the collisional processes involving tungsten, which has required a revitalization of the relativistic R-matrix approach. The implementation of these codes on massively parallel supercomputers has facilitated the progression to models involving thousands of levels in the close-coupling expansion required by the open d and f sub-shell systems of mid Z tungsten. This work also complements the electron-impact excitation of Fe-Peak elements required by astrophysics, in particular the near neutral species, which offer similar atomic structure challenges. Although electron-impact excitation work is our primary focus in terms of fusion application, the single photon photoionisation codes are also being developed in tandem, and benefit greatly from this ongoing work.

Generalized collisional radiative model for light elements: C: Data for the B isonuclear sequence

A first stage collision database is assembled which contains electron-impact excitation, ionization,\r and recombination rate coefficients for B, B + , B 2+ , B 3+ , and B 4+ . The first stage database\r is constructed using the R-matrix with pseudostates, time-dependent close-coupling, and perturbative\r distorted-wave methods. A second stage collision database is then assembled which contains\r generalized collisional-radiative ionization, recombination, and power loss rate coefficients as a\r function of both temperature and density. The second stage database is constructed by solution of\r the collisional-radiative equations in the quasi-static equilibrium approximation using the first\r stage database. Both collision database stages reside in electronic form at the IAEA Labeled Atomic\r Data Interface (ALADDIN) database and the Atomic Data Analysis Structure (ADAS) open database.

Generalized collisional radiative model for light elements. Part B: Data for the Be isonuclear sequence

A first-stage collision database is assembled which contains electron-impact excitation, ionization, and recombination rate coefficients for Be, Be+, Be2+, and Be3+. The first-stage database is constructed using the R-matrix with pseudo-states, time-dependent close-coupling, and perturbative, distorted-wave methods. A second-stage collision database is then assembled which contains generalized collisional-radiative and radiated power loss coefficients. The second-stage database is constructed by solution of collisional-radiative equations in the quasi-static equilibrium approximation using the first-stage database. Both collision database stages reside in electronic form at the ORNL Controlled Fusion Atomic Data Center and in the ADAS database, and are easily accessed over the worldwide internet. © 2007 Elsevier Inc. All rights reserved.

Generalised collisional-radiative model for light elements. A: Data for the Li isonuclear sequence

A first stage collision database is assembled which contains electron-impact effective collision strengths, and ionization and recombination rate coefficients for Li, Li+, and Li2+. The first stage database is constructed using the R-matrix with pseudo-states, time-dependent close-coupling, converged close-coupling, and perturbative distorted-wave methods. A second stage collision database is then assembled which contains generalized collisional-radiative and radiated power loss coefficients. The second stage database is constructed by solution of collisional-radiative equations in the quasi-static equilibrium approximation using the first stage database. Both collision database stages reside in electronic form at the ORNL Controlled Fusion Atomic Data Center and in the ADAS database, and are easily accessed over the worldwide internet. ?? 2006 Elsevier Inc. All rights reserved.

Sparse matrix multiplication on a reconfigurable many-core architecture

Sparse matrix-vector multiplication (SMVM) is a fundamental operation in many scientific and engineering applications. In many cases sparse matrices have thousands of rows and columns where most of the entries are zero, while non-zero data is spread over the matrix. This sparsity of data locality reduces the effectiveness of data cache in general-purpose processors quite reducing their performance efficiency when compared to what is achieved with dense matrix multiplication. In this paper, we propose a parallel processing solution for SMVM in a many-core architecture. The architecture is tested with known benchmarks using a ZYNQ-7020 FPGA. The architecture is scalable in the number of core elements and limited only by the available memory bandwidth. It achieves performance efficiencies up to almost 70% and better performances than previous FPGA designs.

Investigations into the determinations of hydride-forming elements

Analytical methods for the determination of trace amounts of germanium, tin and arsenic were established using hydride generation coupled with direct current plasma atomic emission spectrometry. A continuous gas flowing batch system for the hydride generation was investigated and was applied to the determination of germanium(Ge), tin(Sn), antimony(Sb) and lead(Pb) (Preliminary results suggest that it is also applicable to arsenic)As) ). With this system, the reproducibility of signals was improved and the determination was speeded up, compared with the conventional batch type hydride generation system. Each determination was complete within one minute. Interferences from a number of transition metal ions, especially from Pd(II), Pt(IV), Ni(II), Cu(II), Co(II), and Fe(II, III), have proven to be very serious under normal conditions, in the determination of germanium, tin, and arsenic. These interference effects were eliminated or significantly reduced in the presence of L-cystine or L-cysteine. Thus, a 10-1000 fold excess of Ni(II), Cu(II), Co(II), Fe(II), Pt(IV), Pd(II), etc. can be tolerated without interference, In the presence of L-cystine or L-cysteine, compared with absence of interference reducing agent. The methods for the determination of Ge, Sn, and As were examined by the analyses of standard reference materials. Interference effects from the sample matrix, for example, in transition metal-rich samples, copper, iron and steel, were eliminated by L-cystine (for As and Sn) and by LI cysteine (for Ge). The analysis of a number of standard reference materials gave excellent results of As and Sn contents in agreement with the certified values, showing there was no systematic interference. The detection limits for both germanium and tin were 20 pg ml- I . Preliminary studies were carried out for the determination of antimony and lead. Antimony was found to react with NaBH4, remaInIng from the previous determinations, giving an analytical signal. A reversed injection manner, i.e., injection of the NaBH4 solution prior to the analyte solution was used to avoid uncertainty caused by residual NaBH4 present and to ensure that an excess of NaB H4 was available. A solution of 0.4% L-cysteine was found to reduce the interference from selected transition metal ions, Co(II), Cu(II), Ni(II) and Pt(IV). Hydrochloric acid - hydrogen peroxide, nitric acid - ammonium persulphate, and potassium dichromate malic acid reaction systems for lead hydride generation were compared. The potassium dichromate - malic acid reaction medium proved to be the best with respect to reproducibility and minimal interference. Cu(II), Ni(II), and Fe(II) caused strong interference In lead determinations, which was not reduced by L-cysteine or Lcystine. Sodium citrate, ascorbic acid, dithizone, thiosemicarbazide and penicillamine reduced interferences to some extent. Further interference reduction studies were carried out uSIng a number of amino acids, glycine, alanine, valine, leucine and histidine, as possible interference reducing agents in the determination of germanium. From glycine, alanine, valine to leucine, the interference reduction effect in germanium determinations decreased. Histidine II was found to be very promising In the reduction of interference. In fact, histidine proved more efficient than L-cystine and L-cysteine In the reduction of interference from Ni(II) in the determination of germanium. Signal enhancement by easily ionized elements (EIEs), usually regarded as an interference effect in analysis by DCP-AES, was studied and successfully applied to advantage in improving the sensitivity and detection limit in the determination of As, Ge, Sn, Sb, and Pb. The effect of alkali and alkaline-earth elements on the determination of the above five hydride forming elements was studied. With the appropriate EIE, a signal enhancement of 40-115% was achieved. Linear calibration and good reproducibility were also obtained in the presence of EIEs. III

Influence-matrix diagnostic of a data assimilation system

The influence matrix is used in ordinary least-squares applications for monitoring statistical multiple-regression analyses. Concepts related to the influence matrix provide diagnostics on the influence of individual data on the analysis - the analysis change that would occur by leaving one observation out, and the effective information content (degrees of freedom for signal) in any sub-set of the analysed data. In this paper, the corresponding concepts have been derived in the context of linear statistical data assimilation in numerical weather prediction. An approximate method to compute the diagonal elements of the influence matrix (the self-sensitivities) has been developed for a large-dimension variational data assimilation system (the four-dimensional variational system of the European Centre for Medium-Range Weather Forecasts). Results show that, in the boreal spring 2003 operational system, 15% of the global influence is due to the assimilated observations in any one analysis, and the complementary 85% is the influence of the prior (background) information, a short-range forecast containing information from earlier assimilated observations. About 25% of the observational information is currently provided by surface-based observing systems, and 75% by satellite systems. Low-influence data points usually occur in data-rich areas, while high-influence data points are in data-sparse areas or in dynamically active regions. Background-error correlations also play an important role: high correlation diminishes the observation influence and amplifies the importance of the surrounding real and pseudo observations (prior information in observation space). Incorrect specifications of background and observation-error covariance matrices can be identified, interpreted and better understood by the use of influence-matrix diagnostics for the variety of observation types and observed variables used in the data assimilation system. Copyright © 2004 Royal Meteorological Society

Advanced source apportionment of size-resolved trace elements at multiple sites in London during winter

Trace element measurements in PM10–2.5, PM2.5–1.0 and PM1.0–0.3 aerosol were performed with 2 h time resolution at kerbside, urban background and rural sites during the ClearfLo winter 2012 campaign in London. The environment-dependent variability of emissions was characterized using the Multilinear Engine implementation of the positive matrix factorization model, conducted on data sets comprising all three sites but segregated by size. Combining the sites enabled separation of sources with high temporal covariance but significant spatial variability. Separation of sizes improved source resolution by preventing sources occurring in only a single size fraction from having too small a contribution for the model to resolve. Anchor profiles were retrieved internally by analysing data subsets, and these profiles were used in the analyses of the complete data sets of all sites for enhanced source apportionment. A total of nine different factors were resolved (notable elements in brackets): in PM10–2.5, brake wear (Cu, Zr, Sb, Ba), other traffic-related (Fe), resuspended dust (Si, Ca), sea/road salt (Cl), aged sea salt (Na, Mg) and industrial (Cr, Ni); in PM2.5–1.0, brake wear, other traffic-related, resuspended dust, sea/road salt, aged sea salt and S-rich (S); and in PM1.0–0.3, traffic-related (Fe, Cu, Zr, Sb, Ba), resuspended dust, sea/road salt, aged sea salt, reacted Cl (Cl), S-rich and solid fuel (K, Pb). Human activities enhance the kerb-to-rural concentration gradients of coarse aged sea salt, typically considered to have a natural source, by 1.7–2.2. These site-dependent concentration differences reflect the effect of local resuspension processes in London. The anthropogenically influenced factors traffic (brake wear and other traffic-related processes), dust and sea/road salt provide further kerb-to-rural concentration enhancements by direct source emissions by a factor of 3.5–12.7. The traffic and dust factors are mainly emitted in PM10–2.5 and show strong diurnal variations with concentrations up to 4 times higher during rush hour than during night-time. Regionally influenced S-rich and solid fuel factors, occurring primarily in PM1.0–0.3, have negligible resuspension influences, and concentrations are similar throughout the day and across the regions.

The compositional and configurational heterogeneity of matrix habitats shape woodland carabid communities in wooded-agricultural landscapes

Landscape heterogeneity (the composition and configuration of matrix habitats) plays a major role in shaping species communities in wooded-agricultural landscapes. However, few studies consider the influence of different types of semi-natural and linear habitats in the matrix, despite their known ecological value for biodiversity. Objective To investigate the importance of the composition and configuration of matrix habitats for woodland carabid communities and identify whether specific landscape features can help to maintain long-term populations in wooded-agricultural environments. Methods Carabids were sampled from woodlands in 36 tetrads of 4 km2 across southern Britain. Landscape heterogeneity including an innovative representation of linear habitats was quantified for each tetrad. Carabid community response was analysed using ordination methods combined with variation partitioning and additional response trait analyses. Results Woodland carabid community response was trait-specific and better explained by simultaneously considering the composition and configuration of matrix habitats. Semi-natural and linear features provided significant refuge habitat and functional connectivity. Mature hedgerows were essential for slow-dispersing carabids in fragmented landscapes. Species commonly associated with heathland were correlated with inland water and woodland patches despite widespread heathland conversion to agricultural land, suggesting that species may persist for some decades when elements representative of the original habitat are retained following landscape modification. Conclusions Semi-natural and linear habitats have high biodiversity value. Landowners should identify features that can provide additional resources or functional connectivity for species relative to other habitat types in the landscape matrix. Agri-environment options should consider landscape heterogeneity to identify the most efficacious changes for biodiversity.

A Time Efficient Optical Model for GATE Simulation of a LYSO Scintillation Matrix Used in PET Applications

A time efficient optical model is proposed for GATE simulation of a LYSO scintillation matrix coupled to a photomultiplier. The purpose is to avoid the excessively long computation time when activating the optical processes in GATE. The usefulness of the model is demonstrated by comparing the simulated and experimental energy spectra obtained with the dual planar head equipment for dosimetry with a positron emission tomograph ( DoPET). The procedure to apply the model is divided in two steps. Firstly, a simplified simulation of a single crystal element of DoPET is used to fit an analytic function that models the optical attenuation inside the crystal. In a second step, the model is employed to calculate the influence of this attenuation in the energy registered by the tomograph. The use of the proposed optical model is around three orders of magnitude faster than a GATE simulation with optical processes enabled. A good agreement was found between the experimental and simulated data using the optical model. The results indicate that optical interactions inside the crystal elements play an important role on the energy resolution and induce a considerable degradation of the spectra information acquired by DoPET. Finally, the same approach employed by the proposed optical model could be useful to simulate a scintillation matrix coupled to a photomultiplier using single or dual readout scheme.

Microgranitic Enclaves as Products of Self-mixing Events: a Study of Open-system Processes in the Maua Granite, Sao Paulo, Brazil, Based on in situ Isotopic and Trace Elements in Plagioclase

Felsic microgranular enclaves with structures indicating that they interacted in a plastic state with their chemically similar host granite are abundant in the Maua Pluton, SE Brazil. Larger plagioclase xenocrysts are in textural disequilibrium with the enclave groundmass and show complex zoning patterns with partially resorbed An-rich cores (locally with patchy textures) surrounded by more sodic rims. In situ laser ablation-(multi-collector) inductively coupled plasma mass spectrometry trace element and Sr isotopic analyses performed on the plagioclase xenocrysts indicate open-system crystallization; however, no evidence of derivation from more primitive basic melts is observed. The An-rich cores have more radiogenic initial Sr isotopic ratios that decrease towards the outermost part of the rims, which are in isotopic equilibrium with the matrix plagioclase. These profiles may have been produced by either (1) diffusional re-equilibration after rim crystallization from the enclave-forming magma, as indicated by relatively short calculated residence times, or (2) episodic contamination with a decrease of the contaminant ratio proportional to the extent to which the country rocks were isolated by the crystallization front. Profiles of trace elements with high diffusion coefficients would require unrealistically long residence times, and can be modeled in terms of fractional crystallization. A combination of trace element and Sr isotope data suggests that the felsic microgranular enclaves from the Maua Pluton are the products of interaction between end-member magmas that had similar compositions, thus recording `self-mixing` events.

Matrix Representation of the Stationary Measure for the Multispecies TASEP

In this work we construct the stationary measure of the N species totally asymmetric simple exclusion process in a matrix product formulation. We make the connection between the matrix product formulation and the queueing theory picture of Ferrari and Martin. In particular, in the standard representation, the matrices act on the space of queue lengths. For N > 2 the matrices in fact become tensor products of elements of quadratic algebras. This enables us to give a purely algebraic proof of the stationary measure which we present for N=3.

Experimental damage detection grammian matrix

Structural Health Monitoring (SHM) has diverse potential applications, and many groups work in the development of tools and techniques for monitoring structural performance. These systems use arrays of sensors and can be integrated with remote or local computers. There are several different approaches that can be used to obtain information about the existence, location and extension of faults by non destructive tests. In this paper an experimental technique is proposed for damage location based on an observability grammian matrix. The dynamic properties of the structure are identified through experimental data using the eigensystem realization algorithm (ERA). Experimental tests were carried out in a structure through varying the mass of some elements. Output signals were obtained using accelerometers.

Prostate carcinogenesis induced by N-methyl-N-nitrosourea (mnu) in gerbils: Histopathological diagnosis and potential invasiveness mediated by extracellular matrix components

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)