Many extreme-scale simulation codes encompass multiphysics components in multiple spatial and length scales. The resulting discretized sparse linear systems can be highly indefinite, nonsymmetric and extremely ill-conditioned. For such problems, factorization based algorithms are often the most robust algorithmic choices among many alternatives. We present our recent research on novel parallel factorization algorithms that are efficient for solving such problems. From algorithm side, we incorporate data-sparse low-rank structures, such as hierarchical matrix algebra, to achieve lower arithmetic and communication complexity as well as robust preconditioner. From parallelization side, we exploit sparse data structures represented by DAGs and trees to schedule coarse-grained tasks and SIMD/SIMT for fine-grained parallelism. We will illustrate both theoretical and practical aspects of the methods, and demonstrate their performance on manycore architectures including GPU clusters and the latest Intel Xeon Phi KNL platforms, using a variety of real world problems.
Sherry Li is a Senior Scientist at Lawrence Berkeley National Laboratory. She has worked on diverse problems in high performance scientific computations, including parallel computing, sparse matrix computations, high precision arithmetic, and combinatorial scientific computing. She has (co)authored over 90 publications. She is the lead developer of SuperLU sparse direct solver library, and has contributed to several other widely-used mathematical libraries, including LAPACK, ARPREC, XBLAS, and STRUMPACK. She received Ph.D. in Computer Science from UC Berkeley in 1996. She is a SIAM Fellow and an ACM Senior Member.