def test_output_shapes(method, shape, data):
    t, y, dy = data
    freq = np.asarray(np.random.rand(*shape))
    freq.flat = np.arange(1, freq.size + 1)
    PLS = lombscargle(t, y, frequency=freq,
                      fit_bias=False, method=method)
    assert_equal(PLS.shape, shape)