3 * $Id: gfx.cal,v 1.1 2000/10/08 16:01:37 mdw Exp $
5 * Testbed for %$\gf{2}$% poltnomial arithmetic
7 * (c) 2000 Straylight/Edgeware
10 /*----- Licensing notice --------------------------------------------------*
12 * This file is part of Catacomb.
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.
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.
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,
30 /*----- Revision history --------------------------------------------------*
33 * Revision 1.1 2000/10/08 16:01:37 mdw
34 * Prototypes of various bits of code.
38 /*----- Object types ------------------------------------------------------*/
42 /*----- Static variables --------------------------------------------------*/
44 static obj gf example_gf_object;
46 /*----- Main code ---------------------------------------------------------*/
48 dummy = config("lib_debug", -1);
59 if (istype(x, example_gf_object))
65 define gf_add(x, y) = gf(xor(gfint(x), gfint(y)));
66 define gf_sub(x, y) = gf(xor(gfint(x), gfint(y)));
71 local a = gfint(x), b = gfint(y), z = 0, i, bits = highbit(a);
72 for (i = 0; i <= bits; i++) {
79 define gfx_div(rx, dx)
81 local r = gfint(rx), d = gfint(dx), i;
82 local q = 0, dbits = highbit(d), rbits = highbit(r);
83 for (i = rbits - dbits; i >= 0; i--) {
84 if (bit(r, i + dbits)) {
94 local l = gfx_div(x, y);
100 local l = gfx_div(x, y);
104 /*----- That's all, folks -------------------------------------------------*/