extern void pgpt1_c(float *,float *,fint *);