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)