@set_fft_params<EOL>def psd(timeseries, method_func, *args, **kwargs):
<EOL>return _psdn(timeseries, method_func, *args, **kwargs)<EOL>
Generate a PSD using a method function

    All arguments are presumed to be given in physical units

    Parameters
    ----------
    timeseries : `~gwpy.timeseries.TimeSeries`, `tuple`
        the data to process, or a 2-tuple of series to correlate

    method_func : `callable`
        the function that will be called to perform the signal processing

    *args, **kwargs
        other arguments to pass to ``method_func`` when calling