pone
u = pone(y); v = qone(y);
u = pone(x); v = qone(x);
static __inline double pone(double), qone(double);