mpfc_mul()

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

Definition at line 67 of file complex.c.

{
  mpfr_t tmp;
  mpfr_init2(tmp, mpfr_get_prec(a->re));
 
  // result->re = a->re*b->re - a->im*b->im;
  mpfr_mul(tmp, a->re, b->re, MPFR_RNDN);
  mpfr_set(result->re, tmp, MPFR_RNDN);
  mpfr_mul(tmp, a->im, b->im, MPFR_RNDN);
  mpfr_sub(result->re, result->re, tmp, MPFR_RNDN);
 
  // result->im = a->re*b->im + a->im*b->re;
  mpfr_mul(tmp, a->re, b->im, MPFR_RNDN);
  mpfr_set(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().