mult() [1/2]

void M2::ARingCCC::mult ( ElementType & res, const ElementType & a, const ElementType & b ) const

inline

Definition at line 348 of file aring-CCC.hpp.

  {
    mpfr_t tmp;
    Element result(*this);
    mpfr_init2(tmp, get_precision());
 
    // &result.re = &a.re*&b.re - &a.im*&b.im;
    mpfr_mul(tmp, &a.re, &b.re, MPFR_RNDN);
    mpfr_set(&result.value().re, tmp, MPFR_RNDN);
    mpfr_mul(tmp, &a.im, &b.im, MPFR_RNDN);
    mpfr_sub(&result.value().re, &result.value().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.value().im, tmp, MPFR_RNDN);
    mpfr_mul(tmp, &a.im, &b.re, MPFR_RNDN);
    mpfr_add(&result.value().im, &result.value().im, tmp, MPFR_RNDN);
 
    set(res, result);
    mpfr_clear(tmp);
  }

References get_precision(), cc_struct::im, cc_struct::re, result(), and set().