DESY Zeuthen May 22 27 2005 Denis PerretGallix

• • Institutions: DESY, Zeuthen, Bielefeld Univ. , Helmholtz-Gemeinschaft, EPS Local organizing committee:

AIHENP: “Artificial Intelligence in High Energy and Nuclear Physics” • Lyon 1990, la Londes-Les-Maures

Computin g modern research relies on the virtuous triangle Theor y Experiment Putting thinking,

• Great plenary sessions (quantum computer, symbolic languages, Neural nets, Berlin physics history,

• 2 years before the start of LHC, – Where are the major

1. “BOINC”: Global computation 2. The “Wiki. Way”: Better, simpler communication 3. High precision

Berkeley Open Infrastructure for Network Computing “Free” power Average: 1. 5 Gflops system Þ

• A web site is “self ”-updated. • The site becomes a cooperative

• Security • Project management • Collaborative document writing: • Notes, meeting report,

• First time in this workshop. • Quadruple/octuple precision is needed. • Correct

Simple Example by J. Fujimoto (kek) f = 333. 75 b 6 + a

Actual application By J. Fujimoto KEK t Quadruple precision is required in some phase

Double precision Re. J[1] Re. J[x] Re. J[y] Re. J[w**2] Re. J[w*x] Re. J[x*y]

Minimization algorithms, … Single Double Quad. Gambolatin et al. DESY, Zeuthen May 22 -27,

CPU bandwidth (Ghz) Interconnection bandwidth (Ghz) Memory size (Gbytes) Floating point precision (4 -32

• High precision libraries quadruple/octuple (Hitatchi) • Double-double, quad-double (Arprec) • Multi-precision lib.

• Many short cycles: plan-code-test, • frequent releases • Close collaboration: users, programmers,

More involvement of the community, Discussion of innovative ideas, new approaches, contradictory and hot

• Computing Technology and Environment for Physics Research – – – – •

The Netherlands - Amsterdam - NIKHEF Amsterdam - Nihkef April 2007, Chair. April of

Amsterdam April 2007 DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P

Slides: 23

Download presentation

DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

• • 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/IN 2 P 3

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/IN 2 P 3

Computin g modern research relies on the virtuous triangle Theor y Experiment Putting thinking, “intelligence” into Systems DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

• 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/IN 2 P 3

• 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/IN 2 P 3

1. “BOINC”: Global computation 2. The “Wiki. Way”: Better, simpler communication 3. High precision floating point computation 4. XP: e. Xtreme Programming DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

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/IN 2 P 3

• 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/IN 2 P 3

DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

• 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/IN 2 P 3

• 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/IN 2 P 3

Simple Example by J. Fujimoto (kek) f = 333. 75 b 6 + a 2(11 a 2 b 2 - b 6 - 121 b 4 – 2) + 5. 5 b 8 + a/2 b 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, H 3 Lib: f = - 0. 827396059946821368141165095479816 lost bits = 121 DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

Actual application By J. Fujimoto KEK t Quadruple precision is required in some phase space points due to the Gram determinant mass of photon DESY, Zeuthen May 22 -27, 2005 happens in the reduction algorithm. Denis Perret-Gallix LAPP/IN 2 P 3

Double precision Re. J[1] Re. J[x] Re. J[y] Re. J[w**2] Re. J[w*x] Re. J[x*y] … Re. J[w**3] … = -1. 49368718239238 = - 6. 86111482424926 E-0002 = - 6. 86785270067264 E-0002 = - 1. 39799775179174 = - 1. 36472026946296 = - 2. 708863236843683 E-0002 = - 3. 048903558925384 E-0002 Blow up !! = 93763. 26727997246 Quadruple precision Re. J[1] Re. J[x] Re. J[y] Re. J[w**2] Re. J[w*x] Re. J[x*y] … Re. J[w**3] = - 1. 49368718238777512062307539882045 = - 6. 861114708877389206553392789958382 E-0002 = - 6. 867852585600575199171661642779842 E-0002 = - 1. 39799775496536042464289674154150 = - 1. 3474634674219073562764119128 = - 3. 334744118868393382280835719751654 E-0002 = - 2. 822377826411337874789947823777159 E-0002 = - 1. 60389378482142986480454883491878 … DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3 J. Fujimoto

Minimization algorithms, … Single Double Quad. Gambolatin et al. DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3

CPU bandwidth (Ghz) Interconnection bandwidth (Ghz) Memory size (Gbytes) Floating point precision (4 -32 bytes) DESY, Zeuthen May 22 -27, 2005 Instruction size (64 bits) Denis Perret-Gallix LAPP/IN 2 P 3

• High precision libraries quadruple/octuple (Hitatchi) • Double-double, quad-double (Arprec) • Multi-precision lib. (1000 dig and more) • Interval arithmetic • Exact arithmetic (XR, i. RRAM) • 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/IN 2 P 3

• 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/IN 2 P 3

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/IN 2 P 3

• 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: e. Xtreme 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/IN 2 P 3

The Netherlands - Amsterdam - NIKHEF Amsterdam - Nihkef April 2007, Chair. April 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/IN 2 P 3

Amsterdam April 2007 DESY, Zeuthen May 22 -27, 2005 Denis Perret-Gallix LAPP/IN 2 P 3