RK4

 

RKF

 

TR-

BDF2

Rosen-

brock

IMEX

 

JVODE

ADAMS

JVODE

BDF

JLSODE

 

Protocol 1 (single pulse)

Tol1

1.0

5.14e-4

1.83e-3

5.24e-4

5.21e-4

5.41e-4

8.36e-4

5.25e-4

Tol2

1.0

5.36e-4

3.93e-3

5.84e-4

5.63e-4

8.14e-4

9.54e-3

7.2e-4

ΔP1 Tol1/Tol2

-

2.2e-5

2.1e-3

6e-5

4.2e-5

2.73e-4

8.7e-3

1.95e-4

Protocol 2 (repeated pulses)

Tol1

1.0

5.29e-4

4.55e-3

5.64e-4

5.57e-4

6.08e-4

1.23e-3

6.28e-4

Tol2

1.0

5.79e-4

9.06e-3

7.1e-4

6.58e-4

1.13e-3

1.63e-2

1.03e-3

ΔP2 Tol1/Tol2

-

2.1e-5

9.41e-3

5.06e-4

1.99e-4

1.8e-3

1.47e-2

1.1e-3

ΔTol1 P1/P2

-

1.5e-5

2.72e-3

4e-5

3.6e-5

6.7e-5

3.9e-4

1.03e-4

ΔTol2 P1/P2

-

4.3e-5

5.13e-3

1.26e-4

9.5e-5

8.14e-4

6.76e-3

3.8e-4
Tol1 and Tol2 represent the ||p|| norm value for each algorithm with P1 and P2 protocols. ΔP1 Tol1 / Tol2 represents the absolute difference between Tol1 and Tol2 for P1, whereas ΔP2 Tol1/ Tol2 represents the same absolute difference for P2. ΔTol1 P1/P2 represents the absolute difference between P1 and P2 for Tol1 parameters, whereas ΔTol2 P1/P2 represents the same absolute difference for Tol2
Table 4: Quantification and comparison of normalized solver performances for the NMDA7 model.