refactor: splits supporting EIS methods from BaseModel, combines simu…

…late() and simulateEIS(), general BaseModel refactoring

examples/scripts/battery_parameterisation/eis_fitting.py

@@ -16,13 +16,14 @@
 )
 
 # Create synthetic data for parameter inference
-sim = model.simulateEIS(
+sim = model.simulate(
     inputs={
         "Negative electrode active material volume fraction": 0.531,
         "Positive electrode active material volume fraction": 0.732,
     },
-    f_eval=f_eval,
+    eval=f_eval,
     initial_state=initial_state,
+    eis=True,
 )
 
 # Fitting parameters

examples/scripts/comparison_examples/grouped_SPMe.py

@@ -99,7 +99,7 @@
     impedances = 1j * np.zeros((Nfreq, NSOC))
     for ii, SOC in enumerate(SOCs):
         model.set_initial_state({"Initial SoC": SOC})
-        simulation = model.simulateEIS(inputs=None, f_eval=frequencies)
+        simulation = model.simulate(inputs=None, eval=frequencies, eis=True)
         impedances[:, ii] = simulation["Impedance"]
 
         if i == 0:

examples/scripts/comparison_examples/jaxified-idaklu-benchmarks.py

@@ -64,7 +64,7 @@ def inputs():
 
     start_time = time.time()
     for _i in range(n):
-        out = problem.model.simulate(inputs=inputs(), t_eval=t_eval)
+        out = problem.model.simulate(inputs=inputs(), eval=t_eval)
     print(f"({solver.name}) Time model.simulate: {time.time() - start_time:.4f}")
 
     start_time = time.time()