X-Git-Url: https://git.distorted.org.uk/u/mdw/catacomb/blobdiff_plain/62f34e465125f07961ae8cc8bff167f9cdb6f80a..3563e36580c7dad68cd6d3f7eb82eef570fc0c76:/genprimes.c

diff --git a/genprimes.c b/genprimes.c
index 2961621..2ee4115 100644
--- a/genprimes.c
+++ b/genprimes.c
@@ -1,6 +1,6 @@
 /* -*-c-*-
  *
- * $Id: genprimes.c,v 1.4 2001/03/04 13:08:10 mdw Exp $
+ * $Id: genprimes.c,v 1.6 2004/04/02 01:03:49 mdw Exp $
  *
  * Generate prime number table
  *
@@ -30,6 +30,16 @@
 /*----- Revision history --------------------------------------------------* 
  *
  * $Log: genprimes.c,v $
+ * Revision 1.6  2004/04/02 01:03:49  mdw
+ * Miscellaneous constification.
+ *
+ * Revision 1.5  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.4  2001/03/04 13:08:10  mdw
  * Use @DA_LAST@ to determine @MAXPRIME@, now that it exists.
  *
@@ -71,13 +81,14 @@ int main(int argc, char *argv[])
   char *header = "primetab.h";
   char *source = "primetab.c";
   char *name = "primetab";
+  char *sym = 0;
   intv p = DA_INIT;
   int i;
 
   ego(argv[0]);
 
   for (;;) {
-    int i = getopt(argc, argv, "h:c:i:n:m:t:");
+    int i = getopt(argc, argv, "h:c:i:n:m:t:s:");
     if (i < 0)
       break;
     switch (i) {
@@ -101,6 +112,9 @@ int main(int argc, char *argv[])
       case 't':
 	type = optarg;
 	break;
+      case 's':
+	sym = optarg;
+	break;
       default:
 	pquis(stderr, "Usage: $ [-n nprimes] [-m maxprime] [-t type]\n");
 	exit(EXIT_FAILURE);
@@ -131,15 +145,18 @@ int main(int argc, char *argv[])
     char *q;
     if (!fp)
       die(EXIT_FAILURE, "couldn't write `%s': %s", header, strerror(errno));
-    for (q = header; *q; q++) {
-      int ch = (unsigned char)*q;
-      if (isalnum(ch))
-	ch = toupper(ch);
-      else
-	ch = '_';
-      DPUTC(&d, ch);
+    if (!sym) {
+      for (q = header; *q; q++) {
+	int ch = (unsigned char)*q;
+	if (isalnum(ch))
+	  ch = toupper(ch);
+	else
+	  ch = '_';
+	DPUTC(&d, ch);
+      }
+      DPUTZ(&d);
+      sym = d.buf;
     }
-    DPUTZ(&d);
     fprintf(fp, "\
 /* -*-c-*-\n\
  *\n\
@@ -153,11 +170,11 @@ int main(int argc, char *argv[])
 #define MAXPRIME %uu\n\
 \n\
 typedef %s smallprime;\n\
-extern smallprime %s[];\n\
+extern const smallprime %s[];\n\
 \n\
 #endif\n\
 ",
-	    d.buf, d.buf,
+	    sym, sym,
 	    (unsigned long)DA_LEN(&p),
 	    DA_LAST(&p),
 	    type, name);
@@ -181,7 +198,7 @@ extern smallprime %s[];\n\
 \n\
 #include \"%s\"\n\
 \n\
-%s %s[] = {",
+const %s %s[] = {",
 	    header, type, name);
     for (i = 0; i < DA_LEN(&p); i++) {
       if (i % 8 == 0)