diff options
author | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2004-10-11 20:16:00 +0000 |
---|---|---|
committer | diener <diener@f69c4107-0314-4c4f-9ad4-17e986b73f4a> | 2004-10-11 20:16:00 +0000 |
commit | 785d14d7d111a1120b9595385c38dde6267ade95 (patch) | |
tree | 87b09dee9a4dff8e3520a02fb188f185690a27c4 /src/Get_Coeff.F90 | |
parent | 820cf085c7fba79fa34163a535e5d7a37f1a2697 (diff) |
Added 4-3 Full restricted norm operators. Added parameter to choose the type
of differences (diagonal, full restricted). Added routine to return the
coefficient of the norm at the boundary. Added routine to return the
mask for the norm calculation.
git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@10 f69c4107-0314-4c4f-9ad4-17e986b73f4a
Diffstat (limited to 'src/Get_Coeff.F90')
-rw-r--r-- | src/Get_Coeff.F90 | 32 |
1 files changed, 32 insertions, 0 deletions
diff --git a/src/Get_Coeff.F90 b/src/Get_Coeff.F90 new file mode 100644 index 0000000..7b3fcf5 --- /dev/null +++ b/src/Get_Coeff.F90 @@ -0,0 +1,32 @@ +#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 + end select + end if +end function GetCoeff |