Commit 7f460245 authored by Florian Goth's avatar Florian Goth
Browse files

add the hermitian sixth order method

parent 8dc6f672
Pipeline #7475 passed with stage
in 1 minute and 5 seconds
......@@ -523,11 +523,16 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
! The hermitian fifth order method
call createsplit(splits(40), 10, 10)
splits(40)%Tcoeffs(1)=CMPLX(0.048475520387300861784614942150004872732496275960676605704430, 0.004320853677325454041666651926455479652308652682908398349854853, dp)
splits(40)%Tcoeffs(2)=CMPLX(0.150635519695238479617295632454033652021145596960357218146613, -0.06666435676733963925962676292404991504857683461279976061086538, dp)
splits(40)%Tcoeffs(3)=CMPLX(0.062114931585048261469230020641774248967970727006873004431295, 0.1283252146778869882406120485156541698482809759157383070985385043, dp)
splits(40)%Tcoeffs(4)=CMPLX(0.03942548316294538780808735118455718625805068025688477432907649315, -0.1513147050159361009137685144878063897080171204514607870202000492, dp)
splits(40)%Tcoeffs(5)=CMPLX(0.1993485451694670093207720535696300400203367198152083973885855424, 0.09111504040429731652115356386419287657273864797837757794153296069, dp)
splits(40)%Tcoeffs(1)=CMPLX(0.048475520387300861784614942150004872732496275960676605704430, &
& 0.004320853677325454041666651926455479652308652682908398349854853, dp)
splits(40)%Tcoeffs(2)=CMPLX(0.150635519695238479617295632454033652021145596960357218146613, &
& -0.06666435676733963925962676292404991504857683461279976061086538, dp)
splits(40)%Tcoeffs(3)=CMPLX(0.062114931585048261469230020641774248967970727006873004431295, &
& 0.1283252146778869882406120485156541698482809759157383070985385043, dp)
splits(40)%Tcoeffs(4)=CMPLX(0.03942548316294538780808735118455718625805068025688477432907649315, &
& -0.1513147050159361009137685144878063897080171204514607870202000492, dp)
splits(40)%Tcoeffs(5)=CMPLX(0.1993485451694670093207720535696300400203367198152083973885855424, &
& 0.09111504040429731652115356386419287657273864797837757794153296069, dp)
splits(40)%Tcoeffs(6)=conjg(splits(40)%Tcoeffs(5))
splits(40)%Tcoeffs(7)=conjg(splits(40)%Tcoeffs(4))
splits(40)%Tcoeffs(8)=conjg(splits(40)%Tcoeffs(3))
......@@ -550,23 +555,50 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
! The hermitian sixth order method
call createsplit(splits(41), 16, 16)
splits(30)%Tcoeffs(1)=CMPLX(
splits(30)%Tcoeffs(2)=CMPLX(
splits(30)%Tcoeffs(3)=CMPLX(
splits(30)%Tcoeffs(4)=CMPLX(
splits(30)%Tcoeffs(5)=CMPLX(
splits(30)%Tcoeffs(6)=CMPLX(
splits(30)%Tcoeffs(7)=CMPLX(
splits(30)%Tcoeffs(8)=CMPLX(
splits(30)%Tcoeffs(9)=CMPLX(
splits(30)%Tcoeffs(10) = splits(30)%Tcoeffs(8)
splits(30)%Tcoeffs(11) = splits(30)%Tcoeffs(7)
splits(30)%Tcoeffs(12) = splits(30)%Tcoeffs(6)
splits(30)%Tcoeffs(13) = splits(30)%Tcoeffs(5)
splits(30)%Tcoeffs(14) = splits(30)%Tcoeffs(4)
splits(30)%Tcoeffs(15) = splits(30)%Tcoeffs(3)
splits(30)%Tcoeffs(16) = splits(30)%Tcoeffs(2)
splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
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(9) = conjg(splits(41)%Tcoeffs(8))
splits(41)%Tcoeffs(10) = conjg(splits(41)%Tcoeffs(7))
splits(41)%Tcoeffs(11) = conjg(splits(41)%Tcoeffs(6))
splits(41)%Tcoeffs(12) = conjg(splits(41)%Tcoeffs(5))
splits(41)%Tcoeffs(13) = conjg(splits(41)%Tcoeffs(4))
splits(41)%Tcoeffs(14) = conjg(splits(41)%Tcoeffs(3))
splits(41)%Tcoeffs(15) = conjg(splits(41)%Tcoeffs(2))
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(4)+splits(41)%Vcoeffs(3)+splits(41)%Vcoeffs(2))
splits(41)%Vcoeffs(10)= splits(40)%Vcoeffs(8)
splits(41)%Vcoeffs(11)= splits(40)%Vcoeffs(7)
splits(41)%Vcoeffs(12)= splits(40)%Vcoeffs(6)
splits(41)%Vcoeffs(13)= splits(40)%Vcoeffs(5)
splits(41)%Vcoeffs(14)= splits(40)%Vcoeffs(4)
splits(41)%Vcoeffs(15)= splits(40)%Vcoeffs(3)
splits(41)%Vcoeffs(16)= splits(40)%Vcoeffs(2)
end subroutine
end module splittings
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment