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