Gloria Ortega

Improving the Performance of the Sparse Matrix Vector Product with GPUs

Sparse matrices are involved in linear systems, eigensystems and partial differential equations from a wide spectrum of scientific and engineering disciplines. Hence, sparse matrix vector product (SpMV) is considered as key operation in engineering and scientific computing. For these applications the optimization of the sparse matrix vector product (SpMV) is very relevant. However, the irregular computation involved in SpMV prevents the optimum exploitation of computational architectures when the sparse matrices are very large. Graphics Processing Units (GPUs) have recently emerged as platforms that yield outstanding acceleration factors. SpMV implementations for GPUs have already appeared on the scene. This work proposes and evaluates new implementations of SpMV for GPUs called ELLR-T. They are based on the format ELLPACK-R, which allows storage of the sparse matrix in a regular manner. A comparative evaluation against a variety of storage formats previously proposed has been carried out based on a representative set of test matrices. The results show that: (1) the SpMV is highly accelerated with GPUs; (2) the performance strongly depends on the specific pattern of the matrix; and (3) the implementations ELLR-T achieve higher overall performance. Consequently, the new implementations of SpMV, ELLR-T, described in this paper can help to exploit the GPUs, because, they achieve high performance and they can be easily joined in the engineering and scientific computing.

Fuzzy Content-Based Image Retrieval for Oceanic Remote Sensing

The detection of mesoscale oceanic structures, such as upwellings or eddies, from satellite images has significance for marine environmental studies, coastal resource management, and ocean dynamics studies. Nevertheless, there is a lack of tools that allow us to retrieve automatically relevant mesoscale structures from large satellite image databases. This paper focuses on the development and validation of a content-based image retrieval system to classify and retrieve oceanic structures from satellite images. The images were obtained from the National Oceanic and Atmospheric Administration satellites Advanced Very High Resolution Radiometer sensor. The study area is about W2° - 21°, N19° - 45°. This system conducts labeling and retrieval of the most relevant and typical mesoscale oceanic structures, such as upwellings, eddies, and island wakes located in the Canary Islands area and in the Mediterranean and Cantabrian seas. Our work is based on several soft computing technologies such as fuzzy logic and neurofuzzy systems.

Performance evaluation of kernel fusion BLAS routines on the GPU: iterative solvers as case study

Programmers usually implement iterative methods that solve partial differential equations by expressing them using a sequence of basic kernels from libraries optimized for the graphics processing unit (GPU). The global runtime of the resulting combination is often penalized by the smallest and most inefficient vector operations. To improve the GPU exploitation, we identify and analyze the potential kernels to be fused according to the data dependence, data type and size, and GPU resources. This paper provides an extensive analysis of the impact of fusing vector operations [level 1 of Basic Linear Algebra Subprograms (BLAS)] on the performance of the GPU. The experimental evaluation shows that this optimization provides noticeable improvement especially for kernels with lower memory requirements and on more modern GPUs. It is worth noting that the fused BLAS operations can be very useful to help programmers efficiently code iterative methods to solve large linear systems of equations for the GPU. Iterative methods such as biconjugate gradient method (BCG) are one of the examples that can benefit from this optimization strategy. Indeed, kernel fusion of vector routines makes the most efficient GPU implementation of BCG run between

The BiConjugate gradient method on GPUs

1.09\times

Non-dominated sorting procedure for Pareto dominance ranking on multicore CPU and/or GPU

1.09× and

High performance computing for Optical Diffraction Tomography

1.27\times