Use -pc_type lu -pc_factor_mat_solver_type scalapack to use this direct solver
| -mat_type scalapack | - sets the matrix type to "scalapack" during a call to MatSetFromOptions() | |
| -mat_scalapack_grid_height | - sets Grid Height for 2D cyclic ordering of internal matrix | |
| -mat_scalapack_block_sizes | - size of the blocks to use (one or two integers separated by comma) |