divide() [1/2]

void M2::ARingCC::divide ( ElementType & res, const ElementType & a, const ElementType & b ) const

inline

Definition at line 327 of file aring-CC.hpp.

  {
    RealElementType p, denom;  // double
 
    ElementType result;
    if (fabs(b.re) >= fabs(b.im))
      {
        p = b.im / b.re;
        denom = b.re + p * b.im;
        result.re = (a.re + p * a.im) / denom;
        result.im = (a.im - p * a.re) / denom;
      }
    else
      {
        p = b.re / b.im;
        denom = b.im + p * b.re;
        result.re = (a.im + p * a.re) / denom;
        result.im = (p * a.im - a.re) / denom;
      }
    set(res, result);
  }

References cc_doubles_struct::im, p, cc_doubles_struct::re, result(), and set().