for norm divide by sqrt(hd)

git-svn-id: http://svn.cactuscode.org/arrangements/CactusNumerical/Norms/trunk@5 f4913095-0e4f-0410-abea-a123d184f1f3

1 files changed, 12 insertions, 2 deletions

diff --git a/src/Compute_Norm.c b/src/Compute_Norm.c
index 083f498..482e57f 100644
--- a/src/Compute_Norm.c
+++ b/src/Compute_Norm.c
@@ -1,3 +1,6 @@
+
+
+
 /* XXX cleanup the include files ! */
 
 #include <assert.h>
@@ -175,8 +178,15 @@ void Norms_Compute_Norms (CCTK_ARGUMENTS)
   *norm=sqrt(*norm);
 
   /* multiply with h^d - see discussion above Eq. (12) in
-   * gr-qc/XXX */
-  hd=dx*dy*dz;
+   * gr-qc/0503056 [abs, ps, pdf, other] :
+   *
+   *  Title: Numerical stability for finite difference 
+   *         approximations of Einstein's equations
+   *  Authors: Gioel Calabrese, Ian Hinder, Sascha Husa
+   *  Comments: 18 pages, 5 figures
+   *
+   */
+  hd=sqrt(dx*dy*dz);
   *norm=*norm*hd;
 
   CCTK_VInfo(CCTK_THORNSTRING,"norm= %16.8g",*norm);