Macaulay2 Engine
Loading...
Searching...
No Matches

◆ invert()

void M2::ARingCC::invert ( ElementType & res,
const ElementType & a ) const
inline

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

236 {
238 if (fabs(a.re) >= fabs(a.im))
239 {
240 double p = a.im / a.re;
241 double denom = a.re + p * a.im;
242 result.re = 1.0 / denom;
243 result.im = -p / denom;
244 }
245 else
246 {
247 double p = a.re / a.im;
248 double denom = a.im + p * a.re;
249 result.re = p / denom;
250 result.im = -1.0 / denom;
251 }
252 set(res, result);
253 }
void set(ElementType &result, const ElementType &a) const
Definition aring-CC.hpp:157
elem ElementType
Definition aring-CC.hpp:81
int p
VALGRIND_MAKE_MEM_DEFINED & result(result)
const mpreal fabs(const mpreal &x, mp_rnd_t r=mpreal::get_default_rnd())
Definition mpreal.h:2293

References p, result(), and set().

Referenced by power().