Use -pc_type lu -pc_factor_mat_solver_type elemental to use this direct solver
-mat_type elemental | - sets the matrix type to "elemental" during a call to MatSetFromOptions() | |
-mat_elemental_grid_height | - sets Grid Height for 2D cyclic ordering of internal matrix |