5 * (c) 2004 Straylight/Edgeware
8 /*----- Licensing notice --------------------------------------------------*
10 * This file is part of the Python interface to Catacomb.
12 * Catacomb/Python is free software; you can redistribute it and/or modify
13 * it under the terms of the GNU General Public License as published by
14 * the Free Software Foundation; either version 2 of the License, or
15 * (at your option) any later version.
17 * Catacomb/Python is distributed in the hope that it will be useful,
18 * but WITHOUT ANY WARRANTY; without even the implied warranty of
19 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
20 * GNU General Public License for more details.
22 * You should have received a copy of the GNU General Public License
23 * along with Catacomb/Python; if not, write to the Free Software Foundation,
24 * Inc., 59 Temple Place - Suite 330, Boston, MA 02111-1307, USA.
27 /*----- Header files ------------------------------------------------------*/
29 #include "catacomb-python.h"
31 /*----- Various utilities -------------------------------------------------*/
33 PyTypeObject
*field_pytype
;
34 PyTypeObject
*primefield_pytype
;
35 PyTypeObject
*niceprimefield_pytype
;
36 PyTypeObject
*binfield_pytype
;
37 PyTypeObject
*binpolyfield_pytype
;
38 PyTypeObject
*binnormfield_pytype
;
39 PyTypeObject
*fe_pytype
;
41 static PyObject
*fe_pynew(PyTypeObject
*ty
, PyObject
*arg
, PyObject
*kw
)
45 static const char *const kwlist
[] = { "x", 0 };
47 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "O:fe", KWLIST
, &x
))
49 if (FE_PYCHECK(x
) && FE_F(x
) == FIELD_F(ty
)) RETURN_OBJ(x
);
50 if ((z
= getmp(x
)) == 0) return (0);
51 z
= F_IN(FIELD_F(ty
), z
, z
);
52 return (fe_pywrap((PyObject
*)ty
, z
));
55 static PyObject
*field_dopywrap(PyTypeObject
*ty
, field
*f
)
57 field_pyobj
*fobj
= newtype(ty
, 0, f
->ops
->name
);
59 fobj
->ty
.ht_type
.tp_basicsize
= sizeof(fe_pyobj
);
60 fobj
->ty
.ht_type
.tp_base
= fe_pytype
;
62 fobj
->ty
.ht_type
.tp_flags
= (Py_TPFLAGS_DEFAULT
|
64 Py_TPFLAGS_CHECKTYPES
|
66 fobj
->ty
.ht_type
.tp_alloc
= PyType_GenericAlloc
;
67 fobj
->ty
.ht_type
.tp_free
= 0;
68 fobj
->ty
.ht_type
.tp_new
= fe_pynew
;
69 typeready(&fobj
->ty
.ht_type
);
70 return ((PyObject
*)fobj
);
73 PyObject
*field_pywrap(field
*f
)
77 if (STRCMP(F_NAME(f
), ==, "prime")) ty
= primefield_pytype
;
78 else if (STRCMP(F_NAME(f
), ==, "niceprime")) ty
= niceprimefield_pytype
;
79 else if (STRCMP(F_NAME(f
), ==, "binpoly")) ty
= binpolyfield_pytype
;
80 else if (STRCMP(F_NAME(f
), ==, "binnorm")) ty
= binnormfield_pytype
;
82 return (field_dopywrap(ty
, f
));
85 field
*field_copy(field
*f
)
87 if (STRCMP(F_NAME(f
), ==, "prime"))
88 f
= field_prime(f
->m
);
89 else if (STRCMP(F_NAME(f
), ==, "niceprime"))
90 f
= field_niceprime(f
->m
);
91 else if (STRCMP(F_NAME(f
), ==, "binpoly"))
92 f
= field_binpoly(f
->m
);
93 else if (STRCMP(F_NAME(f
), ==, "binnorm")) {
94 fctx_binnorm
*fc
= (fctx_binnorm
*)f
;
95 f
= field_binnorm(f
->m
, fc
->ntop
.r
[fc
->ntop
.n
- 1]);
101 PyObject
*fe_pywrap(PyObject
*fobj
, mp
*x
)
103 fe_pyobj
*z
= PyObject_New(fe_pyobj
, (PyTypeObject
*)fobj
);
104 z
->f
= FIELD_F(fobj
);
107 return ((PyObject
*)z
);
110 static mp
*implicitfe(field
*f
, PyObject
*o
)
114 if (FE_PYCHECK(o
) && FE_F(o
) == f
) return (MP_COPY(FE_X(o
)));
115 switch (f
->ops
->ty
) {
116 case FTY_PRIME
: x
= implicitmp(o
); break;
117 case FTY_BINARY
: x
= implicitgf(o
); break;
118 default: assert(!"huh?");
121 return (F_IN(f
, MP_NEW
, x
));
124 /*----- Field elements ----------------------------------------------------*/
126 static int febinop(PyObject
*x
, PyObject
*y
,
127 field
**f
, PyObject
**fobj
, mp
**xx
, mp
**yy
)
129 if (FE_PYCHECK(x
)) *fobj
= FE_FOBJ(x
);
130 else if (FE_PYCHECK(y
)) *fobj
= FE_FOBJ(y
);
133 if ((*xx
= implicitfe(*f
, x
)) == 0)
135 if ((*yy
= implicitfe(*f
, y
)) == 0) {
142 #define BINOP(name) \
143 static PyObject *fe_py##name(PyObject *x, PyObject *y) { \
147 if (febinop(x, y, &ff, &fobj, &xx, &yy)) RETURN_NOTIMPL; \
148 zz = ff->ops->name(ff, MP_NEW, xx, yy); \
149 MP_DROP(xx); MP_DROP(yy); \
150 return (fe_pywrap(fobj, zz)); \
157 static PyObject
*fe_pydiv(PyObject
*x
, PyObject
*y
)
163 if (febinop(x
, y
, &ff
, &fobj
, &xx
, &yy
)) RETURN_NOTIMPL
;
164 if (F_ZEROP(ff
, yy
)) ZDIVERR("division by zero");
165 yy
= F_INV(ff
, yy
, yy
);
166 z
= fe_pywrap(fobj
, F_MUL(ff
, MP_NEW
, xx
, yy
));
168 MP_DROP(xx
); MP_DROP(yy
);
172 static PyObject
*fe_pyexp(PyObject
*x
, PyObject
*y
, PyObject
*z
)
177 if (z
!= Py_None
|| !FE_PYCHECK(x
) || (yy
= tomp(y
)) == 0)
179 ff
= FE_F(x
); xx
= FE_X(x
); MP_COPY(xx
);
180 if (MP_NEGP(yy
) && F_ZEROP(ff
, xx
)) ZDIVERR("division by zero");
181 z
= fe_pywrap(FE_FOBJ(x
), field_exp(ff
, MP_NEW
, xx
, yy
));
183 MP_DROP(xx
); MP_DROP(yy
);
187 static PyObject
*fe_pyneg(PyObject
*x
)
189 return fe_pywrap(FE_FOBJ(x
), FE_F(x
)->ops
->neg(FE_F(x
), MP_NEW
, FE_X(x
)));
192 static PyObject
*fe_pyid(PyObject
*x
) { RETURN_OBJ(x
); }
194 static int fe_pynonzerop(PyObject
*x
) { return !F_ZEROP(FE_F(x
), FE_X(x
)); }
196 static PyObject
*fe_pyrichcompare(PyObject
*x
, PyObject
*y
, int op
)
204 if (febinop(x
, y
, &ff
, &fobj
, &xx
, &yy
)) RETURN_NOTIMPL
;
206 case Py_EQ
: b
= MP_EQ(xx
, yy
); break;
207 case Py_NE
: b
= !MP_EQ(xx
, yy
); break;
208 default: TYERR("field elements are unordered");
212 MP_DROP(xx
); MP_DROP(yy
);
216 static Py_hash_t
fe_pyhash(PyObject
*me
)
217 { return (mphash(FE_X(me
))); }
220 static int fe_pycoerce(PyObject
**x
, PyObject
**y
)
224 if (FE_PYCHECK(*y
)) {
225 if (FE_F(*x
) != FE_F(*y
) && !field_samep(FE_F(*x
), FE_F(*y
)))
226 TYERR("field mismatch");
227 Py_INCREF(*x
); Py_INCREF(*y
);
230 if ((z
= implicitfe(FE_F(*x
), *y
)) != 0) {
232 *y
= fe_pywrap(FE_FOBJ(*x
), z
);
242 static PyObject
*fe_pyint(PyObject
*x
)
246 mp
*xx
= F_OUT(FE_F(x
), MP_NEW
, FE_X(x
));
247 if (!mp_tolong_checked(xx
, &l
, 0)) rc
= PyInt_FromLong(l
);
248 else rc
= mp_topylong(xx
);
254 static PyObject
*fe_pylong(PyObject
*x
)
256 mp
*xx
= F_OUT(FE_F(x
), MP_NEW
, FE_X(x
));
257 PyObject
*rc
= mp_topylong(xx
);
263 #define BASEOP(name, radix, pre) \
264 static PyObject *fe_py##name(PyObject *x) { \
265 mp *xx = F_OUT(FE_F(x), MP_NEW, FE_X(x)); \
266 PyObject *rc = mp_topystring(xx, radix, 0, pre, 0); \
273 BASEOP(hex
, 16, "0x");
276 static void fe_pydealloc(PyObject
*me
)
278 Py_DECREF(FE_FOBJ(me
));
283 #define UNOP(name, check) \
284 static PyObject *femeth_##name(PyObject *me) { \
285 field *f = FE_F(me); \
287 if (!f->ops->name) TYERR(#name " not supported for this field"); \
289 x = f->ops->name(f, MP_NEW, x); \
290 if (!x) RETURN_NONE; \
291 return (fe_pywrap(FE_FOBJ(me), x)); \
295 UNOP(inv
, if (F_ZEROP(f
, x
)) ZDIVERR("division by zero"); )
305 static PyObject
*feget_value(PyObject
*me
, void *hunoz
)
307 mp
*x
= F_OUT(FE_F(me
), MP_NEW
, FE_X(me
));
308 if (F_TYPE(FE_F(me
)) == FTY_BINARY
)
309 return (gf_pywrap(x
));
311 return (mp_pywrap(x
));
314 static PyObject
*feget__value(PyObject
*me
, void *hunoz
)
318 if (F_TYPE(FE_F(me
)) == FTY_BINARY
)
319 return (gf_pywrap(x
));
321 return (mp_pywrap(x
));
324 static const PyMemberDef fe_pymembers
[] = {
325 #define MEMBERSTRUCT fe_pyobj
326 MEMRNM(field
, T_OBJECT
, PY23(ob_type
, ob_base
.ob_type
), READONLY
,
327 "X.field -> field containing X")
332 static const PyGetSetDef fe_pygetset
[] = {
333 #define GETSETNAME(op, name) fe##op##_##name
334 GET (value
, "X.value -> `natural' MP/GF representation of X")
335 GET (_value
, "X._value -> internal MP/GF representation of X")
340 static const PyMethodDef fe_pymethods
[] = {
341 #define METHNAME(func) femeth_##func
342 NAMETH(inv
, "X.inv() -> X^{-1}")
343 NAMETH(sqr
, "X.sqr() -> X^2")
344 NAMETH(sqrt
, "X.sqrt() -> sqrt(X)")
345 NAMETH(quadsolve
, "X.quadsolve() -> Y where Y^2 + Y = X (binary only)")
346 NAMETH(dbl
, "X.dbl() -> 2 * X (prime only)")
347 NAMETH(tpl
, "X.tpl() -> 3 * X (prime only)")
348 NAMETH(qdl
, "X.qdl() -> 4 * X (prime only)")
349 NAMETH(hlv
, "X.hlv() -> X/2 (prime only)")
354 static const PyNumberMethods fe_pynumber
= {
355 fe_pyadd
, /* @nb_add@ */
356 fe_pysub
, /* @nb_subtract@ */
357 fe_pymul
, /* @nb_multiply@ */
359 fe_pydiv
, /* @nb_divide@ */
361 0, /* @nb_remainder@ */
363 fe_pyexp
, /* @nb_power@ */
364 fe_pyneg
, /* @nb_negative@ */
365 fe_pyid
, /* @nb_positive@ */
366 0, /* @nb_absolute@ */
367 fe_pynonzerop
, /* @nb_nonzero@ */
375 fe_pycoerce
, /* @nb_coerce@ */
377 fe_pyint
, /* @nb_int@ */
378 PY23(fe_pylong
, 0), /* @nb_long@ */
379 0 /* meaningless */, /* @nb_float@ */
381 fe_pyoct
, /* @nb_oct@ */
382 fe_pyhex
, /* @nb_hex@ */
385 0, /* @nb_inplace_add@ */
386 0, /* @nb_inplace_subtract@ */
387 0, /* @nb_inplace_multiply@ */
389 0, /* @nb_inplace_divide@ */
391 0, /* @nb_inplace_remainder@ */
392 0, /* @nb_inplace_power@ */
393 0, /* @nb_inplace_lshift@ */
394 0, /* @nb_inplace_rshift@ */
395 0, /* @nb_inplace_and@ */
396 0, /* @nb_inplace_xor@ */
397 0, /* @nb_inplace_or@ */
399 0, /* @nb_floor_divide@ */
400 fe_pydiv
, /* @nb_true_divide@ */
401 0, /* @nb_inplace_floor_divide@ */
402 0, /* @nb_inplace_true_divide@ */
404 fe_pyint
, /* @nb_index@ */
407 static const PyTypeObject fe_pytype_skel
= {
408 PyVarObject_HEAD_INIT(0, 0) /* Header */
409 "FE", /* @tp_name@ */
410 sizeof(fe_pyobj
), /* @tp_basicsize@ */
411 0, /* @tp_itemsize@ */
413 fe_pydealloc
, /* @tp_dealloc@ */
415 0, /* @tp_getattr@ */
416 0, /* @tp_setattr@ */
417 0, /* @tp_compare@ */
419 PYNUMBER(fe
), /* @tp_as_number@ */
420 0, /* @tp_as_sequence@ */
421 0, /* @tp_as_mapping@ */
422 fe_pyhash
, /* @tp_hash@ */
424 fe_pyhex
, /* @tp_str@ */
425 0, /* @tp_getattro@ */
426 0, /* @tp_setattro@ */
427 0, /* @tp_as_buffer@ */
428 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
429 Py_TPFLAGS_CHECKTYPES
|
433 "Finite field elements, abstract base class.",
435 0, /* @tp_traverse@ */
437 fe_pyrichcompare
, /* @tp_richcompare@ */
438 0, /* @tp_weaklistoffset@ */
440 0, /* @tp_iternext@ */
441 PYMETHODS(fe
), /* @tp_methods@ */
442 PYMEMBERS(fe
), /* @tp_members@ */
443 PYGETSET(fe
), /* @tp_getset@ */
446 0, /* @tp_descr_get@ */
447 0, /* @tp_descr_set@ */
448 0, /* @tp_dictoffset@ */
450 PyType_GenericAlloc
, /* @tp_alloc@ */
451 abstract_pynew
, /* @tp_new@ */
456 /*----- Fields ------------------------------------------------------------*/
458 static PyObject
*field_pyrichcompare(PyObject
*x
, PyObject
*y
, int op
)
460 int b
= field_samep(FIELD_F(x
), FIELD_F(y
));
464 default: TYERR("can't order fields");
471 static PyObject
*fmeth_rand(PyObject
*me
, PyObject
*arg
, PyObject
*kw
)
473 static const char *const kwlist
[] = { "rng", 0 };
474 grand
*r
= &rand_global
;
476 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "|O&:rand", KWLIST
,
479 return (fe_pywrap(me
, F_RAND(FIELD_F(me
), MP_NEW
, r
)));
482 static PyObject
*fmeth__adopt(PyObject
*me
, PyObject
*arg
)
485 if (!PyArg_ParseTuple(arg
, "O&:_adopt", convmp
, &xx
)) return (0);
486 return (fe_pywrap(me
, xx
));
489 static PyObject
*fmeth_parse(PyObject
*me
, PyObject
*arg
)
496 if (!PyArg_ParseTuple(arg
, "s:parse", &p
)) goto end
;
498 if ((f
= field_parse(&qd
)) == 0) VALERR(qd
.e
);
499 rc
= Py_BuildValue("(Ns)", field_pywrap(f
), qd
.p
);
504 static void field_pydealloc(PyObject
*me
)
506 F_DESTROY(FIELD_F(me
));
507 PyType_Type
.tp_dealloc(me
);
510 static PyObject
*fget_zero(PyObject
*me
, void *hunoz
)
511 { return (fe_pywrap(me
, MP_COPY(FIELD_F(me
)->zero
))); }
513 static PyObject
*fget_one(PyObject
*me
, void *hunoz
)
514 { return (fe_pywrap(me
, MP_COPY(FIELD_F(me
)->one
))); }
516 static PyObject
*fget_q(PyObject
*me
, void *hunoz
)
517 { return (mp_pywrap(MP_COPY(FIELD_F(me
)->q
))); }
519 static PyObject
*fget_nbits(PyObject
*me
, void *hunoz
)
520 { return (PyInt_FromLong(FIELD_F(me
)->nbits
)); }
522 static PyObject
*fget_noctets(PyObject
*me
, void *hunoz
)
523 { return (PyInt_FromLong(FIELD_F(me
)->noctets
)); }
525 static PyObject
*fget_name(PyObject
*me
, void *hunoz
)
526 { return (TEXT_FROMSTR(F_NAME(FIELD_F(me
)))); }
528 static PyObject
*fget_type(PyObject
*me
, void *hunoz
)
529 { return (PyInt_FromLong(F_TYPE(FIELD_F(me
)))); }
531 static const PyGetSetDef field_pygetset
[] = {
532 #define GETSETNAME(op, name) f##op##_##name
533 GET (zero
, "F.zero -> field additive identity")
534 GET (one
, "F.one -> field multiplicative identity")
535 GET (q
, "F.q -> number of elements in field")
536 GET (nbits
, "F.nbits -> bits needed to represent element")
537 GET (noctets
, "F.noctets -> octetss needed to represent element")
538 GET (name
, "F.name -> name of this kind of field")
539 GET (type
, "F.type -> type code of this kind of field")
544 static const PyMethodDef field_pymethods
[] = {
545 #define METHNAME(name) fmeth_##name
546 METH (_adopt
, "F._adopt(X) -> FE")
547 KWMETH(rand
, "F.rand([rng = rand]) -> FE, uniformly distributed")
548 SMTH (parse
, "parse(STR) -> F, REST")
553 static const PyTypeObject field_pytype_skel
= {
554 PyVarObject_HEAD_INIT(0, 0) /* Header */
555 "Field", /* @tp_name@ */
556 sizeof(field_pyobj
), /* @tp_basicsize@ */
557 0, /* @tp_itemsize@ */
559 field_pydealloc
, /* @tp_dealloc@ */
561 0, /* @tp_getattr@ */
562 0, /* @tp_setattr@ */
563 0, /* @tp_compare@ */
565 0, /* @tp_as_number@ */
566 0, /* @tp_as_sequence@ */
567 0, /* @tp_as_mapping@ */
571 0, /* @tp_getattro@ */
572 0, /* @tp_setattro@ */
573 0, /* @tp_as_buffer@ */
574 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
578 "An abstract field. This is an abstract type.",
580 0, /* @tp_traverse@ */
582 field_pyrichcompare
, /* @tp_richcompare@ */
583 0, /* @tp_weaklistoffset@ */
585 0, /* @tp_iternext@ */
586 PYMETHODS(field
), /* @tp_methods@ */
587 0, /* @tp_members@ */
588 PYGETSET(field
), /* @tp_getset@ */
591 0, /* @tp_descr_get@ */
592 0, /* @tp_descr_set@ */
593 0, /* @tp_dictoffset@ */
595 PyType_GenericAlloc
, /* @tp_alloc@ */
596 abstract_pynew
, /* @tp_new@ */
601 /*----- Prime fields ------------------------------------------------------*/
603 static PyObject
*primefield_pynew(PyTypeObject
*ty
,
604 PyObject
*arg
, PyObject
*kw
)
608 static const char *const kwlist
[] = { "p", 0 };
610 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "O&:primefield", KWLIST
,
613 if ((f
= field_prime(xx
)) == 0)
614 VALERR("bad prime for primefield");
616 return (field_dopywrap(ty
, f
));
622 static PyObject
*pfget_p(PyObject
*me
, void *hunoz
)
623 { return (mp_pywrap(MP_COPY(FIELD_F(me
)->m
))); }
625 static const PyGetSetDef primefield_pygetset
[] = {
626 #define GETSETNAME(op, name) pf##op##_##name
627 GET (p
, "F.p -> prime field characteristic")
631 static const PyTypeObject primefield_pytype_skel
= {
632 PyVarObject_HEAD_INIT(0, 0) /* Header */
633 "PrimeField", /* @tp_name@ */
634 sizeof(field_pyobj
), /* @tp_basicsize@ */
635 0, /* @tp_itemsize@ */
637 field_pydealloc
, /* @tp_dealloc@ */
639 0, /* @tp_getattr@ */
640 0, /* @tp_setattr@ */
641 0, /* @tp_compare@ */
643 0, /* @tp_as_number@ */
644 0, /* @tp_as_sequence@ */
645 0, /* @tp_as_mapping@ */
649 0, /* @tp_getattro@ */
650 0, /* @tp_setattro@ */
651 0, /* @tp_as_buffer@ */
652 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
656 "PrimeField(P): prime fields.",
658 0, /* @tp_traverse@ */
660 0, /* @tp_richcompare@ */
661 0, /* @tp_weaklistoffset@ */
663 0, /* @tp_iternext@ */
664 0, /* @tp_methods@ */
665 0, /* @tp_members@ */
666 PYGETSET(primefield
), /* @tp_getset@ */
669 0, /* @tp_descr_get@ */
670 0, /* @tp_descr_set@ */
671 0, /* @tp_dictoffset@ */
673 PyType_GenericAlloc
, /* @tp_alloc@ */
674 primefield_pynew
, /* @tp_new@ */
679 static PyObject
*niceprimefield_pynew(PyTypeObject
*ty
,
680 PyObject
*arg
, PyObject
*kw
)
684 static const char *const kwlist
[] = { "p", 0 };
686 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "O&:niceprimefield",
687 KWLIST
, convmp
, &xx
))
689 if ((f
= field_niceprime(xx
)) == 0)
690 VALERR("bad prime for niceprimefield");
692 return (field_dopywrap(ty
, f
));
698 static const PyTypeObject niceprimefield_pytype_skel
= {
699 PyVarObject_HEAD_INIT(0, 0) /* Header */
700 "NicePrimeField", /* @tp_name@ */
701 sizeof(field_pyobj
), /* @tp_basicsize@ */
702 0, /* @tp_itemsize@ */
704 field_pydealloc
, /* @tp_dealloc@ */
706 0, /* @tp_getattr@ */
707 0, /* @tp_setattr@ */
708 0, /* @tp_compare@ */
710 0, /* @tp_as_number@ */
711 0, /* @tp_as_sequence@ */
712 0, /* @tp_as_mapping@ */
716 0, /* @tp_getattro@ */
717 0, /* @tp_setattro@ */
718 0, /* @tp_as_buffer@ */
719 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
723 "NicePrimeField(P): prime field using Solinas reduction.",
725 0, /* @tp_traverse@ */
727 0, /* @tp_richcompare@ */
728 0, /* @tp_weaklistoffset@ */
730 0, /* @tp_iternext@ */
731 0, /* @tp_methods@ */
732 0, /* @tp_members@ */
736 0, /* @tp_descr_get@ */
737 0, /* @tp_descr_set@ */
738 0, /* @tp_dictoffset@ */
740 PyType_GenericAlloc
, /* @tp_alloc@ */
741 niceprimefield_pynew
, /* @tp_new@ */
746 /*----- Binary fields -----------------------------------------------------*/
748 static PyObject
*bfget_m(PyObject
*me
, void *hunoz
)
749 { return (PyInt_FromLong(FIELD_F(me
)->nbits
)); }
751 static PyObject
*bfget_p(PyObject
*me
, void *hunoz
)
752 { return (gf_pywrap(MP_COPY(FIELD_F(me
)->m
))); }
754 static const PyGetSetDef binfield_pygetset
[] = {
755 #define GETSETNAME(op, name) bf##op##_##name
756 GET (m
, "F.m -> field polynomial degree")
757 GET (p
, "F.p -> field polynomial")
762 static const PyTypeObject binfield_pytype_skel
= {
763 PyVarObject_HEAD_INIT(0, 0) /* Header */
764 "BinField", /* @tp_name@ */
765 sizeof(field_pyobj
), /* @tp_basicsize@ */
766 0, /* @tp_itemsize@ */
768 field_pydealloc
, /* @tp_dealloc@ */
770 0, /* @tp_getattr@ */
771 0, /* @tp_setattr@ */
772 0, /* @tp_compare@ */
774 0, /* @tp_as_number@ */
775 0, /* @tp_as_sequence@ */
776 0, /* @tp_as_mapping@ */
780 0, /* @tp_getattro@ */
781 0, /* @tp_setattro@ */
782 0, /* @tp_as_buffer@ */
783 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
787 "Binary fields. Abstract class.",
789 0, /* @tp_traverse@ */
791 0, /* @tp_richcompare@ */
792 0, /* @tp_weaklistoffset@ */
794 0, /* @tp_iternext@ */
795 0, /* @tp_methods@ */
796 0, /* @tp_members@ */
797 PYGETSET(binfield
), /* @tp_getset@ */
800 0, /* @tp_descr_get@ */
801 0, /* @tp_descr_set@ */
802 0, /* @tp_dictoffset@ */
804 PyType_GenericAlloc
, /* @tp_alloc@ */
805 abstract_pynew
, /* @tp_new@ */
810 static PyObject
*binpolyfield_pynew(PyTypeObject
*ty
,
811 PyObject
*arg
, PyObject
*kw
)
815 static const char *const kwlist
[] = { "p", 0 };
817 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "O&:binpolyfield", KWLIST
,
820 if ((f
= field_binpoly(xx
)) == 0) VALERR("bad poly for binpolyfield");
822 return (field_dopywrap(ty
, f
));
828 static const PyTypeObject binpolyfield_pytype_skel
= {
829 PyVarObject_HEAD_INIT(0, 0) /* Header */
830 "BinPolyField", /* @tp_name@ */
831 sizeof(field_pyobj
), /* @tp_basicsize@ */
832 0, /* @tp_itemsize@ */
834 field_pydealloc
, /* @tp_dealloc@ */
836 0, /* @tp_getattr@ */
837 0, /* @tp_setattr@ */
838 0, /* @tp_compare@ */
840 0, /* @tp_as_number@ */
841 0, /* @tp_as_sequence@ */
842 0, /* @tp_as_mapping@ */
846 0, /* @tp_getattro@ */
847 0, /* @tp_setattro@ */
848 0, /* @tp_as_buffer@ */
849 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
853 "BinPolyField(P): binary fields with polynomial basis representation.",
855 0, /* @tp_traverse@ */
857 0, /* @tp_richcompare@ */
858 0, /* @tp_weaklistoffset@ */
860 0, /* @tp_iternext@ */
861 0, /* @tp_methods@ */
862 0, /* @tp_members@ */
866 0, /* @tp_descr_get@ */
867 0, /* @tp_descr_set@ */
868 0, /* @tp_dictoffset@ */
870 PyType_GenericAlloc
, /* @tp_alloc@ */
871 binpolyfield_pynew
, /* @tp_new@ */
876 static PyObject
*binnormfield_pynew(PyTypeObject
*ty
,
877 PyObject
*arg
, PyObject
*kw
)
881 static const char *const kwlist
[] = { "p", "beta", 0 };
883 if (!PyArg_ParseTupleAndKeywords(arg
, kw
, "O&O&:binnormfield",
884 KWLIST
, convgf
, &xx
, convgf
, &yy
))
886 if ((f
= field_binnorm(xx
, yy
)) == 0) VALERR("bad args for binnormfield");
887 MP_DROP(xx
); MP_DROP(yy
);
888 return (field_dopywrap(ty
, f
));
890 mp_drop(xx
); mp_drop(yy
);
894 static PyObject
*bnfget_beta(PyObject
*me
, void *hunoz
)
896 fctx_binnorm
*fc
= (fctx_binnorm
*)FIELD_F(me
);
897 return (gf_pywrap(MP_COPY(fc
->ntop
.r
[fc
->ntop
.n
- 1])));
900 static const PyGetSetDef binnormfield_pygetset
[] = {
901 #define GETSETNAME(op, name) bnf##op##_##name
902 GET (beta
, "F.beta -> conversion factor")
907 static const PyTypeObject binnormfield_pytype_skel
= {
908 PyVarObject_HEAD_INIT(0, 0) /* Header */
909 "BinNormField", /* @tp_name@ */
910 sizeof(field_pyobj
), /* @tp_basicsize@ */
911 0, /* @tp_itemsize@ */
913 field_pydealloc
, /* @tp_dealloc@ */
915 0, /* @tp_getattr@ */
916 0, /* @tp_setattr@ */
917 0, /* @tp_compare@ */
919 0, /* @tp_as_number@ */
920 0, /* @tp_as_sequence@ */
921 0, /* @tp_as_mapping@ */
925 0, /* @tp_getattro@ */
926 0, /* @tp_setattro@ */
927 0, /* @tp_as_buffer@ */
928 Py_TPFLAGS_DEFAULT
| /* @tp_flags@ */
932 "BinNormField(P, BETA): binary fields with normal basis representation.",
934 0, /* @tp_traverse@ */
936 0, /* @tp_richcompare@ */
937 0, /* @tp_weaklistoffset@ */
939 0, /* @tp_iternext@ */
940 0, /* @tp_methods@ */
941 0, /* @tp_members@ */
942 PYGETSET(binnormfield
), /* @tp_getset@ */
945 0, /* @tp_descr_get@ */
946 0, /* @tp_descr_set@ */
947 0, /* @tp_dictoffset@ */
949 PyType_GenericAlloc
, /* @tp_alloc@ */
950 binnormfield_pynew
, /* @tp_new@ */
955 /*----- Setup -------------------------------------------------------------*/
957 static const struct nameval consts
[] = {
958 CONST(FTY_PRIME
), CONST(FTY_BINARY
),
962 void field_pyinit(void)
965 INITTYPE(field
, type
);
966 INITTYPE(primefield
, field
);
967 INITTYPE(niceprimefield
, primefield
);
968 INITTYPE(binfield
, field
);
969 INITTYPE(binpolyfield
, binfield
);
970 INITTYPE(binnormfield
, binfield
);
973 void field_pyinsert(PyObject
*mod
)
975 INSERT("FE", fe_pytype
);
976 INSERT("Field", field_pytype
);
977 INSERT("PrimeField", primefield_pytype
);
978 INSERT("NicePrimeField", niceprimefield_pytype
);
979 INSERT("BinField", binfield_pytype
);
980 INSERT("BinPolyField", binpolyfield_pytype
);
981 INSERT("BinNormField", binnormfield_pytype
);
982 setconstants(mod
, consts
);
985 /*----- That's all, folks -------------------------------------------------*/