def fire(coords, pot, tol=1e-3, nsteps=100000, **kwargs): """ A wrapper function for the pele FIRE implementation """ opt = Fire(coords, pot, **kwargs) res = opt.run(fmax=tol, steps=nsteps) return res