Adaptive IFT of arbitrary Transfer-Function

I wish to compute the Inverse Fourier Transform of an arbitrary transfer function F(frequency) of a dynamical system to obtain its response in the time-domain F(time) where F(time) must meet certain absolute and relative error tolerances. Additionally, evaluating F(frequency) is relatively expensive so I wish to perform the IFT in an adaptive way to minimize the number of evaluations of F(frequency). If there an existing function in matlab to perform this operation? Something along the lines of:

F_time = function(@F_freq, t_vec, abstol, reltol);