Add cyclic group abstraction, with test code. Separate off exponentation
[u/mdw/catacomb] / pgen.h
diff --git a/pgen.h b/pgen.h
index 05d0a2c..1834f03 100644 (file)
--- a/pgen.h
+++ b/pgen.h
@@ -1,6 +1,6 @@
 /* -*-c-*-
  *
- * $Id: pgen.h,v 1.5 2000/02/12 18:21:03 mdw Exp $
+ * $Id: pgen.h,v 1.8 2004/04/01 12:50:09 mdw Exp $
  *
  * Prime generation glue
  *
 /*----- Revision history --------------------------------------------------* 
  *
  * $Log: pgen.h,v $
+ * Revision 1.8  2004/04/01 12:50:09  mdw
+ * Add cyclic group abstraction, with test code.  Separate off exponentation
+ * functions for better static linking.  Fix a buttload of bugs on the way.
+ * Generally ensure that negative exponents do inversion correctly.  Add
+ * table of standard prime-field subgroups.  (Binary field subgroups are
+ * currently unimplemented but easy to add if anyone ever finds a good one.)
+ *
+ * Revision 1.7  2000/08/18 19:16:12  mdw
+ * New event handler for showing in detail sub-prime generation.
+ *
+ * Revision 1.6  2000/06/17 11:52:12  mdw
+ * Add the GCD filter.
+ *
  * Revision 1.5  2000/02/12 18:21:03  mdw
  * Overhaul of key management (again).
  *
@@ -219,6 +232,24 @@ typedef struct pgen_safetestctx {
 
 extern int pgen_safetest(int /*rq*/, pgen_event */*ev*/, void */*p*/);
 
+/*----- Miscellaneous steppers and testers --------------------------------*/
+
+typedef struct pgen_gcdstepctx {
+  pfilt p, jp;                         /* Prime filter and step filter */
+  mp *q, *jq;                          /* %$p - 1$%, and a step value*/
+  mp *r;                               /* Other argument for GCD */
+  mp *g;                               /* GCD output (must be inited) */
+  mp *max;                             /* Maximum permissible GCD */
+} pgen_gcdstepctx;
+
+/* --- @pgen_gcdstep@ --- *
+ *
+ * Steps @p@ and @q@, until @p@ has no small factors, and
+ * %$\gcd(p, r) \le max$%.
+ */
+
+extern int pgen_gcdstep(int /*rq*/, pgen_event */*ev*/, void */*p*/);
+
 /*----- Standard event handlers -------------------------------------------*/
 
 /* --- @pgen_evspin@ --- *
@@ -235,6 +266,13 @@ extern int pgen_evspin(int /*rq*/, pgen_event */*ev*/, void */*p*/);
 
 extern int pgen_ev(int /*rq*/, pgen_event */*ev*/, void */*p*/);
 
+/* --- @pgen_subev@ --- *
+ *
+ * Subsidiary event handler, mainly for Lim-Lee searches and so on.
+ */
+
+extern int pgen_subev(int /*rq*/, pgen_event */*ev*/, void */*p*/);
+
 /*----- The main driver ---------------------------------------------------*/
 
 /* --- @pgen@ --- *
@@ -262,6 +300,16 @@ extern mp *pgen(const char */*name*/, mp */*d*/, mp */*m*/,
                unsigned /*steps*/, pgen_proc */*step*/, void */*sctx*/,
                unsigned /*tests*/, pgen_proc */*test*/, void */*tctx*/);
 
+/* --- @pgen_primep@ --- *
+ *
+ * Arguments:  @mp *p@ = a number to check
+ *             @grand *gr@ = a random number source
+ *
+ * Returns:    Nonzero if @p@ is really prime.
+ */
+
+extern int pgen_primep(mp */*p*/, grand */*gr*/);
+
 /*----- That's all, folks -------------------------------------------------*/
 
 #ifdef __cplusplus