mpfc_invert()

void mpfc_invert	(	gmp_CCmutable	result,
		gmp_CCmutable	v )

Definition at line 86 of file complex.c.

{
  mpfr_t p, denom;
  mpfr_init2(p, mpfr_get_prec(v->re));
  mpfr_init2(denom, mpfr_get_prec(v->re));
 
  if (mpfr_cmpabs(v->re, v->im) >= 0)
    {
      // double p = v->im/v->re;
      // double denom = v->re + p * v->im;
      // result->re = 1.0/denom;
      // result->im = - p/denom;
 
      mpfr_div(p, v->im, v->re, MPFR_RNDN);
      mpfr_mul(denom, p, v->im, MPFR_RNDN);
      mpfr_add(denom, denom, v->re, MPFR_RNDN);
      mpfr_si_div(result->re, 1, denom, MPFR_RNDN);
      mpfr_div(result->im, p, denom, MPFR_RNDN);
      mpfr_neg(result->im, result->im, MPFR_RNDN);
    }
  else
    {
      // double p = v->re/v->im;
      // double denom = v->im + p * v->re;
      // result->re = p/denom;
      // result->im = -1.0/denom;
 
      mpfr_div(p, v->re, v->im, MPFR_RNDN);
      mpfr_mul(denom, p, v->re, MPFR_RNDN);
      mpfr_add(denom, denom, v->im, MPFR_RNDN);
      mpfr_si_div(result->im, 1, denom, MPFR_RNDN);
      mpfr_neg(result->im, result->im, MPFR_RNDN);
      mpfr_div(result->re, p, denom, MPFR_RNDN);
    }
 
  mpfr_clear(p);
  mpfr_clear(denom);
}

References p, and result().