#492 change solver times in examples

Author: valentinsulzer

examples/scripts/SPM_compare_particle_grid.py

@@ -49,7 +49,7 @@
 
 # solve model
 solutions = [None] * len(models)
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 for i, model in enumerate(models):
     solutions[i] = model.default_solver.solve(model, t_eval)
 

examples/scripts/compare-dae-solver.py

@@ -24,7 +24,7 @@
 disc.process_model(model)
 
 # solve model
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 
 casadi_sol = pybamm.CasadiSolver(atol=1e-8, rtol=1e-8).solve(model, t_eval)
 klu_sol = pybamm.IDAKLUSolver(atol=1e-8, rtol=1e-8).solve(model, t_eval)

examples/scripts/compare_SPM_diffusion_models.py

@@ -42,7 +42,7 @@
 
 # solve model
 solutions = [None] * len(models)
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 for i, model in enumerate(models):
     solutions[i] = model.default_solver.solve(model, t_eval)
 

examples/scripts/compare_lithium_ion.py

@@ -16,7 +16,7 @@
     pybamm.set_logging_level("INFO")
 
 # load models
-options = {"thermal": "isothermal", "operating mode": "power"}
+options = {"thermal": "isothermal"}
 models = [
     pybamm.lithium_ion.SPM(options),
     pybamm.lithium_ion.SPMe(options),

examples/scripts/thermal_lithium_ion.py

@@ -38,7 +38,7 @@
 
 # solve model
 solutions = [None] * len(models)
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 for i, model in enumerate(models):
     solver = pybamm.ScipySolver(atol=1e-8, rtol=1e-8)
     solution = solver.solve(model, t_eval)

results/change_settings/change_solver_tolerances.py

@@ -27,7 +27,7 @@
 voltages = [None] * len(tols)
 voltage_rmse = [None] * len(tols)
 labels = [None] * len(tols)
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 for i, tol in enumerate(tols):
     solver = pybamm.ScikitsDaeSolver(rtol=tol[0], atol=tol[1])
     solutions[i] = solver.solve(model, t_eval)

results/change_settings/compare_var_pts.py

@@ -43,7 +43,7 @@
 solutions = [None] * len(models)
 voltages = [None] * len(models)
 voltage_rmse = [None] * len(models)
-t_eval = np.linspace(0, 0.17, 100)
+t_eval = np.linspace(0, 0.25, 100)
 for i, model in enumerate(models):
     solutions[i] = model.default_solver.solve(model, t_eval)
     voltages[i] = pybamm.ProcessedVariable(

tests/integration/test_solvers/test_idaklu.py

@@ -50,7 +50,7 @@ def test_changing_grid(self):
 
         # Calculate time for each solver and each number of grid points
         var = pybamm.standard_spatial_vars
-        t_eval = np.linspace(0, 0.17, 100)
+        t_eval = np.linspace(0, 0.25, 100)
         for npts in [100, 200]:
             # discretise
             var_pts = {