def METHOD_NAME(self, plot):
    mod_current = 20e6
    self.optional_init_solver.set_profiles(
        self.p_prime, self.ff_prime, I_p=mod_current
    )
    self.optional_init_solver.set_mesh(self.mesh)
    self.optional_init_solver.solve(plot=plot)
    assert np.isclose(self.optional_init_solver._curr_target, mod_current)