Using the Idaklu Solver causing a memory leak

[NicholasBudenberg-CAM]

I am running the following code in order to see how the degradation of the cell is effected by different drive cycles. My drive cycle has spiky currents and I have found that using the Casadi Solver fails unless I use a small timestep which is unreasonable for the other steps. I have tried the IDAKLU solver and this solved that issue but seemed to cause another on.

import pybamm
import numpy as np
import time

ExperimentFolder = f"D:/Stack Exchange"

model = pybamm.lithium_ion.DFN()

### This is a Sinusoidal Discharge Current using the below values
DCOffset = 5
ACAmplitude = 5
ACFrequency = 0.01
NumberOfCycles = 10000
Resolution = 100 # This is the number of points per cycle

timeseries = np.linspace(0, 1 / ACFrequency * NumberOfCycles, NumberOfCycles * Resolution)
current = DCOffset + ACAmplitude * (np.sin(2 * np.pi * timeseries * ACFrequency))
Discharge_Cycle = np.column_stack([timeseries, current])

Steps = [
    pybamm.step.current(Discharge_Cycle, termination="2.5 V"),
    "Charge at 0 A for 5 minutes",
    "Charge at 10 A until 4.2 V",
    "Hold at 4.2 V until 50 mA",
    "Charge at 0 A for 5 minutes"
        ]

experiment = pybamm.Experiment(Steps)
Number_Of_Cycles = 500

for i in range(Number_Of_Cycles):
    Start_Time = time.time()
    # sim = pybamm.Simulation(model, experiment=experiment, solver=pybamm.IDAKLUSolver())
    sim = pybamm.Simulation(model, experiment=experiment, solver=pybamm.CasadiSolver(mode="safe", extra_options_setup={"max_step_size": 1 / (ACFrequency * Resolution)}))
    solution = sim.solve()
    model = model.set_initial_conditions_from(solution.last_state, inplace=False)
    solution.save(f"{ExperimentFolder}/Cycle_{i}.pkl")
    Stop_Time = time.time()
    print(f"Performed Cycle {i} out of {Number_Of_Cycles}, taking {Stop_Time - Start_Time}s to solve the last ### Cycle")

Here you can see the Casadi solver:

however when using the IDAKLU solver there seems to be a memory leak, however the simulation runs a lot more smoothly:

Any help as to what I am doing wrong would be much appreciated (or other ways of accomplishing what I am trying to do).

[github-actions[bot]]

Hi! We have now moved our discussions to Discourse. Please post your question there.

[NicholasBudenberg-CAM]

This can now be found here: https://pybamm.discourse.group/t/using-the-idaklu-solver-causing-a-memory-leak/117/4