add 8th order McLachlan extended precision coefficients

LLNL · Jun 28, 2023 · 38dc866 · 38dc866
1 parent 3795041
commit 38dc866
Showing 1 changed file with 8 additions and 8 deletions.
diff --git a/src/arkode/arkode_sprk.c b/src/arkode/arkode_sprk.c
@@ -247,14 +247,14 @@ ARKodeSPRKStorage ARKodeSymplecticMcLachlan8()
   ARKodeSPRKStorage sprk_storage = ARKodeSPRKStorage_Alloc(16);
   sprk_storage->q                = 8;
   sprk_storage->stages           = 16;
-  sprk_storage->a[0]             = SUN_RCONST(0.74167036435061295344822780);
-  sprk_storage->a[1]             = -SUN_RCONST(0.40910082580003159399730010);
-  sprk_storage->a[2]             = SUN_RCONST(0.19075471029623837995387626);
-  sprk_storage->a[3]             = -SUN_RCONST(0.57386247111608226665638773);
-  sprk_storage->a[4]             = SUN_RCONST(0.29906418130365592384446354);
-  sprk_storage->a[5]             = SUN_RCONST(0.33462491824529818378495798);
-  sprk_storage->a[6]             = SUN_RCONST(0.31529309239676659663205666);
-  sprk_storage->a[7]             = -SUN_RCONST(0.79688793935291635401978884);
+  sprk_storage->a[0]             = SUN_RCONST(0.7416703643506129534482278017838063156035);
+  sprk_storage->a[1]             = -SUN_RCONST(0.4091008258000315939973000958935634173099);
+  sprk_storage->a[2]             = SUN_RCONST(0.1907547102962383799538762564503716627355);
+  sprk_storage->a[3]             = -SUN_RCONST(0.5738624711160822666563877266355357421595);
+  sprk_storage->a[4]             = SUN_RCONST(0.2990641813036559238444635406886029882258);
+  sprk_storage->a[5]             = SUN_RCONST(0.3346249182452981837849579798821822886337);
+  sprk_storage->a[6]             = SUN_RCONST(0.3152930923967665966320566638110024309941);
+  sprk_storage->a[7]             = -SUN_RCONST(0.7968879393529163540197888401737330534463);
   sprk_storage->a[8]             = sprk_storage->a[6];
   sprk_storage->a[9]             = sprk_storage->a[5];
   sprk_storage->a[10]            = sprk_storage->a[4];