cwyenberg

import numpy as np

def auto_corr(data, normalize=True):

    # Performs an autocorrelation on data
    # Note that for data size n_data, we autocorrelate to 2*n_data-1 points
    # (i.e. this returns a larger array than passed)

    n_data = data.shape[0]          # Number of data points in the input data
    n_ac = 2*n_data-1               # Autocorrelation array size