extern double vector_inprod_c(double *,double *);