/* LibTomFloat, multiple-precision floating-point library * * LibTomFloat is a library that provides multiple-precision * floating-point artihmetic as well as trigonometric functionality. * * This library requires the public domain LibTomMath to be installed. *  * This library is free for all purposes without any express * gurantee it works * * Tom St Denis, tomstdenis@iahu.ca, http://float.libtomcrypt.org */#include <tomfloat.h> /* compute 1/x by (1/sqrt(x))^2 */int  mpf_inv(mp_float *a, mp_float *b){   int err, sign;    /* get sign of input */   sign = a->mantissa.sign;    /* force to positive */   a->mantissa.sign = MP_ZPOS;    /* compute 1/sqrt(a) */   if ((err = mpf_invsqrt(a, b)) != MP_OKAY) {      return err;   }    /* square 1/sqrt(a) to get 1/a */   err = mpf_sqr(b, b);    /* now restore the sign */   b->mantissa.sign = a->mantissa.sign = sign;    return err;}