SuperLU (Sparse Direct Solver)

More documents

Recommendations

Info

Scalability• 3D KxKxK cubic grids, scale N 2 = K 6 with P for constant work per processor• Achieved 12.5 and 21.2 Gflops on 128 processors• Performance sensitive to communication latency• Cray T3E latency: 3 microseconds ( ~ 2702 flops)• IBM SP latency: 8 microseconds ( ~ 11940 flops )Short Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 200516
SuperLU_DIST Example Program (C)• SuperLU_DIST_2.0/EXAMPLE/pddrive.c• Five basic steps1. Initialize the MPI environment and SuperLU process grid2. Set up the input matrices A and B3. Set the options argument (can modify the default)4. Call SuperLU routine PDGSSVX5. Release the process grid, deallocate memory, and terminate theMPI environmentShort Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 200517
Page 1 and 2: SIAM Conference on Computational Sc
Page 3 and 4: What is SuperLU• Solve general sp
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: 2D Block Cyclic Layout for L and U
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?