Commit 0f8504e2 authored by Florian Goth's avatar Florian Goth
Browse files

bugfix

parent 9a90b8b6
Pipeline #6959 passed with stage
in 1 minute and 22 seconds
...@@ -154,15 +154,15 @@ subroutine init_splittings(split) ...@@ -154,15 +154,15 @@ subroutine init_splittings(split)
splits(9)%Vcoeffs(2) = 2.D0/5 splits(9)%Vcoeffs(2) = 2.D0/5
splits(9)%Vcoeffs(3) = -1.D0/10 splits(9)%Vcoeffs(3) = -1.D0/10
splits(9)%Vcoeffs(4) = 1.D0 - 2*(splits(9)%Vcoeffs(1)+splits(9)%Vcoeffs(2)) splits(9)%Vcoeffs(4) = 1.D0 - 2*(splits(9)%Vcoeffs(1)+splits(9)%Vcoeffs(2))
splits(9)%Vcoeffs(5) = splits(9)%Vcoeffs(2) splits(9)%Vcoeffs(5) = splits(9)%Vcoeffs(3)
splits(9)%Vcoeffs(6) = splits(9)%Vcoeffs(1) splits(9)%Vcoeffs(6) = splits(9)%Vcoeffs(2)
! Symmetric from symmetric ! Symmetric from symmetric
! This coresponds to the classical Forest/Ruth algorithm ! This corresponds to the classical Forest/Ruth algorithm
call createsplit(splits(10), 4, 4) call createsplit(splits(10), 4, 4)
alpha = 1/(2-2**(1.D0/3.D0) ) alpha = 1/(2-2**(1.D0/3.D0) )
beta = 1-2*alpha beta = 1.D0-2.D0*alpha
splits(10)%Tcoeffs(1) = alpha/2 splits(10)%Tcoeffs(1) = alpha/2
splits(10)%Tcoeffs(2) = (alpha + beta)/2 splits(10)%Tcoeffs(2) = (alpha + beta)/2
splits(10)%Tcoeffs(3) = splits(10)%Tcoeffs(2) splits(10)%Tcoeffs(3) = splits(10)%Tcoeffs(2)
...@@ -188,8 +188,8 @@ subroutine init_splittings(split) ...@@ -188,8 +188,8 @@ subroutine init_splittings(split)
splits(11)%Vcoeffs(2) = alpha splits(11)%Vcoeffs(2) = alpha
splits(11)%Vcoeffs(3) = alpha splits(11)%Vcoeffs(3) = alpha
splits(11)%Vcoeffs(4) = beta splits(11)%Vcoeffs(4) = beta
splits(11)%Vcoeffs(5) = splits(11)%Vcoeffs(2) splits(11)%Vcoeffs(5) = splits(11)%Vcoeffs(3)
splits(11)%Vcoeffs(6) = splits(11)%Vcoeffs(1) splits(11)%Vcoeffs(6) = splits(11)%Vcoeffs(2)
! fourth order method from Omelyan, 2001, previously discovered by McLachlan 1995 ! fourth order method from Omelyan, 2001, previously discovered by McLachlan 1995
......
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