diff options
Diffstat (limited to 'src/GeneralizedPolynomial-Uniform/template.c')
| -rw-r--r-- | src/GeneralizedPolynomial-Uniform/template.c | 24 |
1 files changed, 24 insertions, 0 deletions
diff --git a/src/GeneralizedPolynomial-Uniform/template.c b/src/GeneralizedPolynomial-Uniform/template.c index c12a308..454f733 100644 --- a/src/GeneralizedPolynomial-Uniform/template.c +++ b/src/GeneralizedPolynomial-Uniform/template.c @@ -863,6 +863,19 @@ int pt; } } +#ifdef PICKLE + #if (N_DIMS == 1) +printf("pt=%d interp coord %.16g\n", + pt, (double) interp_coords_fp[0]); + #endif + #if (N_DIMS == 3) +printf("pt=%d interp coords %.16g %.16g %.16g\n", + pt, (double) interp_coords_fp[0], + (double) interp_coords_fp[1], + (double) interp_coords_fp[2]); + #endif +#endif + /* * compute position of interpolation molecules with respect to @@ -934,6 +947,17 @@ int pt; const fp z = xyz_temp [Z_AXIS]; #endif +#ifdef PICKLE + #if (N_DIMS == 1) +printf("interp center_i = %d\n", center_i); +printf("interp grid-relative x = %.16g\n", (double) x); + #endif + #if (N_DIMS == 3) +printf("interp center_ijk = %d %d %d\n", center_i, center_j, center_k); +printf("interp grid-relative xyz = %.16g %.16g %.16g\n", + (double) x, (double) y, (double) z); + #endif +#endif /* * compute 1-d position of molecule center in input arrays |