Utilizing machine topology in numerical algorithms Luke Olson Department of Computer Science Computational Science and Engineering University of Illinois at Urbana-Champaign
Overview • I use Blue Waters to solve large sparse linear systems and to study the performance [ ][][] A b x = • Why? Time step other Linearization Assembly solve Solve Adapt Image: David Keyes??
Application: Plasma-coupled Combustion • The Center for Exascale Simulation of Plasma-Coupled Combustion XPACC @ Illinois Oxidizer Fuel O - Ar; Air H ; CH 2 2 4
<latexit sha1_base64="gyIeNogWzLnXIxRGRrbsZ lcO3k=">A CDnicbVA9SwNBEJ3zM8avU0stFkPAKtzZaCMEbSwTMB+QC2Fus0mW7O0du3tC PkFNv4VGxF bK3t7Pwpbj4ETXw 8Hhvhpl5YSK4Np736Swtr6yurWc2sptb2zu7 t5+VcepoqxCYxGreoiaCS5ZxXAjWD1RDKNQsFrYvxr7tVumNI/ljRk rBlhV/IOp2is1HLzgcRQIAloOzYk0LwbYUuRHzXpcXJBui035xW8Ccgi8WckVzx6Kn8BQKnlfgTtmKYRk4YK1Lrhe4lpDlEZTgUbZYNUswRpH7usYanEiOnmcPLOiOSt0iadWNmShkzU3xNDjLQeRKHtjND09Lw3Fv/zGqnpnDeHXCapYZJOF3VSQUxMxtmQNleMGjGwBKni9lZCe6iQGptg1obgz7+8SKqnBd8r+GWbxiVMkYFDOIYT8OEMinANJagAhTt4gGd4ce6dR+fVeZu2LjmzmQP4A+f9G4wNnVc=</latexit> <latexit sha1_base64="6dsFXVh1S8TqVSgnIwp21R0hVj4=">A CDnicbVC7SgNBFJ2Nrxhfq5aKDIaAVdi10UYI2lgmYB6QDcvdySQZMju7zMwKYUlpZeOv2EhQxNbazm/wJ5w8BE08cOFwzr3ce08Qc6a043xamaXl dW17HpuY3Nre8fe3aupKJGEVknEI9kIQFHOBK1qpjltxJ CGHBaD/pXY79+S6VikbjRg5i2QugK1mE tJF8u+AJCDhgj7Qj T3FuiH4Ev+ocY/hC9z17bxTdCbAi8SdkXzpcFT5ujsalX37w2tHJAmp0ISDUk3XiXUrBakZ4XSY8xJFYyB96NKmoQJCqlrp5J0hLhiljTuRNCU0nqi/J1I lRqEgekMQf UvDcW/ Oaie6ct1Im4kRTQa LOgnHOsLjbHCbSUo0HxgCRDJzKyY9kEC0STBnQnDnX14ktdOi6xTdiknjEk2R QfoGJ0gF52hErpGZVRFBN2jR/SMXqwH68l6td6mrRlrNrOP/sB6/wZqUp69</latexit> <latexit sha1_base64="6dsFXVh1S8TqVSgnIwp21R0hVj4=">A CDnicbVC7SgNBFJ2Nrxhfq5aKDIaAVdi10UYI2lgmYB6QDcvdySQZMju7zMwKYUlpZeOv2EhQxNbazm/wJ5w8BE08cOFwzr3ce08Qc6a043xamaXl dW17HpuY3Nre8fe3aupKJGEVknEI9kIQFHOBK1qpjltxJ CGHBaD/pXY79+S6VikbjRg5i2QugK1mE tJF8u+AJCDhgj7Qj T3FuiH4Ev+ocY/hC9z17bxTdCbAi8SdkXzpcFT5ujsalX37w2tHJAmp0ISDUk3XiXUrBakZ4XSY8xJFYyB96NKmoQJCqlrp5J0hLhiljTuRNCU0nqi/J1I lRqEgekMQf UvDcW/ Oaie6ct1Im4kRTQa LOgnHOsLjbHCbSUo0HxgCRDJzKyY9kEC0STBnQnDnX14ktdOi6xTdiknjEk2R QfoGJ0gF52hErpGZVRFBN2jR/SMXqwH68l6td6mrRlrNrOP/sB6/wZqUp69</latexit> <latexit sha1_base64="cA6f1IiS0YtSmUKcNY5x0R5U1ug=">A CDnicbVBNS8NAEN34WetX1KOXxVLwVBIvehGKXjxWsB/QlD ZbNKlm03Y3Qgl9Bd48a948aCIV8/e/Ddu2wja+mDg8d4M /OCjDOlHefLWl dW9/YrGxVt3d29/btg8O SnNJaJukPJW9ABTlTNC2ZprTXiYpJAGn3WB0PfW791Qqlo 7Pc7oI FYsIgR0Eby7bonIOCAPRKmGnuKxQn4Ev+o2ZDhSxz7ds1pODPgZeKWpIZKtHz70wtTkidUaMJBqb7rZHpQgNSMcDqpermiGZARxLRvqICEqkExe2eC60YJcZRKU0Ljmfp7o BEqXESmM4E9FAtelPxP6+f6+hiUDCR5ZoKMl8U5RzrFE+zwSGTlGg+NgSIZOZWTIYg WiTYNWE4C6+vEw6Zw3Xabi3Tq15VcZRQcfoBJ0iF52jJrpBLdRGBD2gJ/SCXq1H69l6s97nrStWOXOE/sD6+AZ8E5sS</latexit> Application: Plasma-coupled Combustion • The Center for Exascale Simulation of Plasma-Coupled Combustion XPACC @ Illinois • Electric conductivity influences the electric field and current density over time r · σ r r φ = g
Application: Plasma-coupled Combustion • The Center for Exascale Simulation of Plasma-Coupled Combustion @ Illinois XPACC • Electric field a key element in the plasma arc • 30+ meshes Image credit: Kyle Mckay @ Illinois (joining LLNL) •
Why Blue Waters? • Sparse matrix operations are communication dominant — performance models play a key role. • Exploiting machine layout plays an important role in addressing bottlenecks in communication. • Blue Waters has enabled us to develop/test/scales codes to address these issues Structured Multigrid Algebraic Multigrid
Sparsity (a.k.a. data relationships )
Recommend
More recommend