Description

Fixes an issue where updating the slider of a quickplot could produce NaN values at the beginning or end of a timespan due to slight extrapolation because of floating point arithmetic issues.

For example, try running this simulation and moving the slider to the final time:

import pybamm

model = pybamm.lithium_ion.SPM()
parameter_values = pybamm.ParameterValues("Chen2020")
parameter_values.update({"Current function [A]": 0})

solver = pybamm.IDAKLUSolver()
sim = pybamm.Simulation(model, parameter_values=parameter_values, solver=solver)

tf = 3600.04

sol = sim.solve([0, tf])

# This breaks because of issues with floating point precision
# when we normalize and rescale the time
# tf < (tf / 3600.0) * 3600.0

sol.plot()

Type of change

Please add a line in the relevant section of CHANGELOG.md to document the change (include PR #)

Important checks:

Please confirm the following before marking the PR as ready for review:

• No style issues: nox -s pre-commit
• All tests pass: nox -s tests
• The documentation builds: nox -s doctests
• Code is commented for hard-to-understand areas
• Tests added that prove fix is effective or that feature works