SuperLU (Sparse Direct Solver)

More documents

Recommendations

Info

SuperLU_DIST Example Program (Fortran 90)• SuperLU_DIST_2.0/FORTRAN/• All SuperLU objects (e.g., LU structure) are opaque for F90• They are allocated, deallocated and operated in the C side and notdirectly accessible from Fortran side.• C objects are accessed via handles that exist in Fortran’suser space• In Fortran, all handles are of type INTEGERShort Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 200520
f_pddrive.f90 (1/2)program f_pddriveuse superlu_modinclude 'mpif.h'! Declarationsinteger(superlu_ptr) :: gridinteger(superlu_ptr) :: optionsinteger(superlu_ptr) :: ScalePermstructinteger(superlu_ptr) :: LUstructinteger(superlu_ptr) :: SOLVEstructinteger(superlu_ptr) :: Ainteger(superlu_ptr) :: stat! Create Fortran handles for the C structures used inSuperLU_DISTcall f_create_gridinfo(grid)call f_create_options(options)call f_create_ScalePermstruct(ScalePermstruct)call f_create_LUstruct(LUstruct)call f_create_SOLVEstruct(SOLVEstruct)call f_create_SuperMatrix(A)call f_create_SuperLUStat(stat)! Initialize MPI environmentcall mpi_init(ierr)! Set up the distributed input matrix Acall f_dcreate_dist_matrix(A, m, n, nnz, values,rowind, colptr, grid)! Set the default solver optionscall f_set_default_options(options)! Initialize ScalePermstruct and LUstructcall get_SuperMatrix(A,nrow=m,ncol=n)call f_ScalePermstructInit(m, n, ScalePermstruct)call f_LUstructInit(m, n, LUstruct)! Initialize the statistics variablescall f_PStatInit(stat)! Call the linear equation solvercall f_pdgssvx(options, A, ScalePermstruct, b,ldb, nrhs, grid, LUstruct, SOLVEstruct, berr,stat, info)Short Course on the DOE ACTS Collection - SIAM CSE05 ConferenceOrlando, FL - February 11, 200521
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 and 16: 2D Block Cyclic Layout for L and U
Page 17 and 18: SuperLU_DIST Example Program (C)•
Page 19: Pddrive.c (2/2)/* Initialize ScaleP
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?