def simulate(a, xi0, theta0, phi0, H, idx):
sim = Simulator(a, xi0, dirichlet=True)
sim.simulate(theta0, phi0, H)
U = sim.conservation()
R = sim.reflectivity()
print(idx)
return [U, R]
def simulate(a, xi0, theta0, phi0, H, idx):
sim = Simulator(a, xi0, dirichlet=False)
r = sim.simulate(theta0, phi0, H)
U = sim.conservation()
R = sim.reflectivity()
return [U, R, r]
