from hopf.simulation import Simulation # Sphere integrator s = Simulation() s.load_initial_conditions('collapse3.mat') # Medium regularization s.sigma = 0.10 s.run_simulation(tmax=500, numpoints=500, sim='sphere-midpoint') s.post_process() s.save_results('data/collapse3_sphere_sigma10_sim_long.mat') # Runge-Kutta, medium regularization s.sigma = 0.10 s.run_simulation(tmax=500, numpoints=500, sim='rk4') s.post_process() s.save_results('data/collapse3_rk4_sigma10_sim_long.mat') # Lie-Poisson, medium regularization s.sigma = 0.10 s.run_simulation(tmax=500, numpoints=500, sim='lie-poisson') s.post_process() s.save_results('data/collapse3_lp_sigma10_sim_long.mat') # Midpoint, medium regularization s.sigma = 0.10 s.run_simulation(tmax=500, numpoints=500, sim='midpoint') s.post_process() s.save_results('data/collapse3_mp_sigma10_sim_long.mat')
from hopf.simulation import Simulation # Sphere integrator s = Simulation() s.load_initial_conditions('svs5_poles.mat') s.sigma = 0.25 s.run_simulation(h=0.3, tmax=1000, sim='sphere-midpoint-momentum') s.post_process() s.save_results('data/svs5_poles_sphere_long_025.mat') # RK4 integrator s = Simulation() s.load_initial_conditions('svs5_poles.mat') s.sigma = 0.25 s.run_simulation(h=0.3, tmax=1000, sim='rk4') s.post_process() s.save_results('data/svs5_poles_rk4_long_025.mat') # MP integrator s = Simulation() s.load_initial_conditions('svs5_poles.mat') s.sigma = 0.25 s.run_simulation(h=0.3, tmax=1000, sim='midpoint') s.post_process() s.save_results('data/svs5_poles_mp_long_025.mat') # LP integrator s = Simulation() s.load_initial_conditions('svs5_poles.mat') s.sigma = 0.25 s.run_simulation(h=0.3, tmax=1000, sim='lie-poisson')
from hopf.simulation import Simulation # Sphere integrator s = Simulation() s.load_initial_conditions('collapse3.mat') # No regularization #s.sigma = 0.0 #s.run_simulation(tmax=15, numpoints=75, sim='sphere-midpoint') #s.post_process() #s.save_results('data/collapse3_sphere_sigma00_sim.mat') # Medium regularization s.sigma = 0.10 s.run_simulation(tmax=15, numpoints=75, sim='sphere-midpoint') s.post_process() s.save_results('data/collapse3_sphere_sigma10_sim.mat') # Large regularization s.sigma = 0.25 s.run_simulation(tmax=15, numpoints=75, sim='sphere-midpoint') s.post_process() s.save_results('data/collapse3_sphere_sigma25_sim.mat') # Runge-Kutta, medium regularization s.sigma = 0.10 s.run_simulation(tmax=15, numpoints=75, sim='rk4') s.post_process() s.save_results('data/collapse3_rk4_sigma10_sim.mat') # Lie-Poisson, medium regularization