A debugging engine for parallel and distributed programs

Dissertação apresentada para a obtenção do Grau de Doutor em Informática pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia.

Parallel hyperspectral compressive sensing method on GPU

Remote hyperspectral sensors collect large amounts of data per flight usually with low spatial resolution. It is known that the bandwidth connection between the satellite/airborne platform and the ground station is reduced, thus a compression onboard method is desirable to reduce the amount of data to be transmitted. This paper presents a parallel implementation of an compressive sensing method, called parallel hyperspectral coded aperture (P-HYCA), for graphics processing units (GPU) using the compute unified device architecture (CUDA). This method takes into account two main properties of hyperspectral dataset, namely the high correlation existing among the spectral bands and the generally low number of endmembers needed to explain the data, which largely reduces the number of measurements necessary to correctly reconstruct the original data. Experimental results conducted using synthetic and real hyperspectral datasets on two different GPU architectures by NVIDIA: GeForce GTX 590 and GeForce GTX TITAN, reveal that the use of GPUs can provide real-time compressive sensing performance. The achieved speedup is up to 20 times when compared with the processing time of HYCA running on one core of the Intel i7-2600 CPU (3.4GHz), with 16 Gbyte memory.

Parallel hyperspectral coded aperture for compressive sensing on GPUs

The application of compressive sensing (CS) to hyperspectral images is an active area of research over the past few years, both in terms of the hardware and the signal processing algorithms. However, CS algorithms can be computationally very expensive due to the extremely large volumes of data collected by imaging spectrometers, a fact that compromises their use in applications under real-time constraints. This paper proposes four efficient implementations of hyperspectral coded aperture (HYCA) for CS, two of them termed P-HYCA and P-HYCA-FAST and two additional implementations for its constrained version (CHYCA), termed P-CHYCA and P-CHYCA-FAST on commodity graphics processing units (GPUs). HYCA algorithm exploits the high correlation existing among the spectral bands of the hyperspectral data sets and the generally low number of endmembers needed to explain the data, which largely reduces the number of measurements necessary to correctly reconstruct the original data. The proposed P-HYCA and P-CHYCA implementations have been developed using the compute unified device architecture (CUDA) and the cuFFT library. Moreover, this library has been replaced by a fast iterative method in the P-HYCA-FAST and P-CHYCA-FAST implementations that leads to very significant speedup factors in order to achieve real-time requirements. The proposed algorithms are evaluated not only in terms of reconstruction error for different compressions ratios but also in terms of computational performance using two different GPU architectures by NVIDIA: 1) GeForce GTX 590; and 2) GeForce GTX TITAN. Experiments are conducted using both simulated and real data revealing considerable acceleration factors and obtaining good results in the task of compressing remotely sensed hyperspectral data sets.

Parallel hyperspectral unmixing method via split augmented lagrangian on GPU

One of the main problems of hyperspectral data analysis is the presence of mixed pixels due to the low spatial resolution of such images. Linear spectral unmixing aims at inferring pure spectral signatures and their fractions at each pixel of the scene. The huge data volumes acquired by hyperspectral sensors put stringent requirements on processing and unmixing methods. This letter proposes an efficient implementation of the method called simplex identification via split augmented Lagrangian (SISAL) which exploits the graphics processing unit (GPU) architecture at low level using Compute Unified Device Architecture. SISAL aims to identify the endmembers of a scene, i.e., is able to unmix hyperspectral data sets in which the pure pixel assumption is violated. The proposed implementation is performed in a pixel-by-pixel fashion using coalesced accesses to memory and exploiting shared memory to store temporary data. Furthermore, the kernels have been optimized to minimize the threads divergence, therefore achieving high GPU occupancy. The experimental results obtained for the simulated and real hyperspectral data sets reveal speedups up to 49 times, which demonstrates that the GPU implementation can significantly accelerate the method's execution over big data sets while maintaining the methods accuracy.

A study on parallel versus sequential relational fuzzy clustering methods

Dissertação para obtenção do Grau de Mestre em Engenharia Informática

Optimization of breast tomosynthesis image reconstruction using parallel computing

Breast cancer is the most common cancer among women, being a major public health problem. Worldwide, X-ray mammography is the current gold-standard for medical imaging of breast cancer. However, it has associated some well-known limitations. The false-negative rates, up to 66% in symptomatic women, and the false-positive rates, up to 60%, are a continued source of concern and debate. These drawbacks prompt the development of other imaging techniques for breast cancer detection, in which Digital Breast Tomosynthesis (DBT) is included. DBT is a 3D radiographic technique that reduces the obscuring effect of tissue overlap and appears to address both issues of false-negative and false-positive rates. The 3D images in DBT are only achieved through image reconstruction methods. These methods play an important role in a clinical setting since there is a need to implement a reconstruction process that is both accurate and fast. This dissertation deals with the optimization of iterative algorithms, with parallel computing through an implementation on Graphics Processing Units (GPUs) to make the 3D reconstruction faster using Compute Unified Device Architecture (CUDA). Iterative algorithms have shown to produce the highest quality DBT images, but since they are computationally intensive, their clinical use is currently rejected. These algorithms have the potential to reduce patient dose in DBT scans. A method of integrating CUDA in Interactive Data Language (IDL) is proposed in order to accelerate the DBT image reconstructions. This method has never been attempted before for DBT. In this work the system matrix calculation, the most computationally expensive part of iterative algorithms, is accelerated. A speedup of 1.6 is achieved proving the fact that GPUs can accelerate the IDL implementation.

Para cá e para lá dos muros: relações familiares na interface entre o interior e o exterior da prisão

Tese de Doutoramento em Sociologia

An illustrated flora of the northern United States, Canada and the British possessions : from Newfoundland to the parallel of the southern boundary of Virginia, and from the Atlantic Ocean westward to the 102d meridian / by Nathaniel Lord Britton and Hon. Addison Brown.

3

An illustrated flora of the northern United States, Canada and the British possessions : from Newfoundland to the parallel of the southern boundary of Virginia, and from the Atlantic Ocean westward to the 102d meridian / by Nathaniel Lord Britton and Hon. Addison Brown.

2

An illustrated flora of the northern United States, Canada and the British possessions : from Newfoundland to the parallel of the southern boundary of Virginia, and from the Atlantic Ocean westward to the 102d meridian / by Nathaniel Lord Britton and Hon. Addison Brown.

1

An illustrated flora of the northern United States, Canada and the British possessions, from Newfoundland to the parallel of the southern boundary of Virginia, and from the Atlantic Ocean westward to the 102d meridian, by Nathaniel Lord Britton and Addison Brown. The descriptive text chiefly prepared by Britton, with the assistance of specialists in several groups; the figures also drawn under his supervision.

2

An illustrated flora of the northern United States, Canada and the British possessions, from Newfoundland to the parallel of the southern boundary of Virginia, and from the Atlantic Ocean westward to the 102d meridian, by Nathaniel Lord Britton and Addison Brown. The descriptive text chiefly prepared by Britton, with the assistance of specialists in several groups; the figures also drawn under his supervision.

3

An atlas of human gene expression from massively parallel signature sequencing (MPSS).

We have used massively parallel signature sequencing (MPSS) to sample the transcriptomes of 32 normal human tissues to an unprecedented depth, thus documenting the patterns of expression of almost 20,000 genes with high sensitivity and specificity. The data confirm the widely held belief that differences in gene expression between cell and tissue types are largely determined by transcripts derived from a limited number of tissue-specific genes, rather than by combinations of more promiscuously expressed genes. Expression of a little more than half of all known human genes seems to account for both the common requirements and the specific functions of the tissues sampled. A classification of tissues based on patterns of gene expression largely reproduces classifications based on anatomical and biochemical properties. The unbiased sampling of the human transcriptome achieved by MPSS supports the idea that most human genes have been mapped, if not functionally characterized. This data set should prove useful for the identification of tissue-specific genes, for the study of global changes induced by pathological conditions, and for the definition of a minimal set of genes necessary for basic cell maintenance. The data are available on the Web at http://mpss.licr.org and http://sgb.lynxgen.com.

Palaeoenvironment and palaeoclimate of the Middle Miocene lake in the Steinheim basin, SW Germany: A reconstruction from C, O, and Sr isotopes of fossil remains

A differentiated reconstruction of palaeolimnologic, -environmental, and -climatic conditions is presented for the Middle Miocene long-term freshwater lake (14.3 to 13.5 Ma) of the Steinheim basin, on the basis of a combined C, 0, and Sr isotope study of sympatric skeletal fossils of aquatic and terrestrial organisms from the lake sediments. The oxygen isotope composition for lake water of the Steinheim basin (delta O-18(H2O) = +2.0 +/- 0.4 parts per thousand VSMOW, n = 6) was reconstructed from measurements of delta O-18(PO4) of aquatic turtle bones. The drinking water calculated from the enamel of large mammals (proboscideans, rhinocerotids, equids, cervids, suids) has delta O-18(H2O) values (delta(OH2O)-O-18 = -5.9 +/- 1.7 parts per thousand VSMOW, n = 31) typical for Middle Miocene meteoric water of the area. This delta O-18(H2O) value corresponds to a mean annual air temperature (MAT) of 18.8 +/- 3.8 degrees C, calculated using a modem-day delta(OH2O)-O-18-MAT relation. Hence, large mammals did not use the lake water as principal drinking water. In contrast, small mammals, especially the then abundant pika Prolagus oeningensis drank from O-18-enriched water sources (delta O-18(H2O) = +2.7 +/- 2.3 parts per thousand VSMOW, n = 7), such as the lake water. Differences in Sr and 0 isotopic compositions between large and small mammal teeth indicate different home ranges and drinking behaviour and support migration of some large mammals between the Swabian Alb plateau and the nearby Molasse basin, while small mammals ingested their food and water locally. Changes in the lake level, water chemistry, and temperature were inferred using isotopic compositions of ostracod and gastropod shells from a composite lake sediment profile. Calcitic ostracod valves (Ilyocypris binocularis; delta O-18 = +1.7 +/- 1.2 parts per thousand VPDB, delta C-18 = -0.5 +/- 0.9 parts per thousand, VPDB, n = 68) and aragonitic, gastropod shells (Gyraulus spp.; delta O-18 = +2.0 +/- 13 parts per thousand VPDB, delta C-13 = -1.1 +/- 1.3 parts per thousand VPDB, n = 89) have delta O-18 and delta C-13 values similar to or even higher than those of marine, carbonates. delta C-13 values:of the biogenic carbonates parallel lake level fluctuations while delta O-18 values scatter around +2 +/- 2 parts per thousand and reflect the short term variability of meteoric water inflow vs. longer term evaporation. Sr-87/Sr-86 ratios of aragonitic Gyraulus spp. gastropod shells parallel the lake level fluctuations, reflecting variable inputs of groundwater and surface waters. Using a water delta O-18(H2O) value of +2.0 parts per thousand VSMOW, water temperatures calculated from skeletal tissue delta O-18 values of ostracods are 16.7 +/- 5.0 degrees C, gastropods 20.6 +/- 5.6 degrees C, otoliths 21.8 +/- 1.4 degrees C, and fish teeth 17.0 +/- 2.7 degrees C. The calculated MAT (similar to 19 degrees C), lake water temperatures (similar to 17 to 22 degrees C) and the O-18-enriched water compositions are indicative of warm-temperate climatic conditions, possibly with a high humidity during this period. Vegetation in the area surrounding the basin was largely of the C-3-type, as indicated by carbon isotopic compositions of tooth enamel from large mammals (delta C-13 = -11.1 +/- 1.1 parts per thousand VPDB, n = 40). (c) 2006 Elsevier B.V. All rights reserved.

Evolution of C(4) phosphoenolpyruvate carboxykinase in grasses, from genotype to phenotype.

C(4) photosynthesis is an adaptation over the classical C(3) pathway that has evolved multiple times independently. These convergences are accompanied by strong variations among the independent C(4) lineages. The decarboxylating enzyme used to release CO(2) around Rubisco particularly differs between C(4) species, a criterion used to distinguish three distinct biochemical C(4) subtypes. The phosphoenolpyruvate carboxykinase (PCK) serves as a primary decarboxylase in a minority of C(4) species. This enzyme is also present in C(3) plants, where it is responsible for nonphotosynthetic functions. The genetic changes responsible for the evolution of C(4)-specific PCK are still unidentified. Using phylogenetic analyses on PCK sequences isolated from C(3) and C(4) grasses, this study aimed at resolving the evolutionary history of C(4)-specific PCK enzymes. Four independent evolutions of C(4)-PCK were shown to be driven by positive selection, and nine C(4)-adaptive sites underwent parallel genetic changes in different C(4) lineages. C(4)-adaptive residues were also observed in C(4) species from the nicotinamide adenine dinucleotide phosphate-malic enzyme (NADP-ME) subtype and particularly in all taxa where a PCK shuttle was previously suggested to complement the NADP-ME pathway. Acquisitions of C(4)-specific PCKs were mapped on a species tree, which revealed that the PCK subtype probably appeared at the base of the Chloridoideae subfamily and was then recurrently lost and secondarily reacquired at least three times. Linking the genotype to subtype phenotype shed new lights on the evolutionary transitions between the different C(4) subtypes.