vt: fix %.*s in sendncars()

[plan9front.git] / sys / include / mp.h

diff --git a/sys/include/mp.h b/sys/include/mp.h
index 31e4ef3d83b80915beafdfcacda05bd2fca22451..0de6db39c8e67f90c4c02cb8a8dc580ec407ee02 100644 (file)
--- a/sys/include/mp.h
+++ b/sys/include/mp.h
@@ -8,7 +8,6 @@
  * mpdigit must be an atomic type.  mpdigit is defined
  * in the architecture specific u.h
  */
-
 typedef struct mpint mpint;
 
 struct mpint
@@ -25,6 +24,7 @@ enum
        MPstatic=       0x01,   /* static constant */
        MPnorm=         0x02,   /* normalization status */
        MPtimesafe=     0x04,   /* request time invariant computation */
+       MPfield=        0x08,   /* this mpint is a field modulus */
 
        Dbytes=         sizeof(mpdigit),        /* bytes per digit */
        Dbits=          Dbytes*8                /* bits per digit */
@@ -62,6 +62,8 @@ uvlong        mptouv(mpint*);                 /* unsigned vlong */
 mpint* uvtomp(uvlong, mpint*);
 vlong  mptov(mpint*);                  /* vlong */
 mpint* vtomp(vlong, mpint*);
+double mptod(mpint*);                  /* double */
+mpint* dtomp(double, mpint*);
 
 /* divide 2 digits by one */
 void   mpdigdiv(mpdigit *dividend, mpdigit divisor, mpdigit *quotient);
@@ -84,6 +86,7 @@ void  mpnot(mpint *b, mpint *res);
 void   mpxor(mpint *b1, mpint *b2, mpint *res);
 void   mptrunc(mpint *b, int n, mpint *res);
 void   mpxtend(mpint *b, int n, mpint *res);
+void   mpasr(mpint *b, int shift, mpint *res);
 
 /* modular arithmetic, time invariant when 0≤b1≤m-1 and 0≤b2≤m-1 */
 void   mpmodadd(mpint *b1, mpint *b2, mpint *m, mpint *sum);   /* sum = b1+b2 % m */
@@ -99,6 +102,9 @@ int  mpcmp(mpint *b1, mpint *b2);
 /* res = s != 0 ? b1 : b2 */
 void   mpsel(int s, mpint *b1, mpint *b2, mpint *res);
 
+/* return n! */
+mpint* mpfactorial(ulong n);
+
 /* extended gcd return d, x, and y, s.t. d = gcd(a,b) and ax+by = d */
 void   mpextendedgcd(mpint *a, mpint *b, mpint *d, mpint *x, mpint *y);
 
@@ -165,5 +171,18 @@ void       crtout(CRTpre*, CRTres*, mpint*);       /* convert residues to mpint */
 void   crtprefree(CRTpre*);
 void   crtresfree(CRTres*);
 
+/* fast field arithmetic */
+typedef struct Mfield  Mfield;
+
+struct Mfield
+{
+       mpint;
+       int     (*reduce)(Mfield*, mpint*, mpint*);
+};
+
+mpint *mpfield(mpint*);
+
+Mfield *gmfield(mpint*);
+Mfield *cnfield(mpint*);
 
 #pragma        varargck        type    "B"     mpint*