Realized that in the CVS version the optimized restricted full operators

were using the non-optimized values for the norm at the boundary (the
element of the norm used in the penalty term). The results in the SBP paper
were all computed with versions of the code using the right values. It
shouldn't make much difference anyway.


git-svn-id: https://svn.cct.lsu.edu/repos/numrel/LSUThorns/SummationByParts/trunk@42 f69c4107-0314-4c4f-9ad4-17e986b73f4a

1 files changed, 15 insertions, 6 deletions

diff --git a/src/Get_Coeff.F90 b/src/Get_Coeff.F90
index 36c0209..3c7a09f 100644
--- a/src/Get_Coeff.F90
+++ b/src/Get_Coeff.F90
@@ -24,11 +24,20 @@ CCTK_REAL function GetCoeff ()
       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
+    if ( CCTK_EQUALS(operator_type,'Minimal Bandwidth') ) then
+      select case (order)
+      case(4)
+        GetCoeff = 3.0_wp/11.0_wp
+      case(6)
+        GetCoeff = 30.0_wp/137.0_wp
+      end select
+    else
+      select case (order)
+      case(4)
+        GetCoeff = 0.2388575707774486064323157210922003533466_wp
+      case(6)
+        GetCoeff = 0.2028105550720356346665604029847379994496_wp
+      end select
+    end if
   end if
 end function GetCoeff