cvals = np.logspace(-1, 3, 20)