In this ECP project we deliver factorization based sparse solvers based on two algorithm variants: supernodal (SuperLU) and multifrontal (STRUMPACK). STRUMPACK is further enhanced with scalable preconditioning functionality using hierarchical matrix algebra. Both libraries are purely algebraic, applicable to a large variety of applications. We address several exascale computing challenges 1) We develop novel approximation algorithms with lower arithmetic and communication complexity 2) We implement new parallelization techniques that reduce communication and improve vectorization 3) We integrate our software into higher level algebraic solvers such as hypre, PETSc, Trilinos, and we collaborate with ECP application teams.