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))
