DanHickstein

import numpy as np
import abel
import matplotlib.pyplot as plt
import time

# nps = np.logspace(1,9,12)
nps = [10,20,2e9]

print nps
times1 = []

DanHickstein

import numpy as np
import matplotlib.pyplot as plt
import abel

sizes = np.logspace(np.log10(50), np.log10(1e4), 20)
n_max_bs   = 2001
n_max_slow = 2001
transform_repeat = 4 # for less than 400

results = {'basex' :         [], 

DanHickstein

import numpy as np
import matplotlib.pyplot as plt
from scipy.linalg import circulant
from scipy.optimize import curve_fit, brentq
from scipy.interpolate import interp1d
from scipy.ndimage import gaussian_filter
import scipy.ndimage as nd

def gaussian(x, a, mu, sigma, c):
    """

DanHickstein

for i in range(1,10):
    print i