log.info("Simulation configuration complete.") for i,system in enumerate(systems): #Construct Actions actionList = [] if args.i: actionList += [actions.Info()] if display_output: actionList += [gp_plotter.Plotter2D( *gnu_plot_settings, frequency = 1, system = system, delay = 0. )] if store_output: actionList += [actions.SimOutput( hdf_file, RKsolvers[i], system, grids[i], output_actions, overwrite = True, name = grids[i].name, cmp_ = grids[i].comparison )] log.info("Starting simulation %i with system %s"%(i,repr(system))) problem = ibvp.IBVP(RKsolvers[i], system, grid = grids[i],\ maxIteration = 1000000, action = actionList) problem.run(tstart, tstop) log.info("Simulation complete") log.info("Calculations complete")
from coffee import ibvp, solvers, grid from coffee.actions import gp_plotter from OneDWave import OneDwave system = OneDwave() solver = solvers.RungeKutta4(system) grid = grid.UniformCart((200,), [(0, 4)]) plotter = gp_plotter.Plotter1D( system, 'set terminal qt', 'set yrange [-1:1]', 'set style data lines', ) problem = ibvp.IBVP(solver, system, grid=grid, action=[plotter]) problem.run(0, 3)