MRKISS Roadmap

Right now everything that is implemented seems to work. Test coverage is poor for the high level solvers, but solid for the one step solvers and tableaux. I have a few functionality improvements I want to implement, but I don't forsee any API changes. I also have a few examples I would like to add.

• Add dy to a solution without it already
• Example showing error V with t_delta size
• Example showing error vs. order
• Solver that takes a t-vec, and computes solution using those points
• Multi-solver
  • Give it multiple IVs and produce a steps_fixed_* multi-solution
  • Give it multiple PARAMs and produce a steps_fixed_* multi-solution
  • Run parallel across cores using OpenMP
• Finalize step-size and error computation in steps_adapt_etab_*t()
• Add tests for high level solvers
• Hermite Interpolation DONE
  • Given scalar t, y0, y1, dy0, and dy1; compute y(t) CANCELED
  • Given vector t & sol, compute a sol matrix with new y values (no dy) DONE
    • Do linear interpolation if no dy in old solution DONE
  • Add example to three_body.f90 DONE
• Given scalar t & sol; compute y(t) CANCELED
• Compute t-vec with fixed steps. Given vector for storage, and two of: from, to, t_delta DONE
• Solver that takes an t_end and n. Computes single solution at t_end using n fixed sized steps DONE
• Consider how to use examples as consistency tests DONE
• Add dy to solver output DONE
• Add ARKODE_VERNER_13_7_8 DONE
  • BT: https://sundials.readthedocs.io/en/v7.4.0/arkode/Butcher_link.html#c.ARKODE_VERNER_13_7_8 DONE
  • Added ref for Verner's site. DONE
• Add ARKODE_VERNER_10_6_7 DONE
  • BT: https://sundials.readthedocs.io/en/v7.4.0/arkode/Butcher_link.html#c.ARKODE_VERNER_10_6_7 DONE
  • Added ref for Verner's site. DONE
• Create sloppy version of steps_condy_stab_*t called steps_sloppy_condy_stab_*t DONE
• Update comments on steps_* solvers to reflect all options and algorithm. DONE
• Add t_min_o and t_max_o options to print_solution() DONE
• Add example documentation pages DONE
  • three_body DONE
  • lorenz DONE
• Add aspect ratio argument to steps_condy_stab_*t() CANCELED*
• Add "sloppy" option to steps_condy_stab_*t CANCELED*
• Add option to not error out on bisection problems in steps_condy_stab_*t() DONE
• Add option to not error out on bisection problems in steps_adapt_etab_*t() DONE
• Consider "verbosity" options for high level solvers. CANCELED
• Add option to print diagnostics to steps_condy_stab_*t() CANCELED
• Add istats docs to each solver with components used. DONE
• Make comment formatting uniform DONE
• Document all function arguments DONE
• Compare steps_fixed_stab_wt() to other files in rk4.f90 DONE
• Use b1 for stage count in one_step_etab_wt(), and b in in one_step_stab_wt() DONE
• Think about how to document b vectors that require fewer stages than in b DONE
  • mrkiss_eerk_bogacki_shampine_3_2.f90 ..... sd1 1 DONE
  • mrkiss_eerk_dormand_prince_5_4.f90 ....... sd1 1 DONE
  • mrkiss_eerk_fehlberg_7_8.f90 ............. sd1 2 DONE
  • mrkiss_eerk_heun_euler_2_1.f90 ........... sd2 1 DONE
  • mrkiss_eerk_sofroniou_spaletta_4_3.f90 ... sd2 1 DONE
  • mrkiss_eerk_verner_1978_6_5.f90 .......... sd2 2 DONE
• Double Check Order Values DONE
  • mrkiss_eerk_bogacki_shampine_3_2.f90 DONE
  • mrkiss_eerk_bogacki_shampine_4_5.f90 DONE
  • mrkiss_eerk_cash_karp_5_4.f90 DONE
  • mrkiss_eerk_dormand_prince_5_4.f90 DONE
  • mrkiss_eerk_dormand_prince_7_8.f90 DONE
  • mrkiss_eerk_fehlberg_4_5.f90 DONE
  • mrkiss_eerk_fehlberg_7_8.f90 DONE
  • mrkiss_eerk_heun_euler_2_1.f90 DONE
  • mrkiss_eerk_sofroniou_spaletta_4_3.f90 DONE
  • mrkiss_eerk_tsitouras_arkode_5_4.f90 DONE
  • mrkiss_eerk_verner_1978_6_5.f90 DONE
  • mrkiss_eerk_verner_2010_6_5.f90 DONE
  • mrkiss_eerk_verner_9_8.f90 DONE
  • mrkiss_erk_euler_1.f90 DONE
  • mrkiss_erk_feagin_10.f90 DONE
  • mrkiss_erk_knoth_wolke_3.f90 DONE
  • mrkiss_erk_kutta_4.f90 DONE
  • mrkiss_erk_kutta_three_eight_4.f90 DONE
  • mrkiss_erk_midpoint_2.f90 DONE
  • mrkiss_erk_ralston_2.f90 DONE
  • mrkiss_erk_ralston_3.f90 DONE
  • mrkiss_erk_ralston_4.f90 DONE