From 49a46b41c04c4d5847fe8f25f40599aab73e6d0f Mon Sep 17 00:00:00 2001 From: cott <> Date: Mon, 24 May 2004 08:50:00 +0000 Subject: Added a GF containing the relative error (makes more sense). darcs-hash:20040524085055-19929-7581401ba73496cb9150f92994ddcdb3ba01acc8.gz --- CarpetExtra/CarpetRegridTest/interface.ccl | 7 ++++--- CarpetExtra/CarpetRegridTest/src/TestGaussian.c | 4 ++-- 2 files changed, 6 insertions(+), 5 deletions(-) (limited to 'CarpetExtra') diff --git a/CarpetExtra/CarpetRegridTest/interface.ccl b/CarpetExtra/CarpetRegridTest/interface.ccl index 115b123f4..4e709662f 100644 --- a/CarpetExtra/CarpetRegridTest/interface.ccl +++ b/CarpetExtra/CarpetRegridTest/interface.ccl @@ -1,11 +1,12 @@ # Interface definition for thorn CarpetRegridTest -# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/CarpetRegridTest/interface.ccl,v 1.1 2004/05/23 15:56:14 cott Exp $ +# $Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/CarpetRegridTest/interface.ccl,v 1.2 2004/05/24 10:50:55 cott Exp $ implements: CarpetRegridTest inherits: grid -CCTK_REAL testscalar type=GF timelevels=3 tags='Prolongation="ENO" +#tags='Prolongation="ENO" +CCTK_REAL testscalar type=GF timelevels=3 { - phi, phi_error + phi, phi_error, phi_relerror } "The test scalar field" diff --git a/CarpetExtra/CarpetRegridTest/src/TestGaussian.c b/CarpetExtra/CarpetRegridTest/src/TestGaussian.c index fd68ff223..476b322e9 100644 --- a/CarpetExtra/CarpetRegridTest/src/TestGaussian.c +++ b/CarpetExtra/CarpetRegridTest/src/TestGaussian.c @@ -4,7 +4,7 @@ #include "cctk_Functions.h" -static const char *rcsid = "$Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/CarpetRegridTest/src/TestGaussian.c,v 1.1 2004/05/23 15:56:14 cott Exp $"; +static const char *rcsid = "$Header: /home/eschnett/C/carpet/Carpet/CarpetExtra/CarpetRegridTest/src/TestGaussian.c,v 1.2 2004/05/24 10:50:55 cott Exp $"; CCTK_FILEVERSION(CarpetExtra_CarpetRegridTest_TestGaussian_c); @@ -34,7 +34,7 @@ void CarpetRegrid_TestGaussian(CCTK_ARGUMENTS) R = sqrt(X*X + Y*Y + Z*Z); phi_error[index] = phi[index] - amplitude*exp( - pow( (R - radius) / sigma, 2.0 ) ); - + phi_relerror[index] = phi_error[index] / (amplitude*exp( - pow( (R - radius) / sigma, 2.0 ) ) ); } } } -- cgit v1.2.3