Top Banner
Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes, Cédric Viou, Dalal Ait-Allal Plan : Pulsar Timing Instrumentations CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs
20

Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Sep 10, 2018

Download

Documents

truongdang
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Real-Time Pulsar TimingSignal Processing on GPUs

LPC2E, CNRS Orléans - FRANCE -Ismaël Cognard,Gilles Theureau,Grégory Desvignes,Cédric Viou,Dalal Ait-Allal

Plan :Pulsar TimingInstrumentations

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Page 2: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Pulsars

a highly magnetizedneutron star

producing two beamsacting like a cosmic lighthouse

pulsars arehighly stable clocks,suffering from dispersion

most fast and stable pulsars are usedto test Gravitation theories and

to probe the Cosmological Gravitational Wave background

Page 3: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Dispersion by the ISM

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

DM=3.139 pc.cm-3

---

1

28M

Hz

@1.

4GH

z

--

-

Page 4: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Pulsar timing

a huge computing power is requiredto coherent dedisperseover a large bandwith

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

(Lorimer & Kramer, 2004)

the coherent dedispersionremoves the dispersion effectof the interstellar mediumdirectly on recorded voltages

Page 5: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

a large radiotelescope : Nançay

a 100m dish equivalent telescope2 receivers 1.1-1.8 and 1.7-3.5GHzTsys 35K efficiency 1.4K.Jy

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Page 6: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

PSR J1909-3744 rms 109ns

PSR B1937+21 rms 387ns

PSR J1713+0747 rms 350ns

PSR J1600-3053 rms 576ns

Five stable pulsars timed at Nançay

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

PSR J1744-1134 rms 343ns

Page 7: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

a bit of history...

1988 swept LO (co 8MHz)~ 0.5μs uncertainty on PSR B1937+21

1996 filterbank NBPP (inco 120MHz)Navy Berkeley Pulsar Processor

2004 coherent dedispersor BON (co 64-128MHz)Berkeley Orleans Nancay

2008 GPU coherent BON (co 128MHz)

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

swept LO

NBPP

Page 8: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

BON : a coherent dedispersion with a cluster of CPUs 64 bi-Athlon 1.2GHzGigabit fiber network

Serendip5 installed in 2001

Page 9: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

in 2008,it was time to think about the upgrade of the 7 years old cluster

and of the 4 years old instrumentation...

following preliminary tests by Paul Demorest, UC Berkeley, in 2007,we tried different types of processors in addition to standard CPUs :

the GPUs with the Nvidia 8800GTX for example (G80 family)

Upgrade of the BON dedispersor

Page 10: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

so, we tested two types of processors in addition to standard CPUs :the IBM Cell (PlayStation3) and the Nvidia 8800GTX

Tests of different engines

PS3 is really goodbut very limited bythe only Gigabit link it has

(15-20MHz bandwidth)

Nvidia 8800GTX is fairly goodand less limited by the PCIe bus

(few Gbs easily)

Page 11: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

a GPU-based pulsarinstrumentation

developped at Nançay Jun

e 20

08

1 motherboardSupermicro X7DWA-N

2 CPUsquad-core Xeon

E5420 (2.5GHz, 12MB cache)4 GB main memory2 PCI-CDa DMA interface (for EDT)2 GeForce 8800GTXand water cooling

+ CUDA librairies + CUFFT

this 5k€ system is able toremove the dispersion effectof the interstellar mediumdirectly on recorded voltagesover a 128MHz bandwidth(the data rate is 2x 2Gbs)

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Page 12: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Architecture

Page 13: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

INPUT at 2Gb/s from Serendip5 spectrometer

OUTPUTon disk

100kB/s

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

the data path

Apr 2008, prototypewithout water cooling

Page 14: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

monitoring the system...

A/C failure

Page 15: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

oldAthlon cluster

newGPUs

probably still the only ones to do GPU dedispersion...

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Page 16: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

iBOB

800MHz clock

PC + 10GbE card (+ GPU)

10GbE Gore cable

Sept 2008

the next coherent pulsar dedispersor at Nançaywill have a 400MHz bandwith

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

iBOB tests conducted on Sept 2008 with P.McMahon

Page 17: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Early 2010 : a 400MHz coherent pulsar dedispersor

Acquisition : ROACH board (CASPER)Processing : 4 PCs with 8 GPUs (GT280)

~ 2012-2013direct sampling of the receivers outputs ?UNIBOARD collaboration with PRISME Univ Orléans (PhD)

(http://www.radionet-eu.org/uniboard)

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Page 18: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

the GT200 family

GT280 is 2.5-3 timesfaster than the 8800GTXfor dedispersion (FFT,chirp_mult,FFT-1)

Page 19: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs

Folding in the GPU

we developed kernel code todo the folding within the GPUand tested it...

Difference

we are now working on the re-orderingof data received from a PFB

(a matrix transposition)

CPU

GPU

Page 20: Real-Time Pulsar Timing Signal Processing on GPUs · Real-Time Pulsar Timing Signal Processing on GPUs LPC2E, CNRS Orléans - FRANCE - Ismaël Cognard, Gilles Theureau, Grégory Desvignes,

Conclusion

we guess GPUs are a good alternativeto do coherent dedispersion of pulsars

they are inexpensive at ~300€ each

the GT200 generation is able to processmore than 100MHz bw each

properly coded, GPUs can do all the job :re-ordering data, dedispersing, folding(maybe doing everything for bw 100MHz)

CASPER workshop, Sept 30, 2009 , Cape Town, SA Real-time Pulsar Timing signal processing on GPUs