![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to gf2n.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / engine|
Return to gf2n.c CVS log | Up to [local] / OpenXM_contrib2 / asir2000 / engine |
1.1 noro 1: /* $OpenXM: OpenXM/src/asir99/engine/gf2n.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
2: #include "ca.h"
3: #include "base.h"
4:
5: extern int lm_lazy;
6:
7: GEN_UP2 current_mod_gf2n;
8:
9: void setmod_gf2n(p)
10: P p;
11: {
12: if ( !current_mod_gf2n ) {
13: current_mod_gf2n = (GEN_UP2)MALLOC(sizeof(struct oGEN_UP2));
14: current_mod_gf2n->id = UP2_SPARSE; /* use sparse rep. by default */
15: }
16: ptoup2(p,¤t_mod_gf2n->dense);
17: ptoup2_sparse(p,¤t_mod_gf2n->sparse);
18: }
19:
20: void getmod_gf2n(p)
21: UP2 *p;
22: {
23: if ( current_mod_gf2n )
24: *p = current_mod_gf2n->dense;
25: else
26: *p = 0;
27: }
28:
29: void simpgf2n(n,r)
30: GF2N n;
31: GF2N *r;
32: {
33: UP2 rem;
34:
35: if ( !n )
36: *r = 0;
37: else if ( NID(n) != N_GF2N )
38: *r = n;
39: else {
40: gen_simpup2(n->body,current_mod_gf2n,&rem);
41: MKGF2N(rem,*r);
42: }
43: }
44:
45: void ptogf2n(q,l)
46: Obj q;
47: GF2N *l;
48: {
49: UP2 q1;
50:
51: if ( !q || (OID(q)==O_N && NID(q)==N_GF2N) ) {
52: *l = (GF2N)q;
53: } else if ( (OID(q)==O_N && NID(q)==N_Q) || OID(q)==O_P ) {
54: ptoup2((P)q,&q1);
55: MKGF2N(q1,*l);
56: } else
57: error("ptogf2n : invalid argument");
58: }
59:
60: void gf2ntop(q,l)
61: GF2N q;
62: P *l;
63: {
64: if ( !q )
65: *l = 0;
66: else
67: up2top(q->body,l);
68: }
69:
70: void gf2ntovect(q,l)
71: GF2N q;
72: VECT *l;
73: {
74: if ( !q )
75: *l = 0;
76: else
77: up2tovect(q->body,l);
78: }
79:
80: #define NZGF2N(a) ((a)&&(OID(a)==O_N)&&(NID(a)==N_GF2N))
81:
82: void addgf2n(a,b,c)
83: GF2N a,b;
84: GF2N *c;
85: {
86: UP2 t,t1;
87: GF2N z;
88:
89: ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z;
90: if ( !a )
91: *c = b;
92: else if ( !b )
93: *c = a;
94: else {
95: addup2(a->body,b->body,&t);
96: gen_simpup2(t,current_mod_gf2n,&t1);
97: MKGF2N(t1,*c);
98: }
99: }
100:
101: void subgf2n(a,b,c)
102: GF2N a,b;
103: GF2N *c;
104: {
105: addgf2n(a,b,c);
106: }
107:
108: void mulgf2n(a,b,c)
109: GF2N a,b;
110: GF2N *c;
111: {
112: UP2 t;
113: GF2N z;
114:
115: ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z;
116: if ( !a || !b )
117: *c = 0;
118: else {
119: mulup2(a->body,b->body,&t);
120: #if 0
121: gen_simpup2(t,current_mod_gf2n,&t1);
122: MKGF2N(t1,*c);
123: #else
124: gen_simpup2_destructive(t,current_mod_gf2n);
125: if ( !t || !t->w )
126: *c = 0;
127: else
128: MKGF2N(t,*c);
129: #endif
130: }
131: }
132:
133: void squaregf2n(a,c)
134: GF2N a;
135: GF2N *c;
136: {
137: UP2 t;
138: GF2N z;
139:
140: ptogf2n((Obj)a,&z); a = z;
141: if ( !a )
142: *c = 0;
143: else {
144: squareup2(a->body,&t);
145: #if 0
146: gen_simpup2(t,current_mod_gf2n,&t1);
147: MKGF2N(t1,*c);
148: #else
149: gen_simpup2_destructive(t,current_mod_gf2n);
150: if ( !t || !t->w )
151: *c = 0;
152: else
153: MKGF2N(t,*c);
154: #endif
155: }
156: }
157:
158: void divgf2n(a,b,c)
159: GF2N a,b;
160: GF2N *c;
161: {
162: UP2 t,i,s;
163: GF2N z;
164:
165: ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z;
166: if ( !b )
167: error("divgf2n: division by 0");
168: else if ( !a )
169: *c = 0;
170: else {
171: gen_invup2(b->body,current_mod_gf2n,&i);
172: mulup2(a->body,i,&t);
173: gen_simpup2(t,current_mod_gf2n,&s);
174: MKGF2N(s,*c);
175: }
176: }
177:
178: void invgf2n(b,c)
179: GF2N b;
180: GF2N *c;
181: {
182: UP2 i;
183: GF2N z;
184:
185: ptogf2n((Obj)b,&z); b = z;
186: if ( !b )
187: error("divgf2n: division by 0");
188: else {
189: gen_invup2(b->body,current_mod_gf2n,&i);
190: MKGF2N(i,*c);
191: }
192: }
193:
194: void chsgngf2n(a,c)
195: GF2N a,*c;
196: {
197: *c = a;
198: }
199:
200: void pwrgf2n(a,b,c)
201: GF2N a;
202: Q b;
203: GF2N *c;
204: {
205: UP2 t;
206: GF2N r;
207:
208: if ( !b ) {
209: MKGF2N(ONEUP2,*c);
210: } else if ( !a )
211: *c = 0;
212: else {
213: gen_pwrmodup2(a->body,b,current_mod_gf2n,&t);
214: MKGF2N(t,r);
215: if ( SGN(b) < 0 )
216: invgf2n(r,c);
217: else
218: *c = r;
219: }
220: }
221:
222: int cmpgf2n(a,b)
223: GF2N a,b;
224: {
225: GF2N z;
226:
227: ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z;
228: if ( !a )
229: if ( !b )
230: return 0;
231: else
232: return -1;
233: else if ( !b )
234: return 1;
235: else
236: return compup2(a->body,b->body);
237: }
238:
239: void randomgf2n(r)
240: GF2N *r;
241: {
242: int i,w,d;
243: unsigned int *tb;
244: UP2 t;
245:
246: if ( !current_mod_gf2n )
247: error("randomgf2n : current_mod_gf2n is not set");
248: w = current_mod_gf2n->dense->w;
249: d = degup2(current_mod_gf2n->dense);
250: NEWUP2(t,w);
251: for ( i = 0, tb = t->b; i < w; i++ )
252: tb[i] = mt_genrand();
253: tb[w-1] &= (1<<(d%BSH))-1;
254: for ( i = w-1; i >= 0 && !tb[i]; i-- );
255: if ( i < 0 )
256: *r = 0;
257: else {
258: t->w = i+1; MKGF2N(t,*r);
259: }
260: }