diff options
| author | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2004-06-25 16:10:19 +0000 |
|---|---|---|
| committer | jthorn <jthorn@0f49ee68-0e4f-0410-9b9c-b2c123ded7ef> | 2004-06-25 16:10:19 +0000 |
| commit | 4223f30b01113a46ac12889355f7c5298a46c77e (patch) | |
| tree | e933bae232883039f8a8360e8f52312aa4569959 | |
| parent | 37d6f6cf743682363e5b3bc700c5a6a65b4c10d6 (diff) | |
change the error message printed if an interpolation point is
either outside the grid, or inside but too close to the grid boundary,
to also include the grid spacings. A sample of the new msg is:
WARNING[L1,P0] (AEILocalInterp):
CCTK_InterpLocalUniform():
interpolation point is either outside the grid,
or inside but too close to the grid boundary!
0-origin interpolation point number pt=691
interpolation point (x,y)=(2.26354,-0.591776)
grid x_min(delta_x)x_max = -0.1683(0.1122)3.30989
grid y_min(delta_y)y_max = -0.448799(0.224399)6.50758
git-svn-id: http://svn.aei.mpg.de/numrel/AEIThorns/AEILocalInterp/trunk@20 0f49ee68-0e4f-0410-9b9c-b2c123ded7ef
| -rw-r--r-- | src/template.c | 30 |
1 files changed, 18 insertions, 12 deletions
diff --git a/src/template.c b/src/template.c index 26be330..3b1670d 100644 --- a/src/template.c +++ b/src/template.c @@ -998,32 +998,38 @@ fflush(stdout); " 0-origin interpolation point number pt=%d\n" #if (N_DIMS == 1) " interpolation point x=%g\n" -" grid [x_min,x_max] = [%g,%g]\n" +" grid x_min(delta_x)x_max = %g(%g)%g\n" , pt, (double) interp_coords_fp[X_AXIS], - (double) grid_min_xyz[X_AXIS], (double) grid_max_xyz[X_AXIS]); + (double) grid_min_xyz[X_AXIS], (double) coord_delta [X_AXIS], + (double) grid_max_xyz[X_AXIS]); #elif (N_DIMS == 2) " interpolation point (x,y)=(%g,%g)\n" -" grid [x_min,x_max] = [%g,%g]\n" -" grid [y_min,y_max] = [%g,%g]\n" +" grid x_min(delta_x)x_max = %g(%g)%g\n" +" grid y_min(delta_y)y_max = %g(%g)%g\n" , pt, (double) interp_coords_fp[X_AXIS], (double) interp_coords_fp[Y_AXIS], - (double) grid_min_xyz[X_AXIS], (double) grid_max_xyz[X_AXIS], - (double) grid_min_xyz[Y_AXIS], (double) grid_max_xyz[Y_AXIS]); + (double) grid_min_xyz[X_AXIS], (double) coord_delta [X_AXIS], + (double) grid_max_xyz[X_AXIS], + (double) grid_min_xyz[Y_AXIS], (double) coord_delta [Y_AXIS], + (double) grid_max_xyz[Y_AXIS]); #elif (N_DIMS == 3) " interpolation point (x,y,z)=(%g,%g,%g)\n" -" grid [x_min,x_max] = [%g,%g]\n" -" grid [y_min,y_max] = [%g,%g]\n" -" grid [z_min,z_max] = [%g,%g]\n" +" grid x_min(delta_x)x_max = %g(%g)%g\n" +" grid y_min(delta_y)y_max = %g(%g)%g\n" +" grid z_min(delta_z)z_max = %g(%g)%g\n" , pt, (double) interp_coords_fp[X_AXIS], (double) interp_coords_fp[Y_AXIS], (double) interp_coords_fp[Z_AXIS], - (double) grid_min_xyz[X_AXIS], (double) grid_max_xyz[X_AXIS], - (double) grid_min_xyz[Y_AXIS], (double) grid_max_xyz[Y_AXIS], - (double) grid_min_xyz[Z_AXIS], (double) grid_max_xyz[Z_AXIS]); + (double) grid_min_xyz[X_AXIS], (double) coord_delta [X_AXIS], + (double) grid_max_xyz[X_AXIS], + (double) grid_min_xyz[Y_AXIS], (double) coord_delta [Y_AXIS], + (double) grid_max_xyz[Y_AXIS], + (double) grid_min_xyz[Z_AXIS], (double) coord_delta [Z_AXIS], + (double) grid_max_xyz[Z_AXIS]); #else #error "N_DIMS may not be > 3!" #endif |