38 double x = ncoords >=1 ? coords[0] : 0.0;
39 double y = ncoords >=2 ? coords[1] : 0.0;
40 double z = ncoords >=3 ? coords[2] : 0.0;
41 double mag = sqrt(x*x+y*y+z*z);
45 if (ncoords >= 1) results[0] = xhat;
46 if (ncoords >= 2) results[1] = yhat;
47 if (ncoords >= 3) results[2] = zhat;
48 return MACSIO_EVAL_OK;
macsio_eval_state_t normal_vector(int ncoords, double const *coords, double *results)