DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
•
•
•
•
Institutions: DESY,Zeuthen, Bielefeld Univ., Helmholtz-Gemeinschaft, EPS
Local organizing committee: Johannes Blümlein, Jochem Fleischer, Wolfgang Friebel,
Ulrich Gensch, Matthias Kasemann, Thomas Naumann, Tord Riemann (Chair),
Peter Wegner
International Advisory Committee members
Conference Secretaries: Martina Mende, Christiane John, Sabine Miethling
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
AIHENP: “Artificial Intelligence in High Energy and Nuclear Physics”
• Lyon 1990, la Londes-Les-Maures 1992, Oberammergau 1993, Pisa
1995, Lausanne 1996, Heraklion 1999
ACAT: “Advanced Computing and Analysis Techniques in Physics”
• Chicago 2000, Moscow 2002, Tsukuba 2003, Zeuthen 2005
•
•
•
•
•
•
WWW: Jan. 92, first presentation to the HEP community (TBL)
Symbolic/automatic calculations: building a community
Neural networks, genetic algorithms
Expert systems  agents
Parallel computing
Analysis tools
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
modern research relies on
the virtuous triangle
Computing
Theor
y
Experiment
Putting thinking, “intelligence”
into Systems
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
• Great plenary sessions (quantum computer, symbolic languages,
Neural nets, Berlin physics history, DESY research)
•
•
•
•
•
High quality presentations (High precision computation)
New topics
New contacts between participants
New collaborations initiated
Important Proceedings in preparation
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
• 2 years before the start of LHC,
– Where are the major question marks?
– Where should we concentrate our effort?
• ILC, … Are we ready for globalization of
research: Technology ? Sociology ?
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
1.
“BOINC”: Global computation
2.
The “WikiWay”: Better, simpler communication
3.
High precision floating point computation
4.
XP: eXtreme Programming
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
Berkeley Open Infrastructure for Network Computing
“Free” power
Average: 1.5 Gflops system
 250 Tflops
 ½ for 80% time 100 Tflops
Top 500 N.1- Blue Gene 70 TFlops
Climateprediction.net: study climate change
Einstein@home:
search for gravitational signals coming from pulsars
LHC@home:
improve the design of the CERN LHC particle accelerator
Predictor@home: investigate protein-related diseases
SETI@home:
Look for radio evidence of extraterrestrial life
Very important outreach
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
• A web site is “self ”-updated.
• The site becomes a cooperative
developing tool
• Stimulate the discussion
• Building a community
• Wikipedia…
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
•Security
•Project management
•Collaborative document writing:
•Notes, meeting report, presentation, draft paper
•Latex, Word interface (plugin)?
•Document indexing for in-text search
•Collaborative Code writing ??
•Aggregation (RSS) of web and application data in a structured way
•Database ..
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
•First time in this workshop.
•Quadruple/octuple precision is needed.
•Correct results.
•Faster algorithms.
•But software implementations are too slow.
New hardware/software development.
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
Simple Example
by J. Fujimoto (kek)
f = 333.75 b6 + a2(11a2b2- b6- 121b4 – 2) + 5.5b8 + a/2b
where a=77617.0, b=33096.0.
(C. Hu, S. Xu and X. Yang)
Quadruple precision result
f = 1.1726039400531786318588349045201801
Analytical result = - 54767/66192
f = - 0.82739605994682136814116509547981370
New Octuple precision library, H3Lib:
f = - 0.827396059946821368141165095479816
lost bits = 121
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
Actual application
By J. Fujimoto KEK
e  e   e  e   
t

s'
'
 : mass of 
 : mass of photon
DESY,Zeuthen
May 22-27, 2005
Quadruple precision is
required in some phase
space points due to the
Gram determinant  0
0
happens in the reduction
0 algorithm.
Denis Perret-Gallix
LAPP/IN2P3
  1.6,   10 20 ,
Double precision
ReJ[1]
ReJ[x]
ReJ[y]
ReJ[w]
ReJ[w**2]
ReJ[w*x]
ReJ[x*y]
…
ReJ[w**3]
…
= -1.49368718239238
= - 6.86111482424926E-0002
= - 6.86785270067264E-0002
= - 1.39799775179174
= - 1.36472026946296
= - 2.708863236843683E-0002
= - 3.048903558925384E-0002
s '  11.08527740873578,
t  1.46311226509359,
  2.007662965116379 10 11
 '  3.041076446800326 10  2
Blow up !!
= 93763.26727997246
Quadruple precision
ReJ[1]
ReJ[x]
ReJ[y]
ReJ[w]
ReJ[w**2]
ReJ[w*x]
ReJ[x*y]
…
ReJ[w**3]
= - 1.49368718238777512062307539882045
= - 6.861114708877389206553392789958382E-0002
= - 6.867852585600575199171661642779842E-0002
= - 1.39799775496536042464289674154150
= - 1.34746346742190735627641191119128
= - 3.334744118868393382280835719751654E-0002
= - 2.822377826411337874789947823777159E-0002
= - 1.60389378482142986480454883491878
…
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
J. Fujimoto
Minimization algorithms, …
Single
Double
Quad.
Gambolatin et al.
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
CPU bandwidth (Ghz)
Interconnection bandwidth (Ghz)
Memory size (Gbytes)
Floating point precision (4-32 bytes)
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
Instruction size (64 bits)
•High precision libraries quadruple/octuple (Hitatchi)
•Double-double, quad-double (Arprec)
•Multi-precision lib. (1000 dig and more)
•Interval arithmetic
•Exact arithmetic (XR, iRRAM)
•Linpack double/quad: 30 times slower
High precision Arithmetic Parallel Processor Yoke
HAPPY
Pulling heavy computations !!!
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
•Many short cycles: plan-code-test,
•frequent releases
•Close collaboration: users, programmers,
•daily meetings
•Standardization
•Language
•Structure
•Interface
•naming convention
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
More involvement of the community,
Discussion of innovative ideas, new approaches,
contradictory and hot issues, cross-experiments exchange
Not only at the next ACAT workshop,
but during the next 18-24 months
• Setting up of a Wiki ACAT web site
– A page for each working group discussion
– A repository for notes, presentations, papers related to the ACAT topics,
code, …
– Up-to-date reference list
– Collaborative writing of the History and trends for some topics
– ...
Security issues ??
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
•
Computing Technology and Environment for Physics Research
–
–
–
–
–
–
–
•
Data Analysis - Algorithms and Tools
–
–
–
–
–
–
•
Parallel and Vector Computing Technologies and Applications
Distributed Computing: Tools for the GRID, BOINC, …
Data Fabric and Data Management
Online Monitoring and Control
Software Engineering: eXtreme programming, …
Graphic User Interfaces, Common Libraries
Communication: web and Wiki’s tools
Neural Networks and Other Pattern Recognition Techniques
Evolutionary and Genetic Algorithms, Statistical learning
Advanced Data Analysis Environments
Statistical Methods
Detector and Accelerator Simulations
Reconstruction Algorithms
Topical sessions
Round-Tables
Simulations and Computations in Theoretical Physics and Phenomenology
–
–
–
–
–
Automatic Computation Systems: from Processes to Event Generators
Multi-loop Calculations and Higher Order Corrections
Multi-dimensional Integration and Event Generators
Intensive High Precision Numerical Computations: Algorithms and Systems
Computer Algebra Techniques and Applications: Intellectics
Quantum Computing
Lattice Computing
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
The
Netherlands
The
Netherlands
- Amsterdam - NIKHEF
Amsterdam
- Nihkef
April 2007,
2007,
ChairApril
of the organizing
committee:
Chair
of the
Jos Vermaseren
organizing
committee:
Candidates for backup and next-to-next:
Jos
Vermaseren
Switzerland( Zurich), US, France (Annecy), Norway,…
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3
Amsterdam April 2007
DESY,Zeuthen
May 22-27, 2005
Denis Perret-Gallix
LAPP/IN2P3