## File: bn_mp_sub.c

package info (click to toggle)
libtommath 0.42.0-1
• links: PTS
• area: main
• in suites: wheezy
• size: 4,764 kB
• sloc: ansic: 14,350; perl: 361; makefile: 255; asm: 30; sh: 1
 file content (59 lines) | stat: -rwxr-xr-x 1,674 bytes parent folder | download | duplicates (2)
 `1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859` ``````#include #ifdef BN_MP_SUB_C /* LibTomMath, multiple-precision integer library -- Tom St Denis * * LibTomMath is a library that provides multiple-precision * integer arithmetic as well as number theoretic functionality. * * The library was designed directly after the MPI library by * Michael Fromberger but has been written from scratch with * additional optimizations in place. * * The library is free for all purposes without any express * guarantee it works. * * Tom St Denis, tomstdenis@gmail.com, http://libtom.org */ /* high level subtraction (handles signs) */ int mp_sub (mp_int * a, mp_int * b, mp_int * c) { int sa, sb, res; sa = a->sign; sb = b->sign; if (sa != sb) { /* subtract a negative from a positive, OR */ /* subtract a positive from a negative. */ /* In either case, ADD their magnitudes, */ /* and use the sign of the first number. */ c->sign = sa; res = s_mp_add (a, b, c); } else { /* subtract a positive from a positive, OR */ /* subtract a negative from a negative. */ /* First, take the difference between their */ /* magnitudes, then... */ if (mp_cmp_mag (a, b) != MP_LT) { /* Copy the sign from the first */ c->sign = sa; /* The first has a larger or equal magnitude */ res = s_mp_sub (a, b, c); } else { /* The result has the *opposite* sign from */ /* the first number. */ c->sign = (sa == MP_ZPOS) ? MP_NEG : MP_ZPOS; /* The second has a larger magnitude */ res = s_mp_sub (b, a, c); } } return res; } #endif /* \$Source\$ */ /* \$Revision: 0.41 \$ */ /* \$Date: 2007-04-18 09:58:18 +0000 \$ */ ``````