sys/src/libmp/port/mpsub.c

   1 #include "os.h"
   2 #include <mp.h>
   3 #include "dat.h"
   4
   5 // diff = abs(b1) - abs(b2), i.e., subtract the magnitudes
   6 void
   7 mpmagsub(mpint *b1, mpint *b2, mpint *diff)
   8 {
   9         int n, m, sign;
  10         mpint *t;
  11
  12         // get the sizes right
  13         if(mpmagcmp(b1, b2) < 0){
  14                 assert(((b1->flags | b2->flags | diff->flags) & MPtimesafe) == 0);
  15                 sign = -1;
  16                 t = b1;
  17                 b1 = b2;
  18                 b2 = t;
  19         } else {
  20                 diff->flags |= (b1->flags | b2->flags) & MPtimesafe;
  21                 sign = 1;
  22         }
  23         n = b1->top;
  24         m = b2->top;
  25         if(m == 0){
  26                 mpassign(b1, diff);
  27                 diff->sign = sign;
  28                 return;
  29         }
  30         mpbits(diff, n*Dbits);
  31
  32         mpvecsub(b1->p, n, b2->p, m, diff->p);
  33         diff->sign = sign;
  34         diff->top = n;
  35         mpnorm(diff);
  36 }
  37
  38 // diff = b1 - b2
  39 void
  40 mpsub(mpint *b1, mpint *b2, mpint *diff)
  41 {
  42         int sign;
  43
  44         if(b1->sign != b2->sign){
  45                 assert(((b1->flags | b2->flags | diff->flags) & MPtimesafe) == 0);
  46                 sign = b1->sign;
  47                 mpmagadd(b1, b2, diff);
  48                 diff->sign = sign;
  49                 return;
  50         }
  51
  52         sign = b1->sign;
  53         mpmagsub(b1, b2, diff);
  54         if(diff->top != 0)
  55                 diff->sign *= sign;
  56 }