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]