SuperLU (Sparse Direct Solver)

More documents

Recommendations

Info

Outline• Overview of the software• Some background of the algorithms• <strong>Sparse</strong> matrix distribution and user interface• Example program (Fortran 90 interface)• ApplicationsShort Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 20052
What is <strong>SuperLU</strong>• Solve general sparse linear system A x = b.• Example: A of dimension 10 5, only 10 ~ 100 nonzeros per row• Algorithm: Gaussian elimination (LU factorization: A =LU), followed by lower/upper triangular solutions.• Store only nonzeros and perform operations only on nonzeros.• Efficient and portable implementation for highperformancearchitectures, flexible interface.Short Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 20053
Page 1: SIAM Conference on Computational Sc
Page 5 and 6: Contents of the SuperLU Library•
Page 7 and 8: Supernode• Exploit dense submatri
Page 9 and 10: Numerical Pivoting• Goal of pivot
Page 11 and 12: Ordering Interface in SuperLU• Li
Page 13 and 14: Numerical Factorization• Sequenti
Page 15 and 16: 2D Block Cyclic Layout for L and U
Page 17 and 18: SuperLU_DIST Example Program (C)•
Page 19 and 20: Pddrive.c (2/2)/* Initialize ScaleP
Page 21 and 22: f_pddrive.f90 (1/2)program f_pddriv
Page 23 and 24: Applications• Used to solve open
Page 25: Summary• Efficient implementation

SuperLU (Sparse Direct Solver)

Create successful ePaper yourself

Delete template?

Save as template?