double *nearest_interp_1d ( int nd, double xd[], double yd[], int ni, 
  double xi[] );