From 4223f30b01113a46ac12889355f7c5298a46c77e Mon Sep 17 00:00:00 2001 From: jthorn Date: Fri, 25 Jun 2004 16:10:19 +0000 Subject: 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 --- src/template.c | 30 ++++++++++++++++++------------ 1 file 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 -- cgit v1.2.3