Commit 9194df49 authored by Florian Goth's avatar Florian Goth
Browse files

fix wrong coefficient.

parent 2ee180f4
Pipeline #10544 passed with stage
in 1 minute and 6 seconds
......@@ -47,14 +47,14 @@ subroutine init_splittings(split, myparam)
real(kind = kind(0.D0)), intent(in) :: myparam
integer, parameter :: dp = kind(1.d0)
Integer :: nrofsplits, n, beg, tmpint, i
character(len=11) :: filepath
complex(kind = kind(0.D0)) :: alpha, beta, xi, lam, chi
real(kind = kind(0.D0)) :: ts
real(kind = kind(0.D0)) :: ts, u1,u2,u3,v1,v2,v3,b1,b2
namelist /MET/ u1, u2, u3, v1, v2, v3, b1, b2
chosensplitting = split
nrofsplits = 41
nrofsplits = 42
Allocate(splits(nrofsplits))
! V-T Euler / NS_1 1
call createsplit(splits(1), 1, 1)
splits(1)%Tcoeffs(1) = 1.0
......@@ -559,7 +559,7 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
! 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.018407829100474781904003178937442, &
splits(41)%Tcoeffs(1)=CMPLX(0.054840017212412411776848, &
&0.060063925504056716987384011897458, dp)
splits(41)%Tcoeffs(2)=CMPLX(0.037640136069562677570517610827614,&
& -0.078183261287778862724569177453338, dp)
......@@ -603,6 +603,23 @@ splits(30)%Tcoeffs(17) = splits(30)%Tcoeffs(1)
splits(41)%Vcoeffs(15)= splits(41)%Vcoeffs(3)
splits(41)%Vcoeffs(16)= splits(41)%Vcoeffs(2)
filepath = 'method.cfg'
open(action='read', file=filepath, iostat=istat, newunit=iunit)
read (unit=iunit, nml=MET, iostat=istat)
write (*,*) u1,u2,u3,v1,v2,v3,b1,b2
call createsplit(splits(42), 6, 6)
splits(42)%Tcoeffs(1) = CMPLX(u1, v1, dp)
splits(42)%Tcoeffs(2) = CMPLX(u2, v2, dp)
splits(42)%Tcoeffs(3) = CMPLX(u3, v3, dp)
splits(42)%Tcoeffs(4) = conjg(splits(42)%Tcoeffs(3))
splits(42)%Tcoeffs(5) = conjg(splits(42)%Tcoeffs(2))
splits(42)%Tcoeffs(6) = conjg(splits(42)%Tcoeffs(1))
splits(42)%Vcoeffs(1) = 0.D0
splits(42)%Vcoeffs(2) = b1
splits(42)%Vcoeffs(3) = b2
splits(42)%Vcoeffs(4) = 1.D0-2.D0*(splits(42)%Vcoeffs(2)+splits(42)%Vcoeffs(3) )
splits(42)%Vcoeffs(5) = splits(42)%Vcoeffs(3)
splits(42)%Vcoeffs(6) = splits(42)%Vcoeffs(2)
end subroutine
end module splittings
Markdown is supported
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