mpfc_sub_mult()

void mpfc_sub_mult	(	gmp_CCmutable	result,
		gmp_CCmutable	a,
		gmp_CCmutable	b )

Definition at line 200 of file complex.c.

{
  // result->re -= a->re*b->re - a->im*b->im;
  // result->im -= a->re*b->im + a->im*b->re;
 
  mpfr_t tmp;
  mpfr_init2(tmp, mpfr_get_prec(a->re));
 
  mpfr_mul(tmp, a->re, b->re, MPFR_RNDN);
  mpfr_add(result->re, result->re, tmp, MPFR_RNDN);
  mpfr_mul(tmp, a->im, b->im, MPFR_RNDN);
  mpfr_sub(result->re, result->re, tmp, MPFR_RNDN);
 
  mpfr_mul(tmp, a->re, b->im, MPFR_RNDN);
  mpfr_add(result->im, result->im, tmp, MPFR_RNDN);
  mpfr_mul(tmp, a->im, b->re, MPFR_RNDN);
  mpfr_add(result->im, result->im, tmp, MPFR_RNDN);
 
  mpfr_clear(tmp);
}

References result().