Construction and comparison of parallel implicit kinetic solvers in three spatial dimensions