def run(A, choen=False):
V_x = barrier(x, A, x1, x2)
sch = Schrodinger(x, Psi_x, V_x, hbar=hbar, m=m, args=x2, t=0)
I = sch.impedencePacket()
while sch.t < finalt:
print("Height :{h} Time:{t}".format(h=A, t=sch.t))
sch.evolve_t(step, dt)
T = sch.barrier_transmition()
return T, I
def run_A(A):
Psi_x = gauss_init(x, k0, x0=x0, d=sig)
V_x = gauss_barrier(x, A, x1, omeg)
sch = Schrodinger(x, Psi_x, V_x, hbar=hbar, m=m, t=0, args=x1)
imp = sch.impedencePacket(tol=10**-11)
time_list = []
Trans_list = []
while sch.t < t_final:
sch.evolve_t(N_steps=step, dt=dt)
time_list.append(sch.t)
Trans_list.append(sch.barrier_transmition())
print("Height {v}, Time {t}".format(v=A / scale, t=sch.t))
return sch.barrier_transmition(), imp