Use the new sixth order method with all positive real parts

Prog/splittings.f90

@@ -554,24 +554,24 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
     
     
     
-    ! CHR_15 6 The hermitian sixth order method
+    ! CHR_15 6 The hermitian sixth order method, now with all-positive real parts
     call createsplit(splits(41), 16, 16)
-    splits(41)%Tcoeffs(1)=CMPLX(0.095194656759747098611056828381842973876432181037008757210040895, &
-    &-0.03705839851161959800834474467108207397222103451188259084113901385692675627401, dp)
-    splits(41)%Tcoeffs(2)=CMPLX(-0.15104982809645175788409611415108720646496746403230715858946175650322203501371,&
-    & 0.05117453376236495306981471912960643425256400016528840676048190434042291332412, dp)
-    splits(41)%Tcoeffs(3)=CMPLX(0.21444360513379835812496101375826321624154191246294877169099223152148831688044, &
-    & 0.07359732051007663675131606858704201625662656959313205914177936096667505786252, dp)
-    splits(41)%Tcoeffs(4)=CMPLX(-0.37730521516220511037719174514074952136575189147542328395542965222449129781211, &
-    & -0.36816955377298596235653518144042308597642054664893410445462793958605313278742, dp)
-    splits(41)%Tcoeffs(5)=CMPLX(0.40447858431312636848006975094315767332042173322914650526876550339201968853901, &
-    & 0.16787162132465120731961550607265408536394314127464993116585034439506196635885, dp)
-    splits(41)%Tcoeffs(6)=CMPLX(0.0758332887820294125867303882030046525792465007278686494263456890795286587039, &
-    & 0.19359832437117728499779816875001896381499448865737772505570186568641242918512, dp)
-    splits(41)%Tcoeffs(7)=CMPLX(0.28739421273943439068655447271890074375674329885392002549567039178249224038604, &
-    & -0.07927378003029551600806762944587427553776341998148942145393739830240004259861, dp)
-    splits(41)%Tcoeffs(8)=CMPLX(-0.04898930446947876022808459471333253194366627080316226654692330248307450386446, &
-    & -0.02347347249741907478797873048501010730909553249059523184166077807099858428049, dp)
+    splits(41)%Tcoeffs(1)=CMPLX(0.018407829100474781904003178937442, &
+    &0.060063925504056716987384011897458, dp)
+    splits(41)%Tcoeffs(2)=CMPLX(0.037640136069562677570517610827614,&
+    & -0.078183261287778862724569177453338, dp)
+    splits(41)%Tcoeffs(3)=CMPLX(0.081602637043566199583289018818188, &
+    & 0.044188227706548767875286813464238, dp)
+    splits(41)%Tcoeffs(4)=CMPLX(0.08163394645009567118137080903697, &
+    & -0.10257156270865030554005352335570, dp)
+    splits(41)%Tcoeffs(5)=CMPLX(0.082577837091366965137963780758131, &
+    & 0.077099467183537374713948506168652, dp)
+    splits(41)%Tcoeffs(6)=CMPLX(0.096825373607680378408689859741856, &
+    & 0.035803350489900133270555458990013, dp)
+    splits(41)%Tcoeffs(7)=CMPLX(0.028447864413379336326027072059858, &
+    & -0.087263167058591631392917553485744, dp)
+    splits(41)%Tcoeffs(8)=CMPLX(0.07286437622387398988813866981994, &
+    &0.14749079765627217175324647517686, dp)
 
     splits(41)%Tcoeffs(9)  = conjg(splits(41)%Tcoeffs(8))
     splits(41)%Tcoeffs(10) = conjg(splits(41)%Tcoeffs(7))
@@ -583,14 +583,14 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
     splits(41)%Tcoeffs(16) = conjg(splits(41)%Tcoeffs(1))
 
     splits(41)%Vcoeffs(1) = 0.0
-    splits(41)%Vcoeffs(2) = 0.03671110579787723640996671077037600088884565405304545274883123346646457750091
-    splits(41)%Vcoeffs(3) = 0.04613662497236318249749919552727510611956046532609799759380515271198742087952
-    splits(41)%Vcoeffs(4) = 0.06768585482545818439443581352088950415355030213002216767635230717099400568149
-    splits(41)%Vcoeffs(5) = 0.00073425375666997071669156265943175774559697103109632953812839071428711046668 
-    splits(41)%Vcoeffs(6) = 0.09716293263662686055508652826515692175297846447931632392502750914892324410087
-    splits(41)%Vcoeffs(7) = 0.10813899252679550041628722441016135854566318598566710383062396269076592454722
-    splits(41)%Vcoeffs(8) = 921.0/79000
-    splits(41)%Vcoeffs(9) = 1-2.D0*(splits(41)%Vcoeffs(8) + splits(41)%Vcoeffs(7) + splits(41)%Vcoeffs(6) + splits(41)%Vcoeffs(5) &
+    splits(41)%Vcoeffs(2) = 0.013599999999999999242272785693331
+    splits(41)%Vcoeffs(3) = 0.10132124738923978821039934776422
+    splits(41)%Vcoeffs(4) = 0.08496255121091424240063136564739
+    splits(41)%Vcoeffs(5) = 0.023858373653897122447214292182343
+    splits(41)%Vcoeffs(6) = 0.12203432116901233344250458684336
+    splits(41)%Vcoeffs(7) = 0.09460538089501579698109829739963
+    splits(41)%Vcoeffs(8) = 0.058633436576191442147232738898912
+    splits(41)%Vcoeffs(9) = 1.D0-2.D0*(splits(41)%Vcoeffs(8) + splits(41)%Vcoeffs(7) + splits(41)%Vcoeffs(6) + splits(41)%Vcoeffs(5) &
     &+splits(41)%Vcoeffs(4)+splits(41)%Vcoeffs(3)+splits(41)%Vcoeffs(2))
     splits(41)%Vcoeffs(10)= splits(41)%Vcoeffs(8)
     splits(41)%Vcoeffs(11)= splits(41)%Vcoeffs(7)