chiark / gitweb /
Add some more vectors, and a whinge about how Skipjack test vectors are.
[catacomb] / factorial.c
1 /* -*-c-*-
2  *
3  * $Id: factorial.c,v 1.1 2000/07/09 21:30:49 mdw Exp $
4  *
5  * Example factorial computation
6  *
7  * (c) 2000 Straylight/Edgeware
8  */
9
10 /*----- Licensing notice --------------------------------------------------* 
11  *
12  * This file is part of Catacomb.
13  *
14  * Catacomb is free software; you can redistribute it and/or modify
15  * it under the terms of the GNU Library General Public License as
16  * published by the Free Software Foundation; either version 2 of the
17  * License, or (at your option) any later version.
18  * 
19  * Catacomb is distributed in the hope that it will be useful,
20  * but WITHOUT ANY WARRANTY; without even the implied warranty of
21  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
22  * GNU Library General Public License for more details.
23  * 
24  * You should have received a copy of the GNU Library General Public
25  * License along with Catacomb; if not, write to the Free
26  * Software Foundation, Inc., 59 Temple Place - Suite 330, Boston,
27  * MA 02111-1307, USA.
28  */
29
30 /*----- Revision history --------------------------------------------------* 
31  *
32  * $Log: factorial.c,v $
33  * Revision 1.1  2000/07/09 21:30:49  mdw
34  * Demo program to compute factorials.
35  *
36  */
37
38 /*----- Header files ------------------------------------------------------*/
39
40 #include <stdio.h>
41 #include <stdlib.h>
42
43 #include <mLib/quis.h>
44 #include <mLib/report.h>
45
46 #include "mpmul.h"
47
48 /*----- Main code ---------------------------------------------------------*/
49
50 int main(int argc, char *argv[])
51 {
52   unsigned long x;
53   char *p;
54   mp *f;
55
56   ego(argv[0]);
57
58   if (argc != 2) {
59     pquis(stderr, "Usage: $ integer\n");
60     exit(EXIT_FAILURE);
61   }
62   x = strtoul(argv[1], &p, 0);
63   if (*p)
64     die(EXIT_FAILURE, "bad integer `%s'", argv[1]);
65   f = mp_factorial(x);
66   mp_writefile(f, stdout, 10);
67   fputc('\n', stdout);
68   mp_drop(f);
69   return (0);
70 }
71
72 /*----- That's all, folks -------------------------------------------------*/