dbl_mp(ax,&mpx,p);  dbl_mp(ay,&mpy,p);  dvd(&mpy,&mpx,&mpz,p); dbl_mp(ue.d,&mpt1,p);   mul(&mpz,&mpt1,&mperr,p); sub(&mpz,&mperr,&mpz2,p);  __mp_dbl(&mpz2,&z,p);
__dbl_mp(ax,&mpx,p);  __dbl_mp(ay,&mpy,p);  __dvd(&mpy,&mpx,&mpz,p); __dbl_mp(ue.d,&mpt1,p);   __mul(&mpz,&mpt1,&mperr,p); __sub(&mpz,&mperr,&mpz2,p);  __mp_dbl(&mpz2,&z,p);
static double  normalized(double ax,double ay,double y, double z)    { int p;      mp_no mpx,mpy,mpz,mperr,mpz2,mpt1;  p=6;  dbl_mp(ax,&mpx,p);  dbl_mp(ay,&mpy,p);  dvd(&mpy,&mpx,&mpz,p);  dbl_mp(ue.d,&mpt1,p);   mul(&mpz,&mpt1,&mperr,p);  sub(&mpz,&mperr,&mpz2,p);  __mp_dbl(&mpz2,&z,p);  return signArctan2(y,z);}