A fast and robust numerical scheme for solving models of charge carrier transport and ion vacancy motion in perovskite solar cells