fix method 6 and 23

Prog/splittings.f90

@@ -103,14 +103,14 @@ subroutine init_splittings(split)
     splits(6)%Tcoeffs(7) = splits(6)%Tcoeffs(1)
 
 
-
-    splits(6)%Vcoeffs(1) = 0.209515106613362
-    splits(6)%Vcoeffs(2) = - 0.143851773179818
-    splits(6)%Vcoeffs(3) = (1-2.D0*(splits(6)%Vcoeffs(1)+splits(6)%Vcoeffs(2) ))/2.D0
-    splits(6)%Vcoeffs(4) = splits(6)%Vcoeffs(3)
-    splits(6)%Vcoeffs(5) = splits(6)%Vcoeffs(2)
-    splits(6)%Vcoeffs(6) = splits(6)%Vcoeffs(1)
-    splits(6)%Vcoeffs(7) = 0.0
+    splits(6)%Vcoeffs(1) = 0.0
+    splits(6)%Vcoeffs(2) = 0.209515106613362
+    splits(6)%Vcoeffs(3) = - 0.143851773179818
+    splits(6)%Vcoeffs(4) = (1-2.D0*(splits(6)%Vcoeffs(2)+splits(6)%Vcoeffs(3) ))/2.D0
+    splits(6)%Vcoeffs(5) = splits(6)%Vcoeffs(4)
+    splits(6)%Vcoeffs(6) = splits(6)%Vcoeffs(3)
+    splits(6)%Vcoeffs(7) = splits(6)%Vcoeffs(2)
+!    splits(6)%Vcoeffs(8) = 0.0
 
     ! S_10 4 from https://arxiv.org/pdf/hep-lat/0505020.pdf
     call createsplit(splits(7), 5, 5)
@@ -253,13 +253,13 @@ subroutine init_splittings(split)
     enddo
 
     ! the variant with m=2 optimized for large step sizes in Blanes 2014, 2. order
-    call createsplit(splits(23), 3, 2)
+    call createsplit(splits(23), 3, 3)
     splits(23)%Tcoeffs(1) = 0.21178
     splits(23)%Tcoeffs(2) = (1.0 - 2.0*splits(23)%Tcoeffs(1))
     splits(23)%Tcoeffs(3) = splits(23)%Tcoeffs(1)
 
-    splits(23)%Vcoeffs(1) = 0.5
-    splits(23)%Vcoeffs(2) = 0.5
+    splits(23)%Vcoeffs = 0.5
+    splits(23)%Vcoeffs(1) = 0.0
 
     ! the variant with m=3 optimized for large step sizes in Blanes 2014, 2. Order
     call createsplit(splits(24), 4, 4)