diff options
Diffstat (limited to 'src/template.c')
-rw-r--r-- | src/template.c | 31 |
1 files changed, 14 insertions, 17 deletions
diff --git a/src/template.c b/src/template.c index 6c46e4e..26be330 100644 --- a/src/template.c +++ b/src/template.c @@ -996,40 +996,37 @@ fflush(stdout); " interpolation point is either outside the grid,\n" " or inside but too close to the grid boundary!\n" " 0-origin interpolation point number pt=%d\n" -" 0-origin coordinate axis=%d\n" #if (N_DIMS == 1) " interpolation point x=%g\n" " grid [x_min,x_max] = [%g,%g]\n" + , + pt, + (double) interp_coords_fp[X_AXIS], + (double) grid_min_xyz[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" + , + 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]); #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" -#else - #error "N_DIMS may not be > 3!" -#endif - , - pt, axis, - #if (N_DIMS == 1) - (double) interp_coords_fp[X_AXIS], - (double) grid_min_xyz[X_AXIS], (double) grid_max_xyz[X_AXIS]); - #elif (N_DIMS == 2) - (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]); - #elif (N_DIMS == 3) + , + 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]); - #else - #error "N_DIMS may not be > 3!" - #endif +#else + #error "N_DIMS may not be > 3!" +#endif } if (this_point_status < return_status) |