def test_compare_with_mathematica2(): mathematica = [-8.20749E-8, -3.17403E-8, -5.23501E-9, -5.90173E-10,\ 2.90532E-10, 6.28546E-10, 1.30025E-9,3.02873E-9] vb = [-1, -.75, -.5, -.25, .25, .50, .75, 1] for i in range(len(vb)): f = models.tunnelmodel_1level_nogate_300K_gauss(vb[i],.0005, 0.02,\ 0.62, 0.7, 0.14) assert abs((f - mathematica[i]) / f) < 0.05
def fxn(vb, gammaL, gammaR, deltaE, eta, gauss): n = 1 I = models.tunnelmodel_1level_nogate_300K_gauss(vb, gammaL, gammaR, deltaE, eta, gauss) return np.log10(abs((n * I) / 2.41E-06))