Application of Maxwell Equations to Human Body Modelling

Size: px
Start display at page:

Download "Application of Maxwell Equations to Human Body Modelling"

Transcription

1 Application of Maxwell Equations to Human Body Modelling Fumie Costen Room E3, E0c at Sackville Street Building, November 6, 00 Fumie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 / 0

2 Outline Research background Research issues umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 / 0

3 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

4 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

5 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

6 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

7 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

8 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

9 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

10 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

11 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

12 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

13 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

14 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

15 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

16 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

17 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

18 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

19 Research background is a common form of heart disease Catheter ablation with RF radiation is a proven technique to remove arrhythmogenic foci umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

20 Research background For emergency, electric shock is applied for defibrillation Not high successful rate 3 The optimum size, location of electrolodes are unknown a Fumie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 4 / 0

21 Research background Brain disorders Amnesia, Alzehimer,Epilepsy, Parkinson s disease, Depression, Schizophrenia Medication and counselling for treatment a Fumie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 5 / 0

22 Research background Brain disorders Alternative treatment: deep brain stimulation Pacemaker implanted for treatment Fumie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 6 / 0

23 Research background Brain disorders Literally Invasive treatment Non-invasive treatment desired Fumie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 7 / 0 a

24 Research background Needs of numerical simulation How ElectroMagnetic (EM) wave propagates from certain location and shape of the electrolode on the torso for successful defibrillation from the skull to do the same effect as invasive treatment( brain stimulation ) a electroconvulsive-therapy-is-like Fumie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 8 / 0

25 Numerical simulation tool development Maxwell s equation finite difference time domain (FDTD) 3 most straightforward, robust, widely applicable 4 high accuracy possible with high computational cost umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

26 Numerical simulation tool development Maxwell s equation finite difference time domain (FDTD) 3 most straightforward, robust, widely applicable 4 high accuracy possible with high computational cost umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

27 Numerical simulation tool development Maxwell s equation finite difference time domain (FDTD) 3 most straightforward, robust, widely applicable 4 high accuracy possible with high computational cost umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

28 Numerical simulation tool development Maxwell s equation finite difference time domain (FDTD) 3 most straightforward, robust, widely applicable 4 high accuracy possible with high computational cost umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

29 Many tissues involved in the computation Digital Human Phantom from mm resolution MRI Fumie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 0 / 0

30 Each tissue is frequency-dependent Measurement and data fitting by US Air Force Fumie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 / 0

31 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

32 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

33 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

34 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

35 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

36 Wave propagation simulation inside human remeshing from mm resolution to 0.3 mm resolution loading of the fine geometrical detail to the Maxwell equation solver 3 allocation of the frequency dependent characteristics to each tissue 4 load-balanced parallel computation on EUgrid 5 data collection from EUgrid to our local cluster 6 visualisation of the data from the Maxwell equation solver Fumie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 / 0

37 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

38 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

39 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

40 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

41 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

42 How to solve partial differencial equations fast Original equation A φ = Bφ; B has spatial discretization t elements A φn+ φ n = B φn+ + φ n t (A tb)φ n+ = (A + tb)φ n Sparse matrix to handle 3 Approximate factorization *example* to handle tridiagonal matrix (I X L )(I Y L )(I Z L )φ n+ = (I + X R )(I + Y R )(I + Z R )φ n+ φ n how to factorise multiplicative solution additive solution umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 3 / 0

43 Reduction of memory requirement by sophisticated boundary condition The normal wave propagation Boundary placed close to the signal source umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 4 / 0

44 Reduction of memory requirement by Subgridding Signal leaking at the interface between the fine mesh and the coarse mesh Filtering out the leaking signal for stability Independent fine mesh FDTD computation only for heart Communication between the coarse and fine mesh at the interface umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 5 / 0

45 Efficient computation on EUgrid The production code requires more than 300 GB of memory for ECG simulation Each node can share upto 3 GB of memory with the latest technology 3 Distributed computing is essential 4 High scalability is aimed at with the sophisticated parallel computing umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 6 / 0

46 Efficient computation on EUgrid The production code requires more than 300 GB of memory for ECG simulation Each node can share upto 3 GB of memory with the latest technology 3 Distributed computing is essential 4 High scalability is aimed at with the sophisticated parallel computing umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 6 / 0

47 Efficient computation on EUgrid The production code requires more than 300 GB of memory for ECG simulation Each node can share upto 3 GB of memory with the latest technology 3 Distributed computing is essential 4 High scalability is aimed at with the sophisticated parallel computing umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 6 / 0

48 Efficient computation on EUgrid The production code requires more than 300 GB of memory for ECG simulation Each node can share upto 3 GB of memory with the latest technology 3 Distributed computing is essential 4 High scalability is aimed at with the sophisticated parallel computing umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 6 / 0

49 Efficient computation on EUgrid The concept of computation on the distributed memory architecture machine OpenMP for speedup machine 4 GB shared memory CPU MPI for increase of memory usage CPU3 4 GB shared memory CPU CPU4 Fumie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 7 / 0

50 Efficient computation on EUgrid The directory information tree User Element site information Computing Element Storage Element Network information between this and other sites status supported protocols file statistics umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 8 / 0

51 Computation on Graphics Processing Units How do we change the computational strategy suitable for GPU? How much computational gain do we get over the computation with MPI? umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

52 Computation on Graphics Processing Units How do we change the computational strategy suitable for GPU? How much computational gain do we get over the computation with MPI? umie Costen Room E3, E0c at Sackville Street Building, Application of Maxwell equations ( ) to Human Body November 6, 00 9 / 0

53 Finally... Welcome to Computation in Electromagnetics Hope you like our research domain/activity ( biomedical modelling, parallel computing in Beowulf clusters, parabolic partial differential equations, programming, shell-scripting ) Please talk to me for the detail with your future plan umie Costen Room E3, E0c at Sackville Street Building, Application fc@cs.man.ac.uk of Maxwell equations ( ) to Human Body November 6, 00 0 / 0

SPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics

SPARSE SOLVERS POISSON EQUATION. Margreet Nool. November 9, 2015 FOR THE. CWI, Multiscale Dynamics SPARSE SOLVERS FOR THE POISSON EQUATION Margreet Nool CWI, Multiscale Dynamics November 9, 2015 OUTLINE OF THIS TALK 1 FISHPACK, LAPACK, PARDISO 2 SYSTEM OVERVIEW OF CARTESIUS 3 POISSON EQUATION 4 SOLVERS

More information

Thermo-elastic Response of Cutaneous and Subcutaneous Tissues to Noninvasive Radiofrequency Heating

Thermo-elastic Response of Cutaneous and Subcutaneous Tissues to Noninvasive Radiofrequency Heating Thermo-elastic Response of Cutaneous and Subcutaneous Tissues to Noninvasive Radiofrequency Heating Joel N. Jiménez Lozano, Paulino Vacas-Jacques, Walfre Franco. Excerpt from the Proceedings of the 2012

More information

Welcome to MCS 572. content and organization expectations of the course. definition and classification

Welcome to MCS 572. content and organization expectations of the course. definition and classification Welcome to MCS 572 1 About the Course content and organization expectations of the course 2 Supercomputing definition and classification 3 Measuring Performance speedup and efficiency Amdahl s Law Gustafson

More information

Extension of the FDTD Huygens Subgridding to Frequency Dependent Media

Extension of the FDTD Huygens Subgridding to Frequency Dependent Media Noname manuscript No. (will be inserted by the editor) Extension of the FDTD Huygens Subgridding to Frequency Dependent Media Fumie Costen Jean-Pierre Bérenger eceived: date / Accepted: date Abstract A

More information

ORBIT Code Review and Future Directions. S. Cousineau, A. Shishlo, J. Holmes ECloud07

ORBIT Code Review and Future Directions. S. Cousineau, A. Shishlo, J. Holmes ECloud07 ORBIT Code Review and Future Directions S. Cousineau, A. Shishlo, J. Holmes ECloud07 ORBIT Code ORBIT (Objective Ring Beam Injection and Transport code) ORBIT is an object-oriented, open-source code started

More information

MPI at MPI. Jens Saak. Max Planck Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory

MPI at MPI. Jens Saak. Max Planck Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory MAX PLANCK INSTITUTE November 5, 2010 MPI at MPI Jens Saak Max Planck Institute for Dynamics of Complex Technical Systems Computational Methods in Systems and Control Theory FOR DYNAMICS OF COMPLEX TECHNICAL

More information

Scalable and Power-Efficient Data Mining Kernels

Scalable and Power-Efficient Data Mining Kernels Scalable and Power-Efficient Data Mining Kernels Alok Choudhary, John G. Searle Professor Dept. of Electrical Engineering and Computer Science and Professor, Kellogg School of Management Director of the

More information

Open-source finite element solver for domain decomposition problems

Open-source finite element solver for domain decomposition problems 1/29 Open-source finite element solver for domain decomposition problems C. Geuzaine 1, X. Antoine 2,3, D. Colignon 1, M. El Bouajaji 3,2 and B. Thierry 4 1 - University of Liège, Belgium 2 - University

More information

Concurrent Divide-and-Conquer Library

Concurrent Divide-and-Conquer Library with Petascale Electromagnetics Applications, Tech-X Corporation CScADS Workshop on Libraries and Algorithms for Petascale Applications, 07/30/2007, Snowbird, Utah Background Particle In Cell (PIC) in

More information

Dense Arithmetic over Finite Fields with CUMODP

Dense Arithmetic over Finite Fields with CUMODP Dense Arithmetic over Finite Fields with CUMODP Sardar Anisul Haque 1 Xin Li 2 Farnam Mansouri 1 Marc Moreno Maza 1 Wei Pan 3 Ning Xie 1 1 University of Western Ontario, Canada 2 Universidad Carlos III,

More information

Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano

Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS. Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano Parallel PIPS-SBB Multi-level parallelism for 2-stage SMIPS Lluís-Miquel Munguia, Geoffrey M. Oxberry, Deepak Rajan, Yuji Shinano ... Our contribution PIPS-PSBB*: Multi-level parallelism for Stochastic

More information

Acceleration of WRF on the GPU

Acceleration of WRF on the GPU Acceleration of WRF on the GPU Daniel Abdi, Sam Elliott, Iman Gohari Don Berchoff, Gene Pache, John Manobianco TempoQuest 1434 Spruce Street Boulder, CO 80302 720 726 9032 TempoQuest.com THE WORLD S FASTEST

More information

An H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages

An H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages PIERS ONLINE, VOL. 6, NO. 7, 2010 679 An H-LU Based Direct Finite Element Solver Accelerated by Nested Dissection for Large-scale Modeling of ICs and Packages Haixin Liu and Dan Jiao School of Electrical

More information

SCITECH Volume 4, Issue 1 RESEARCH ORGANISATION November 09, 2017

SCITECH Volume 4, Issue 1 RESEARCH ORGANISATION November 09, 2017 SCITECH Volume 4, Issue 1 RESEARCH ORGANISATION November 9, 17 Boson Journal of Modern Physics www.scitecresearch.com Numerical Study The Dielectric Properties And Specific Absorption Rate Of Nerve Human

More information

Static-scheduling and hybrid-programming in SuperLU DIST on multicore cluster systems

Static-scheduling and hybrid-programming in SuperLU DIST on multicore cluster systems Static-scheduling and hybrid-programming in SuperLU DIST on multicore cluster systems Ichitaro Yamazaki University of Tennessee, Knoxville Xiaoye Sherry Li Lawrence Berkeley National Laboratory MS49: Sparse

More information

Cactus Tools for Petascale Computing

Cactus Tools for Petascale Computing Cactus Tools for Petascale Computing Erik Schnetter Reno, November 2007 Gamma Ray Bursts ~10 7 km He Protoneutron Star Accretion Collapse to a Black Hole Jet Formation and Sustainment Fe-group nuclei Si

More information

Improvements for Implicit Linear Equation Solvers

Improvements for Implicit Linear Equation Solvers Improvements for Implicit Linear Equation Solvers Roger Grimes, Bob Lucas, Clement Weisbecker Livermore Software Technology Corporation Abstract Solving large sparse linear systems of equations is often

More information

Coupling atomistic and continuum modelling of magnetism

Coupling atomistic and continuum modelling of magnetism Coupling atomistic and continuum modelling of magnetism M. Poluektov 1,2 G. Kreiss 2 O. Eriksson 3 1 University of Warwick WMG International Institute for Nanocomposites Manufacturing 2 Uppsala University

More information

Arches Part 1: Introduction to the Uintah Computational Framework. Charles Reid Scientific Computing Summer Workshop July 14, 2010

Arches Part 1: Introduction to the Uintah Computational Framework. Charles Reid Scientific Computing Summer Workshop July 14, 2010 Arches Part 1: Introduction to the Uintah Computational Framework Charles Reid Scientific Computing Summer Workshop July 14, 2010 Arches Uintah Computational Framework Cluster Node Node Node Node Node

More information

An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors

An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Contemporary Mathematics Volume 218, 1998 B 0-8218-0988-1-03024-7 An Efficient FETI Implementation on Distributed Shared Memory Machines with Independent Numbers of Subdomains and Processors Michel Lesoinne

More information

EE-353 ELECTROMAGNETICS SIMULATION LABORATORY

EE-353 ELECTROMAGNETICS SIMULATION LABORATORY EE-353 ELECTROMAGNETICS SIMULATION LABORATORY DEPARTMENT OF ELECTRICAL ENGINEERING UNIVERSITY OF ENGINEERING & TECHNOLOGY, LAHORE Magnetic Field Around Current Carrying Conductors CONTACT INFORMATION Electrical

More information

TR A Comparison of the Performance of SaP::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems

TR A Comparison of the Performance of SaP::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems TR-0-07 A Comparison of the Performance of ::GPU and Intel s Math Kernel Library (MKL) for Solving Dense Banded Linear Systems Ang Li, Omkar Deshmukh, Radu Serban, Dan Negrut May, 0 Abstract ::GPU is a

More information

1) Electronic Circuits & Laboratory

1) Electronic Circuits & Laboratory ENSEA COURSES TAUGHT IN ENGLISH SPRING Semester 1) Electronic Circuits & Laboratory Lecture : 45 hours Laboratory : 45 hours US Credits : 6 Analysis of integrated amplifiers with bipolar junction transistors

More information

Medical Imaging Physics Spring Quarter Week 9-1

Medical Imaging Physics Spring Quarter Week 9-1 Medical Imaging Physics Spring Quarter Week 9-1 NMR and MRI Davor Balzar balzar@du.edu www.du.edu/~balzar Intro MRI Outline NMR & MRI Guest lecturer fmri Thursday, May 22 Visit to CUHSC It s not mandatory

More information

Sparse Principal Component Analysis via Alternating Maximization and Efficient Parallel Implementations

Sparse Principal Component Analysis via Alternating Maximization and Efficient Parallel Implementations Sparse Principal Component Analysis via Alternating Maximization and Efficient Parallel Implementations Martin Takáč The University of Edinburgh Joint work with Peter Richtárik (Edinburgh University) Selin

More information

MODERN engineering problems grow in size and at the

MODERN engineering problems grow in size and at the 4336 IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION, VOL. 60, NO. 9, SEPTEMBER 2012 Huygens Subgridding for 3-D Frequency-Dependent Finite-Difference Time-Domain Method Maksims Abaļenkovs, Member, IEEE,

More information

Acoustics Analysis of Speaker ANSYS, Inc. November 28, 2014

Acoustics Analysis of Speaker ANSYS, Inc. November 28, 2014 Acoustics Analysis of Speaker 1 Introduction ANSYS 14.0 offers many enhancements in the area of acoustics. In this presentation, an example speaker analysis will be shown to highlight some of the acoustics

More information

User Documentation and Examples (II) in GEANT p01

User Documentation and Examples (II) in GEANT p01 User Documentation and Examples (II) in GEANT 4.9.3-p01 Michael H. Kelsey SLAC National Accelerator Laboratory GEANT4 Tutorial, BUAF Puebla, Mexico 14 Jun 2010 Advanced User Documentation Toolkit developers

More information

Parallelization of the Molecular Orbital Program MOS-F

Parallelization of the Molecular Orbital Program MOS-F Parallelization of the Molecular Orbital Program MOS-F Akira Asato, Satoshi Onodera, Yoshie Inada, Elena Akhmatskaya, Ross Nobes, Azuma Matsuura, Atsuya Takahashi November 2003 Fujitsu Laboratories of

More information

CAEFEM v9.5 Information

CAEFEM v9.5 Information CAEFEM v9.5 Information Concurrent Analysis Corporation, 50 Via Ricardo, Thousand Oaks, CA 91320 USA Tel. (805) 375 1060, Fax (805) 375 1061 email: info@caefem.com or support@caefem.com Web: http://www.caefem.com

More information

Tight-Focusing of Short Intense Laser Beams in Particle-in-Cell Simulations of Laser-Plasma Interaction

Tight-Focusing of Short Intense Laser Beams in Particle-in-Cell Simulations of Laser-Plasma Interaction 28/03/2017, CTU in Prague Tight-Focusing of Short Intense Laser Beams in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej

More information

High Performance Computing for Neutron Tomography Reconstruction

High Performance Computing for Neutron Tomography Reconstruction High Performance Computing for Neutron Tomography Reconstruction A Parallel Approach to Filtered Backprojection (FBP) Zongpu Li 1 Cain Gantt 2 1 Department of Physics and Materials Science City University

More information

A robust multilevel approximate inverse preconditioner for symmetric positive definite matrices

A robust multilevel approximate inverse preconditioner for symmetric positive definite matrices DICEA DEPARTMENT OF CIVIL, ENVIRONMENTAL AND ARCHITECTURAL ENGINEERING PhD SCHOOL CIVIL AND ENVIRONMENTAL ENGINEERING SCIENCES XXX CYCLE A robust multilevel approximate inverse preconditioner for symmetric

More information

Improving weather prediction via advancing model initialization

Improving weather prediction via advancing model initialization Improving weather prediction via advancing model initialization Brian Etherton, with Christopher W. Harrop, Lidia Trailovic, and Mark W. Govett NOAA/ESRL/GSD 15 November 2016 The HPC group at NOAA/ESRL/GSD

More information

FEAST eigenvalue algorithm and solver: review and perspectives

FEAST eigenvalue algorithm and solver: review and perspectives FEAST eigenvalue algorithm and solver: review and perspectives Eric Polizzi Department of Electrical and Computer Engineering University of Masachusetts, Amherst, USA Sparse Days, CERFACS, June 25, 2012

More information

Phased Array Feed Scattering Analysis. Stuart Hay, Raj Mittra, Neng-Tien Huang and Wenhua Wu

Phased Array Feed Scattering Analysis. Stuart Hay, Raj Mittra, Neng-Tien Huang and Wenhua Wu Phased Array Feed Scattering Analysis Stuart Hay, Raj Mittra, Neng-Tien Huang and Wenhua Wu 3 May 200 Outline of talk Motivation Phased array feed (PAF) scattering analysis Two approaches Characteristic

More information

Real time network modulation for intractable epilepsy

Real time network modulation for intractable epilepsy Real time network modulation for intractable epilepsy Behnaam Aazhang! Electrical and Computer Engineering Rice University! and! Center for Wireless Communications University of Oulu Finland Real time

More information

GPU-Enabled Spatiotemporal Model of Stochastic Cardiac Calcium Dynamics and Arrhythmias

GPU-Enabled Spatiotemporal Model of Stochastic Cardiac Calcium Dynamics and Arrhythmias GPU-Enabled Spatiotemporal Model of Stochastic Cardiac Calcium Dynamics and Arrhythmias M. Saleet Jafri Hoang Trong Minh Tuan Department of Bioinformatics and Computational Biology George Mason University

More information

上海超级计算中心 Shanghai Supercomputer Center. Lei Xu Shanghai Supercomputer Center San Jose

上海超级计算中心 Shanghai Supercomputer Center. Lei Xu Shanghai Supercomputer Center San Jose 上海超级计算中心 Shanghai Supercomputer Center Lei Xu Shanghai Supercomputer Center 03/26/2014 @GTC, San Jose Overview Introduction Fundamentals of the FDTD method Implementation of 3D UPML-FDTD algorithm on GPU

More information

INITIAL INTEGRATION AND EVALUATION

INITIAL INTEGRATION AND EVALUATION INITIAL INTEGRATION AND EVALUATION OF SLATE PARALLEL BLAS IN LATTE Marc Cawkwell, Danny Perez, Arthur Voter Asim YarKhan, Gerald Ragghianti, Jack Dongarra, Introduction The aim of the joint milestone STMS10-52

More information

Direct Self-Consistent Field Computations on GPU Clusters

Direct Self-Consistent Field Computations on GPU Clusters Direct Self-Consistent Field Computations on GPU Clusters Guochun Shi, Volodymyr Kindratenko National Center for Supercomputing Applications University of Illinois at UrbanaChampaign Ivan Ufimtsev, Todd

More information

A Finite-Element based Navier-Stokes Solver for LES

A Finite-Element based Navier-Stokes Solver for LES A Finite-Element based Navier-Stokes Solver for LES W. Wienken a, J. Stiller b and U. Fladrich c. a Technische Universität Dresden, Institute of Fluid Mechanics (ISM) b Technische Universität Dresden,

More information

One Optimized I/O Configuration per HPC Application

One Optimized I/O Configuration per HPC Application One Optimized I/O Configuration per HPC Application Leveraging I/O Configurability of Amazon EC2 Cloud Mingliang Liu, Jidong Zhai, Yan Zhai Tsinghua University Xiaosong Ma North Carolina State University

More information

Multiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU

Multiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU Multiphase Flow Simulations in Inclined Tubes with Lattice Boltzmann Method on GPU Khramtsov D.P., Nekrasov D.A., Pokusaev B.G. Department of Thermodynamics, Thermal Engineering and Energy Saving Technologies,

More information

The Quantum Supremacy Experiment

The Quantum Supremacy Experiment The Quantum Supremacy Experiment John Martinis, Google & UCSB New tests of QM: Does QM work for 10 15 Hilbert space? Does digitized error model also work? Demonstrate exponential computing power: Check

More information

Recent Progress of Parallel SAMCEF with MUMPS MUMPS User Group Meeting 2013

Recent Progress of Parallel SAMCEF with MUMPS MUMPS User Group Meeting 2013 Recent Progress of Parallel SAMCEF with User Group Meeting 213 Jean-Pierre Delsemme Product Development Manager Summary SAMCEF, a brief history Co-simulation, a good candidate for parallel processing MAAXIMUS,

More information

Simulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling. Henrik Nordborg HSR University of Applied Sciences Rapperswil

Simulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling. Henrik Nordborg HSR University of Applied Sciences Rapperswil Simulations of Electrical Arcs: Algorithms, Physical Scales, and Coupling Henrik Nordborg HSR University of Applied Sciences Rapperswil What is an electrical arc? 2 Technical applications of arcs and industrial

More information

Projectile Motion Slide 1/16. Projectile Motion. Fall Semester. Parallel Computing

Projectile Motion Slide 1/16. Projectile Motion. Fall Semester. Parallel Computing Projectile Motion Slide 1/16 Projectile Motion Fall Semester Projectile Motion Slide 2/16 Topic Outline Historical Perspective ABC and ENIAC Ballistics tables Projectile Motion Air resistance Euler s method

More information

Stochastic Modelling of Electron Transport on different HPC architectures

Stochastic Modelling of Electron Transport on different HPC architectures Stochastic Modelling of Electron Transport on different HPC architectures www.hp-see.eu E. Atanassov, T. Gurov, A. Karaivan ova Institute of Information and Communication Technologies Bulgarian Academy

More information

Improvements to wave codes

Improvements to wave codes Improvements to wave codes John C. Wright P. T. Bonoli Haru Kohno and Jungpyo Lee 1 1 Plasma Science and Fusion Center Massachusetts Institute of Technology Sep 2010 - PPPL 1 RF SciDAC PPPL 2010 Contents

More information

Model Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University

Model Order Reduction via Matlab Parallel Computing Toolbox. Istanbul Technical University Model Order Reduction via Matlab Parallel Computing Toolbox E. Fatih Yetkin & Hasan Dağ Istanbul Technical University Computational Science & Engineering Department September 21, 2009 E. Fatih Yetkin (Istanbul

More information

Radiofrequency Dosimetry in Subjects Implanted with Metallic Structures Undergoing MRI: a Numerical Study

Radiofrequency Dosimetry in Subjects Implanted with Metallic Structures Undergoing MRI: a Numerical Study American Journal of Biomedical Sciences ISSN: 1937-9080 nwpii.com/ajbms Radiofrequency Dosimetry in Subjects Implanted with Metallic Structures Undergoing MRI: a Numerical Study E Mattei 1, M Triventi

More information

Recommended Courses by ECE Topic Area Graduate Students

Recommended Courses by ECE Topic Area Graduate Students Recommended s by ECE Topic Area Graduate Students The course recommendations below do not represent full plans of study. The courses listed under each heading represent appropriate courses to take if you

More information

Impression Store: Compressive Sensing-based Storage for. Big Data Analytics

Impression Store: Compressive Sensing-based Storage for. Big Data Analytics Impression Store: Compressive Sensing-based Storage for Big Data Analytics Jiaxing Zhang, Ying Yan, Liang Jeff Chen, Minjie Wang, Thomas Moscibroda & Zheng Zhang Microsoft Research The Curse of O(N) in

More information

Computational Electromagnetics Definitions, applications and research

Computational Electromagnetics Definitions, applications and research Computational Electromagnetics Definitions, applications and research Luis E. Tobón Pontificia Universidad Javeriana Seminario de investigación Departamento de Electrónica y Ciencias de la Computación

More information

Assignment on iterative solution methods and preconditioning

Assignment on iterative solution methods and preconditioning Division of Scientific Computing, Department of Information Technology, Uppsala University Numerical Linear Algebra October-November, 2018 Assignment on iterative solution methods and preconditioning 1.

More information

High-performance processing and development with Madagascar. July 24, 2010 Madagascar development team

High-performance processing and development with Madagascar. July 24, 2010 Madagascar development team High-performance processing and development with Madagascar July 24, 2010 Madagascar development team Outline 1 HPC terminology and frameworks 2 Utilizing data parallelism 3 HPC development with Madagascar

More information

Jacobi-Based Eigenvalue Solver on GPU. Lung-Sheng Chien, NVIDIA

Jacobi-Based Eigenvalue Solver on GPU. Lung-Sheng Chien, NVIDIA Jacobi-Based Eigenvalue Solver on GPU Lung-Sheng Chien, NVIDIA lchien@nvidia.com Outline Symmetric eigenvalue solver Experiment Applications Conclusions Symmetric eigenvalue solver The standard form is

More information

Some notes on efficient computing and setting up high performance computing environments

Some notes on efficient computing and setting up high performance computing environments Some notes on efficient computing and setting up high performance computing environments Andrew O. Finley Department of Forestry, Michigan State University, Lansing, Michigan. April 17, 2017 1 Efficient

More information

Crossing the Chasm. On the Paths to Exascale: Presented by Mike Rezny, Monash University, Australia

Crossing the Chasm. On the Paths to Exascale: Presented by Mike Rezny, Monash University, Australia On the Paths to Exascale: Crossing the Chasm Presented by Mike Rezny, Monash University, Australia michael.rezny@monash.edu Crossing the Chasm meeting Reading, 24 th October 2016 Version 0.1 In collaboration

More information

TIME DEPENDENCE OF SHELL MODEL CALCULATIONS 1. INTRODUCTION

TIME DEPENDENCE OF SHELL MODEL CALCULATIONS 1. INTRODUCTION Mathematical and Computational Applications, Vol. 11, No. 1, pp. 41-49, 2006. Association for Scientific Research TIME DEPENDENCE OF SHELL MODEL CALCULATIONS Süleyman Demirel University, Isparta, Turkey,

More information

Behavioral Simulations in MapReduce

Behavioral Simulations in MapReduce Behavioral Simulations in MapReduce Guozhang Wang, Marcos Vaz Salles, Benjamin Sowell, Xun Wang, Tuan Cao, Alan Demers, Johannes Gehrke, Walker White Cornell University 1 What are Behavioral Simulations?

More information

CRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel?

CRYSTAL in parallel: replicated and distributed (MPP) data. Why parallel? CRYSTAL in parallel: replicated and distributed (MPP) data Roberto Orlando Dipartimento di Chimica Università di Torino Via Pietro Giuria 5, 10125 Torino (Italy) roberto.orlando@unito.it 1 Why parallel?

More information

Numerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes. Kai Diethelm. International Symposium on Fractional PDEs June 3 5, 2013

Numerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes. Kai Diethelm. International Symposium on Fractional PDEs June 3 5, 2013 Numerical Fractional Calculus Using Methods Based on Non-Uniform Step Sizes Gesellschaft für numerische Simulation mbh Braunschweig AG Numerik Institut Computational Mathematics Technische Universität

More information

Using R for Iterative and Incremental Processing

Using R for Iterative and Incremental Processing Using R for Iterative and Incremental Processing Shivaram Venkataraman, Indrajit Roy, Alvin AuYoung, Robert Schreiber UC Berkeley and HP Labs UC BERKELEY Big Data, Complex Algorithms PageRank (Dominant

More information

Utilisation de la compression low-rank pour réduire la complexité du solveur PaStiX

Utilisation de la compression low-rank pour réduire la complexité du solveur PaStiX Utilisation de la compression low-rank pour réduire la complexité du solveur PaStiX 26 Septembre 2018 - JCAD 2018 - Lyon Grégoire Pichon, Mathieu Faverge, Pierre Ramet, Jean Roman Outline 1. Context 2.

More information

Contents. Preface... xi. Introduction...

Contents. Preface... xi. Introduction... Contents Preface... xi Introduction... xv Chapter 1. Computer Architectures... 1 1.1. Different types of parallelism... 1 1.1.1. Overlap, concurrency and parallelism... 1 1.1.2. Temporal and spatial parallelism

More information

Numerically Solving Partial Differential Equations

Numerically Solving Partial Differential Equations Numerically Solving Partial Differential Equations Michael Lavell Department of Applied Mathematics and Statistics Abstract The physics describing the fundamental principles of fluid dynamics can be written

More information

Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation

Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Introduction Solving Time-Harmonic Scattering Problems by the Ultra Weak Variational Formulation Plane waves as basis functions Peter Monk 1 Tomi Huttunen 2 1 Department of Mathematical Sciences University

More information

Jacobi-Davidson Eigensolver in Cusolver Library. Lung-Sheng Chien, NVIDIA

Jacobi-Davidson Eigensolver in Cusolver Library. Lung-Sheng Chien, NVIDIA Jacobi-Davidson Eigensolver in Cusolver Library Lung-Sheng Chien, NVIDIA lchien@nvidia.com Outline CuSolver library - cusolverdn: dense LAPACK - cusolversp: sparse LAPACK - cusolverrf: refactorization

More information

31545 Medical Imaging systems

31545 Medical Imaging systems Simulation of ultrasound systems and non-linear imaging 545 Medical Imaging systems Lecture 9: Simulation of ultrasound systems and non-linear imaging Jørgen Arendt Jensen Department of Electrical Engineering

More information

Performance Evaluation of Scientific Applications on POWER8

Performance Evaluation of Scientific Applications on POWER8 Performance Evaluation of Scientific Applications on POWER8 2014 Nov 16 Andrew V. Adinetz 1, Paul F. Baumeister 1, Hans Böttiger 3, Thorsten Hater 1, Thilo Maurer 3, Dirk Pleiter 1, Wolfram Schenck 4,

More information

ERROR CONVERGENCE ANALYSIS FOR LOCAL HYPERTHERMIA APPLICATIONS

ERROR CONVERGENCE ANALYSIS FOR LOCAL HYPERTHERMIA APPLICATIONS Journal of Engineering Science and Technology Vol. 11, No. 1 (2016) 060-067 School of Engineering, Taylor s University ERROR CONVERGENCE ANALYSIS FOR LOCAL HYPERTHERMIA APPLICATIONS NEERU MALHOTRA 1, *,

More information

Toward Binary Black Hole Simulations in Numerical Relativity

Toward Binary Black Hole Simulations in Numerical Relativity Toward Binary Black Hole Simulations in Numerical Relativity Frans Pretorius California Institute of Technology BIRS Workshop on Numerical Relativity Banff, April 19 2005 Outline generalized harmonic coordinates

More information

arxiv:physics/ v1 19 Feb 2004

arxiv:physics/ v1 19 Feb 2004 Finite Difference Time Domain (FDTD) Simulations of Electromagnetic Wave Propagation Using a Spreadsheet David W. Ward and Keith A. Nelson Department of Chemistry Massachusetts Institute of Technology,

More information

Computation of Large Sparse Aggregated Areas for Analytic Database Queries

Computation of Large Sparse Aggregated Areas for Analytic Database Queries Computation of Large Sparse Aggregated Areas for Analytic Database Queries Steffen Wittmer Tobias Lauer Jedox AG Collaborators: Zurab Khadikov Alexander Haberstroh Peter Strohm Business Intelligence and

More information

Practical Combustion Kinetics with CUDA

Practical Combustion Kinetics with CUDA Funded by: U.S. Department of Energy Vehicle Technologies Program Program Manager: Gurpreet Singh & Leo Breton Practical Combustion Kinetics with CUDA GPU Technology Conference March 20, 2015 Russell Whitesides

More information

Parallel programming using MPI. Analysis and optimization. Bhupender Thakur, Jim Lupo, Le Yan, Alex Pacheco

Parallel programming using MPI. Analysis and optimization. Bhupender Thakur, Jim Lupo, Le Yan, Alex Pacheco Parallel programming using MPI Analysis and optimization Bhupender Thakur, Jim Lupo, Le Yan, Alex Pacheco Outline l Parallel programming: Basic definitions l Choosing right algorithms: Optimal serial and

More information

S8241 VERSIONING GPU- ACCLERATED WRF TO Jeff Adie, 26 March, 2018 (Presented by Stan Posey, NVIDIA)

S8241 VERSIONING GPU- ACCLERATED WRF TO Jeff Adie, 26 March, 2018 (Presented by Stan Posey, NVIDIA) S8241 VERSIONING GPU- ACCLERATED WRF TO 3.7.1 Jeff Adie, 26 March, 2018 (Presented by Stan Posey, NVIDIA) 1 ACKNOWLEDGEMENT The work presented here today would not have been possible without the efforts

More information

Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction

Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction 16/05/2017, CTU in Prague Tight-Focusing of Short Intense Laser Pulses in Particle-in-Cell Simulations of Laser-Plasma Interaction Bc. Petr Valenta (petr.valenta@eli-beams.eu) Supervisors: doc. Ing. Ondrej

More information

EMF PENETRATION IN BIOLOGICAL TISSUE WHEN EXPOSED IN THE NEAR FIELD OF A MOBILE PHONE ANTENNA

EMF PENETRATION IN BIOLOGICAL TISSUE WHEN EXPOSED IN THE NEAR FIELD OF A MOBILE PHONE ANTENNA EMF PENETRATION IN BIOLOGICAL TISSUE WHEN EXPOSED IN THE NEAR FIELD OF A MOBILE PHONE ANTENNA Mihaela Morega, Alina Machedon POLITEHNICA University of Bucharest, mihaela@iem.pub.ro Abstract. The paper

More information

Technology Computer Aided Design (TCAD) Laboratory. Lecture 2, A simulation primer

Technology Computer Aided Design (TCAD) Laboratory. Lecture 2, A simulation primer Technology Computer Aided Design (TCAD) Laboratory Lecture 2, A simulation primer [Source: Synopsys] Giovanni Betti Beneventi E-mail: gbbeneventi@arces.unibo.it ; giobettibeneventi@gmail.com Office: Engineering

More information

Open-Source Parallel FE Software : FrontISTR -- Performance Considerations about B/F (Byte per Flop) of SpMV on K-Supercomputer and GPU-Clusters --

Open-Source Parallel FE Software : FrontISTR -- Performance Considerations about B/F (Byte per Flop) of SpMV on K-Supercomputer and GPU-Clusters -- Parallel Processing for Energy Efficiency October 3, 2013 NTNU, Trondheim, Norway Open-Source Parallel FE Software : FrontISTR -- Performance Considerations about B/F (Byte per Flop) of SpMV on K-Supercomputer

More information

N-body Simulations. On GPU Clusters

N-body Simulations. On GPU Clusters N-body Simulations On GPU Clusters Laxmikant Kale Filippo Gioachin Pritish Jetley Thomas Quinn Celso Mendes Graeme Lufkin Amit Sharma Joachim Stadel Lukasz Wesolowski James Wadsley Edgar Solomonik Fabio

More information

@SoyGema GEMA PARREÑO PIQUERAS

@SoyGema GEMA PARREÑO PIQUERAS @SoyGema GEMA PARREÑO PIQUERAS WHAT IS AN ARTIFICIAL NEURON? WHAT IS AN ARTIFICIAL NEURON? Image Recognition Classification using Softmax Regressions and Convolutional Neural Networks Languaje Understanding

More information

A Fast, Parallel Potential Flow Solver

A Fast, Parallel Potential Flow Solver Advisor: Jaime Peraire December 16, 2012 Outline 1 Introduction to Potential FLow 2 The Boundary Element Method 3 The Fast Multipole Method 4 Discretization 5 Implementation 6 Results 7 Conclusions Why

More information

Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems

Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Scalable Domain Decomposition Preconditioners For Heterogeneous Elliptic Problems Pierre Jolivet, F. Hecht, F. Nataf, C. Prud homme Laboratoire Jacques-Louis Lions Laboratoire Jean Kuntzmann INRIA Rocquencourt

More information

MUMPS. The MUMPS library: work done during the SOLSTICE project. MUMPS team, Lyon-Grenoble, Toulouse, Bordeaux

MUMPS. The MUMPS library: work done during the SOLSTICE project. MUMPS team, Lyon-Grenoble, Toulouse, Bordeaux The MUMPS library: work done during the SOLSTICE project MUMPS team, Lyon-Grenoble, Toulouse, Bordeaux Sparse Days and ANR SOLSTICE Final Workshop June MUMPS MUMPS Team since beg. of SOLSTICE (2007) Permanent

More information

RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS

RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS RECENT DEVELOPMENTS IN COMPUTATIONAL REACTOR ANALYSIS Dean Wang April 30, 2015 24.505 Nuclear Reactor Physics Outline 2 Introduction and Background Coupled T-H/Neutronics Safety Analysis Numerical schemes

More information

Coupling Physics. Tomasz Stelmach Senior Application Engineer

Coupling Physics. Tomasz Stelmach Senior Application Engineer Coupling Physics Tomasz Stelmach Senior Application Engineer Agenda Brief look @ Multiphysics solution What is new in R18 Fluent Maxwell coupling wireless power transfer Brief look @ ANSYS Multiphysics

More information

CONTEMPORARY ANALYTICAL ECOSYSTEM PATRICK HALL, SAS INSTITUTE

CONTEMPORARY ANALYTICAL ECOSYSTEM PATRICK HALL, SAS INSTITUTE CONTEMPORARY ANALYTICAL ECOSYSTEM PATRICK HALL, SAS INSTITUTE Copyright 2013, SAS Institute Inc. All rights reserved. Agenda (Optional) History Lesson 2015 Buzzwords Machine Learning for X Citizen Data

More information

Optimising PICCANTE an Open Source Particle-in-Cell Code for Advanced Simulations on Tier-0 Systems

Optimising PICCANTE an Open Source Particle-in-Cell Code for Advanced Simulations on Tier-0 Systems Available online at www.prace-ri.eu Partnership for Advanced Computing in Europe Optimising PICCANTE an Open Source Particle-in-Cell Code for Advanced Simulations on Tier-0 Systems A. Sgattoni a, L. Fedeli

More information

Modeling 3-D Calcium Waves from Stochastic Calcium sparks in a Sarcomere Using COMSOL

Modeling 3-D Calcium Waves from Stochastic Calcium sparks in a Sarcomere Using COMSOL Modeling 3-D Calcium Waves from Stochastic Calcium sparks in a Sarcomere Using COMSOL Zana Coulibaly 1, Leighton T. Izu 2 and Bradford E. Peercy 1 1 University of Maryland, Baltimore County 2 University

More information

BSEE REQUIREMENTS

BSEE REQUIREMENTS College of Engineering The Klipsch School of Electrical and Computer Engineering BSEE REQUIREMENTS 2014-2015 This document presents a summary of the requirements for earning a Bachelor of Science degree

More information

Impedance Imaging of the Thorax: Why it s difficult, and what we are doing about it?

Impedance Imaging of the Thorax: Why it s difficult, and what we are doing about it? Impedance Imaging of the Thorax: Why it s difficult, and what we are doing about it? Biomedical Engineering Research Centre (BIRC) Western Univerity, London, ON, 6 May 2015 Andy Adler Professor & Canada

More information

Heterogeneous programming for hybrid CPU-GPU systems: Lessons learned from computational chemistry

Heterogeneous programming for hybrid CPU-GPU systems: Lessons learned from computational chemistry Heterogeneous programming for hybrid CPU-GPU systems: Lessons learned from computational chemistry and Eugene DePrince Argonne National Laboratory (LCF and CNM) (Eugene moved to Georgia Tech last week)

More information

Research. Ji Chen Department of Electrical and Computer Engineering University of Houston Houston, TX 77204

Research. Ji Chen Department of Electrical and Computer Engineering University of Houston Houston, TX 77204 EMC/EMI Issues in Biomedical Research Ji Chen Department of Electrical and Computer Engineering University of Houston Houston, TX 77204 Email: jchen18@uh.eduedu UH: close to downtown of Houston 37,000

More information

Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations

Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations Robust Preconditioned Conjugate Gradient for the GPU and Parallel Implementations Rohit Gupta, Martin van Gijzen, Kees Vuik GPU Technology Conference 2012, San Jose CA. GPU Technology Conference 2012,

More information

Biomedical Instrumentation System

Biomedical Instrumentation System BME/EECS 458 - Biomedical Instrumentation & Design Matt O Donnell I.0 Introduction What is a biomedical instrument? To many it s an EKG machine, to others it s a chemical biosensor, and to some it s a

More information

ZFS - Flow Solver AIA NVIDIA Workshop

ZFS - Flow Solver AIA NVIDIA Workshop ZFS - Flow Solver AIA NVIDIA Workshop Andreas Lintermann, Matthias Meinke, Wolfgang Schröder Institute of Aerodynamics and Chair of Fluid Mechanics RWTH Aachen University Outline Softwareproject ZFS (applications

More information