CHPID: Calcul Hautes Performances et Informatique Distribuée

Table of Contents

1 [10/12/2009] GPU Computing, the numerical methods point of view

1.1 Invitation

You are invited to participate to the workshop on "GPU computing, the numerical methods point of view" that will be held on December 10th 2009, from 1 pm to 5pm in the room 316 of UFR IMAG on the campus of Grenoble.

The purpose of this workshop is to discuss the impact of GPU architectures in the development of numerical methods and software platforms in scientific computing. The intent is to present on going work for specific methods and then to discuss the benefit and drawback on a larger scale. The goal is to come out of the workshop with some answers for applied mathematicians on how GPU can be used to improve their codes. The workshop will mainly address the mathematical aspect of GPU computing.

During this workshop, four talks on GPU numerical computing will be given. It will conclude by a round table to discuss different aspects of GPU computing.

If you wish to participate to this workshop, please send an email to Christophe.Picard@imag.fr. The inscription is free, but we need to have a rough approximation on the number of people who will come to organize the coffee break.

Christophe Picard,
Christophe Prud'homme,
Laboratoire Jean Kuntzmann

1.2 Organization

This workshop is organized by C. Picard and C. Prud'homme with the support of the project CHPID of the research cluster ISLE.

  • Where: Grenoble, Université Joseph Fourier, Room 316 of UFR IMAG.
  • When: December 10, 2009 from 13:30 to 17:00

1.3 Talks

The presentations are about 45 min max, 30-35 min presentations and 5-10 min questions.

Time table for the workshop
TimeSpeaker
13:30-14:15Diego Rossinelli
14:15-15:00Vittoria Rezzonico
15:00-15:45Thomas GUIGNON & Jean-Marc Gratien
15:45-16:15Luigi Genovese
16:15-17:00Discussions
Round table

1.3.1 Diego Rossinelli

Email: diegor@inf.ethz.ch

  • Title

GPU accelerated simulations of bluff body flows using penalization and vortex methods

  • Abstract

We consider a state-of-the-art implementation of hybrid particle methods on GPUs. We present a GPU solver and a GPU-assisted solver for simulating bluff body flows using a remeshed vortex particle method. The efficiency of the method relies on fast and accurate particle-grid interpolations on GPUs for the remeshing of the particles and the computation of the field operators. The GPU implementation uses hardware framebuffer operations to perform efficient particle-grid operations and a CUFFT-based solver for the Poisson Equation with unbounded boundary conditions.

1.3.2 Vittoria Rezzonico

Email: Vittoria.Rezzonico@epfl.ch

  • Title

GPGPU @ EPFL: application showcase

  • Abstract

After a short presentation of the GPU architecture and of the available tools for efficiently programming GPUs, we will talk about several applications at EPFL which take advantage of the computational power of GPUs. These applications include Coarse-Grained potential for molecular simulation of proteins and multi-scale Hemodynamics.

  • Slides

The slides are here.

1.3.3 Thomas GUIGNON and Jean-Marc Gratien

Emails: Thomas.GUIGNON@ifp.fr, Jean-Marc.Gratien@ifp.fr

  • Title

Numerical methods on GPU to improve the efficiency of finite volume applications in Geosciences

Méthodes numériques sur GPU pour améliorer les performances des applications de types volumes finis en Géosciences.

  • Abstract

Afin d'optimiser et de déterminer le plan de développement d'un gisement pétrolier, les compagnies pétrolières ont besoin de prédire l'historique de production d'un gisement en fonction de divers scénarios d'implémentation des puits. Dans ce but, elles utilisent des simulateurs de réservoirs qui modélisent numériquement l'exploitation d'un champ pétrolier à partir des équations d'écoulement des fluides dans un milieu poreux. L'IFP poursuit des travaux de recherche dans le domaine de la simulation de réservoir dont la thématique s'inscrit dans une de ses priorités stratégiques visant à repousser les limites de l'exploration pétrolière.

La simulation de réservoir est basée sur des modèles d'écoulement en milieu poreux utilisant la loi de darcy résolus principalement avec des méthodes Volumes finis qui peuvent être multi échelles. Avec ce genre de méthode, l'essentiel du temps des application est réa parti dans : (1) l'utilisation intensive de solveurs itératifs creux de type biCGStab couplés à des préconditionneurs iLU ou multigrilles ; (2) la résolution de petits systèmes non linéaires locaux représentant les équilibres chimiques dans chaque maille (algorithme de flash) ; (3) la construction des données des maillages à l'échelle de la simulation à partir de données géophysiques obtenues à une échelle plus fine (algorithme d'Upscaling); (4) la construction de fonctions de base à l'échelle fine permettant de reconstruire les solutions échelle fine à partir de celles obtenues à une échelle grossière (méthode VF multi niveau).

L'IFP dans le cadre du projet ANR PARA (Parallelisme et Amélioration du Rendement des Applications) explore depuis 2 ans l'utilisation d'accélérateurs graphiques afin de déporter l'exécution de ce type d'algorithmes et ainsi réduire les temps d'exécution.

Nous présenterons l'état de nos travaux sur l'utilisation d'accélérateurs graphiques tout en présentant dans un premier temps les méthodes numériques utilisées, la méthodologie suivie afin d'implémenter le solveur creux efficacement sur un GPU puis nous présenterons les résultats obtenus.

  • Slides

The slides are here and here.

1.3.4 Luigi Genovese

Email: Luigi.Genovese@esrf.fr

  • Title

Efficient use of hybrid computing clusters for nanosciences

  • Abstract

In this talk, we will study the programming and the utilization of hybrid clusters in the field of computational physics. These massively parallel computers are composed of a fast network (Infiniband) connecting classical nodes with multicore Intel processors and accelerators. In our case, the accelerators used are GPUs from NVIDIA.

First, we will analyze some ways to use with efficiency CPUs cores and GPUs together in a code (BiGDFT, http://inac.cea.fr/L_Sim/BigDFT) without hotspot routines. Starting from this analyze, we have designed a new library: SGPU, used to share GPUs between the CPU cores of a node. The implementation and the usage of SGPU will be described in this talk.

Then, we will evaluate and compare performances between SGPU and others approaches to share GPUs with CPUs. This performance evaluation will be based on BigDFT, an ab-initio simulation software designed to take advantage of massively parallel clusters as the titane cluster (CCRT). Our experiments will be performed, on a first time, on one hybrid node, and on a second time, on a large number of nodes of the hybrid cluster.

  • Slides

The slides are here.

1.4 Roundtable

Here are a few question for the round table:

  • GPU computing : hype or hope?
  • What is the status of the development library? Is opencl mature enough?
  • How to integrate the various technologies? GPU, Clusters or Core?
  • Do we have to rethink the numerical approaches or do we leave the work to computer scientists?

1.5 Participants

PrenomNomAffiliationEmail
WalidBen RomdhaneUJF/LSPwalid.benromdhane@gmail.com
Dan-GabrielCalugaruUFC/LMBdan.calugaru@univ-fcomte.fr
Cerasela IlianaiCalugaruENS Lyoncerasela.iliana.calugaru@ens-lyon.fr
VincentChabannesUJF/LJKVincent.Chabannes@imag.fr
StéphaneChaillolUJF/LPGstephane.chaillol@obs.ujf-grenoble.fr
Georges-HenriCottetUJF/LJKGeorges-Henri.Cottet@imag.fr
ThomasDufaudUCBL/ICJThomas.Dufaud@cdcsp.univ-lyon1.fr
FrançoisFaureINRIAFrancois.Faure@imag.fr
PatrickFulconisUJF/LGITPatrick.Fulconis@ujf-grenoble.fr
LuigiGenoveseESRFLuigi.Genovese@esrf.fr
Pierre-YvesGiresUJF/LSPpgires@ujf-grenoble.fr
Jean-MarcGratienIFPJean-Marc.Gratien@ifp.fr
DavidGuibertUCBL/ICdavid.guibert@cdcsp.univ-lyon1.fr
ThomasGuignonIFPThomas.Guignon@ifp.fr
MarcHonnoratUJF/LJKmarc.honnorat@imag.fr
MouradIsmailUJFismail@ujf-grenoble.fr
StéphaneLabbéLJK/UJFStephane.labbe@imag.fr
FabienLe mentecIMAGfabien.le-mentec@imag.fr
PatriceLinelUCBL/ICJPatrice.Linel@cdcsp.univ-lyon1.fr
ViolaineLouvetUCBLlouvet@math.univ-lyon1.fr
AdrienMagniUJF/LJKadrien.magni@imag.fr
RolandMassonIFProland.masson@ifp.fr
NicolasMeilardUCBL/CRALnicolas.meilard@obs.univ-lyon1.fr
FabienMomeyUCBL/CRALfabien.momey@obs.univ-lyon1.fr
FranckPérignonLJK/UJFfranck.perignon@imag.fr
ChristophePicardINPG/LJKchristophe.picard@imag.fr
ChristophePrud'hommeUJF/LJKchristophe.prudhomme@imag.fr
EmmanuelQuémenerENS/LIPemmanuel.quemener@ens-lyon.fr
VittoriaRezzonicoEPFLVittoria.Rezzonico@epfl.ch
DiegoRossinelliETZHdiegor@inf.ethz.ch
GeraldSimeoniHPGerald.Simeoni@hp.com
DamienTromeur-DervoutUCBL/ICJDamien.Tromeur-Dervout@cdcsp.univ-lyon1.fr
MathiasVautrinINSERMvautrin@esrf.fr
FlorentVielfaureIMAGflorent.vielfaure@gmail.com
LaurenceViryUJF/LJKlaurence.Viry@imag.fr
WangYiqingINPG/LJKyqwang9@gmail.com

Author: E. Caron, C. Prud'homme

Date: <2009-10-19 Mon>

HTML generated by org-mode 7.3 in emacs 23