blob: 36c0209ef78b42d160cc35c96eb4d68d7f035ca8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
|
#include "cctk.h"
#include "cctk_Functions.h"
#include "cctk_Parameters.h"
CCTK_REAL function GetCoeff ()
implicit none
DECLARE_CCTK_FUNCTIONS
DECLARE_CCTK_PARAMETERS
CCTK_REAL, parameter :: zero = 0.0
integer, parameter :: wp = kind(zero)
if ( CCTK_EQUALS(norm_type,'Diagonal') ) then
select case (order)
case (2)
GetCoeff = 1.0_wp / 2.0_wp
case (4)
GetCoeff = 17.0_wp/48.0_wp
case (6)
GetCoeff = 13649.0_wp/43200.0_wp
case (8)
GetCoeff = 1498139.0_wp/5080320.0_wp
end select
else
select case (order)
case(4)
GetCoeff = 3.0_wp/11.0_wp
case(6)
GetCoeff = 30.0_wp/137.0_wp
end select
end if
end function GetCoeff
|