double *pwl_basis_1d ( int nd, double xd[], int k, int ni, double xi[] );
double *pwl_interp_1d ( int nd, double xd[], double yd[], int ni, double xi[] );