version 1.4, 2001/10/09 01:36:11 |
version 1.5, 2018/03/29 01:32:52 |
|
|
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. |
* |
* |
* $OpenXM: OpenXM_contrib2/asir2000/engine/gf2n.c,v 1.3 2000/08/22 05:04:05 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/gf2n.c,v 1.4 2001/10/09 01:36:11 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 56 GEN_UP2 current_mod_gf2n; |
|
Line 56 GEN_UP2 current_mod_gf2n; |
|
|
|
void setmod_gf2n(P p) |
void setmod_gf2n(P p) |
{ |
{ |
if ( !current_mod_gf2n ) { |
if ( !current_mod_gf2n ) { |
current_mod_gf2n = (GEN_UP2)MALLOC(sizeof(struct oGEN_UP2)); |
current_mod_gf2n = (GEN_UP2)MALLOC(sizeof(struct oGEN_UP2)); |
current_mod_gf2n->id = UP2_SPARSE; /* use sparse rep. by default */ |
current_mod_gf2n->id = UP2_SPARSE; /* use sparse rep. by default */ |
} |
} |
ptoup2(p,¤t_mod_gf2n->dense); |
ptoup2(p,¤t_mod_gf2n->dense); |
ptoup2_sparse(p,¤t_mod_gf2n->sparse); |
ptoup2_sparse(p,¤t_mod_gf2n->sparse); |
} |
} |
|
|
void getmod_gf2n(UP2 *p) |
void getmod_gf2n(UP2 *p) |
{ |
{ |
if ( current_mod_gf2n ) |
if ( current_mod_gf2n ) |
*p = current_mod_gf2n->dense; |
*p = current_mod_gf2n->dense; |
else |
else |
*p = 0; |
*p = 0; |
} |
} |
|
|
void simpgf2n(GF2N n,GF2N *r) |
void simpgf2n(GF2N n,GF2N *r) |
{ |
{ |
UP2 rem; |
UP2 rem; |
|
|
if ( !n ) |
if ( !n ) |
*r = 0; |
*r = 0; |
else if ( NID(n) != N_GF2N ) |
else if ( NID(n) != N_GF2N ) |
*r = n; |
*r = n; |
else { |
else { |
gen_simpup2(n->body,current_mod_gf2n,&rem); |
gen_simpup2(n->body,current_mod_gf2n,&rem); |
MKGF2N(rem,*r); |
MKGF2N(rem,*r); |
} |
} |
} |
} |
|
|
void ptogf2n(Obj q,GF2N *l) |
void ptogf2n(Obj q,GF2N *l) |
{ |
{ |
UP2 q1; |
UP2 q1; |
|
|
if ( !q || (OID(q)==O_N && NID(q)==N_GF2N) ) { |
if ( !q || (OID(q)==O_N && NID(q)==N_GF2N) ) { |
*l = (GF2N)q; |
*l = (GF2N)q; |
} else if ( (OID(q)==O_N && NID(q)==N_Q) || OID(q)==O_P ) { |
} else if ( (OID(q)==O_N && NID(q)==N_Q) || OID(q)==O_P ) { |
ptoup2((P)q,&q1); |
ptoup2((P)q,&q1); |
MKGF2N(q1,*l); |
MKGF2N(q1,*l); |
} else |
} else |
error("ptogf2n : invalid argument"); |
error("ptogf2n : invalid argument"); |
} |
} |
|
|
void gf2ntop(GF2N q,P *l) |
void gf2ntop(GF2N q,P *l) |
{ |
{ |
if ( !q ) |
if ( !q ) |
*l = 0; |
*l = 0; |
else |
else |
up2top(q->body,l); |
up2top(q->body,l); |
} |
} |
|
|
void gf2ntovect(GF2N q,VECT *l) |
void gf2ntovect(GF2N q,VECT *l) |
{ |
{ |
if ( !q ) |
if ( !q ) |
*l = 0; |
*l = 0; |
else |
else |
up2tovect(q->body,l); |
up2tovect(q->body,l); |
} |
} |
|
|
#define NZGF2N(a) ((a)&&(OID(a)==O_N)&&(NID(a)==N_GF2N)) |
#define NZGF2N(a) ((a)&&(OID(a)==O_N)&&(NID(a)==N_GF2N)) |
|
|
void addgf2n(GF2N a,GF2N b,GF2N *c) |
void addgf2n(GF2N a,GF2N b,GF2N *c) |
{ |
{ |
UP2 t,t1; |
UP2 t,t1; |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
if ( !a ) |
if ( !a ) |
*c = b; |
*c = b; |
else if ( !b ) |
else if ( !b ) |
*c = a; |
*c = a; |
else { |
else { |
addup2(a->body,b->body,&t); |
addup2(a->body,b->body,&t); |
gen_simpup2(t,current_mod_gf2n,&t1); |
gen_simpup2(t,current_mod_gf2n,&t1); |
MKGF2N(t1,*c); |
MKGF2N(t1,*c); |
} |
} |
} |
} |
|
|
void subgf2n(GF2N a,GF2N b,GF2N *c) |
void subgf2n(GF2N a,GF2N b,GF2N *c) |
{ |
{ |
addgf2n(a,b,c); |
addgf2n(a,b,c); |
} |
} |
|
|
void mulgf2n(GF2N a,GF2N b,GF2N *c) |
void mulgf2n(GF2N a,GF2N b,GF2N *c) |
{ |
{ |
UP2 t; |
UP2 t; |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
if ( !a || !b ) |
if ( !a || !b ) |
*c = 0; |
*c = 0; |
else { |
else { |
mulup2(a->body,b->body,&t); |
mulup2(a->body,b->body,&t); |
#if 0 |
#if 0 |
gen_simpup2(t,current_mod_gf2n,&t1); |
gen_simpup2(t,current_mod_gf2n,&t1); |
MKGF2N(t1,*c); |
MKGF2N(t1,*c); |
#else |
#else |
gen_simpup2_destructive(t,current_mod_gf2n); |
gen_simpup2_destructive(t,current_mod_gf2n); |
if ( !t || !t->w ) |
if ( !t || !t->w ) |
*c = 0; |
*c = 0; |
else |
else |
MKGF2N(t,*c); |
MKGF2N(t,*c); |
#endif |
#endif |
} |
} |
} |
} |
|
|
void squaregf2n(GF2N a,GF2N *c) |
void squaregf2n(GF2N a,GF2N *c) |
{ |
{ |
UP2 t; |
UP2 t; |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)a,&z); a = z; |
ptogf2n((Obj)a,&z); a = z; |
if ( !a ) |
if ( !a ) |
*c = 0; |
*c = 0; |
else { |
else { |
squareup2(a->body,&t); |
squareup2(a->body,&t); |
#if 0 |
#if 0 |
gen_simpup2(t,current_mod_gf2n,&t1); |
gen_simpup2(t,current_mod_gf2n,&t1); |
MKGF2N(t1,*c); |
MKGF2N(t1,*c); |
#else |
#else |
gen_simpup2_destructive(t,current_mod_gf2n); |
gen_simpup2_destructive(t,current_mod_gf2n); |
if ( !t || !t->w ) |
if ( !t || !t->w ) |
*c = 0; |
*c = 0; |
else |
else |
MKGF2N(t,*c); |
MKGF2N(t,*c); |
#endif |
#endif |
} |
} |
} |
} |
|
|
void divgf2n(GF2N a,GF2N b,GF2N *c) |
void divgf2n(GF2N a,GF2N b,GF2N *c) |
{ |
{ |
UP2 t,i,s; |
UP2 t,i,s; |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
if ( !b ) |
if ( !b ) |
error("divgf2n: division by 0"); |
error("divgf2n: division by 0"); |
else if ( !a ) |
else if ( !a ) |
*c = 0; |
*c = 0; |
else { |
else { |
gen_invup2(b->body,current_mod_gf2n,&i); |
gen_invup2(b->body,current_mod_gf2n,&i); |
mulup2(a->body,i,&t); |
mulup2(a->body,i,&t); |
gen_simpup2(t,current_mod_gf2n,&s); |
gen_simpup2(t,current_mod_gf2n,&s); |
MKGF2N(s,*c); |
MKGF2N(s,*c); |
} |
} |
} |
} |
|
|
void invgf2n(GF2N b,GF2N *c) |
void invgf2n(GF2N b,GF2N *c) |
{ |
{ |
UP2 i; |
UP2 i; |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)b,&z); b = z; |
ptogf2n((Obj)b,&z); b = z; |
if ( !b ) |
if ( !b ) |
error("divgf2n: division by 0"); |
error("divgf2n: division by 0"); |
else { |
else { |
gen_invup2(b->body,current_mod_gf2n,&i); |
gen_invup2(b->body,current_mod_gf2n,&i); |
MKGF2N(i,*c); |
MKGF2N(i,*c); |
} |
} |
} |
} |
|
|
void chsgngf2n(GF2N a,GF2N *c) |
void chsgngf2n(GF2N a,GF2N *c) |
{ |
{ |
*c = a; |
*c = a; |
} |
} |
|
|
void pwrgf2n(GF2N a,Q b,GF2N *c) |
void pwrgf2n(GF2N a,Q b,GF2N *c) |
{ |
{ |
UP2 t; |
UP2 t; |
GF2N r; |
GF2N r; |
|
|
if ( !b ) { |
if ( !b ) { |
MKGF2N(ONEUP2,*c); |
MKGF2N(ONEUP2,*c); |
} else if ( !a ) |
} else if ( !a ) |
*c = 0; |
*c = 0; |
else { |
else { |
gen_pwrmodup2(a->body,b,current_mod_gf2n,&t); |
gen_pwrmodup2(a->body,b,current_mod_gf2n,&t); |
MKGF2N(t,r); |
MKGF2N(t,r); |
if ( SGN(b) < 0 ) |
if ( SGN(b) < 0 ) |
invgf2n(r,c); |
invgf2n(r,c); |
else |
else |
*c = r; |
*c = r; |
} |
} |
} |
} |
|
|
int cmpgf2n(GF2N a,GF2N b) |
int cmpgf2n(GF2N a,GF2N b) |
{ |
{ |
GF2N z; |
GF2N z; |
|
|
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
ptogf2n((Obj)a,&z); a = z; ptogf2n((Obj)b,&z); b = z; |
if ( !a ) |
if ( !a ) |
if ( !b ) |
if ( !b ) |
return 0; |
return 0; |
else |
else |
return -1; |
return -1; |
else if ( !b ) |
else if ( !b ) |
return 1; |
return 1; |
else |
else |
return compup2(a->body,b->body); |
return compup2(a->body,b->body); |
} |
} |
|
|
void randomgf2n(GF2N *r) |
void randomgf2n(GF2N *r) |
{ |
{ |
int i,w,d; |
int i,w,d; |
unsigned int *tb; |
unsigned int *tb; |
UP2 t; |
UP2 t; |
|
|
if ( !current_mod_gf2n ) |
if ( !current_mod_gf2n ) |
error("randomgf2n : current_mod_gf2n is not set"); |
error("randomgf2n : current_mod_gf2n is not set"); |
w = current_mod_gf2n->dense->w; |
w = current_mod_gf2n->dense->w; |
d = degup2(current_mod_gf2n->dense); |
d = degup2(current_mod_gf2n->dense); |
NEWUP2(t,w); |
NEWUP2(t,w); |
for ( i = 0, tb = t->b; i < w; i++ ) |
for ( i = 0, tb = t->b; i < w; i++ ) |
tb[i] = mt_genrand(); |
tb[i] = mt_genrand(); |
tb[w-1] &= (1<<(d%BSH))-1; |
tb[w-1] &= (1<<(d%BSH))-1; |
for ( i = w-1; i >= 0 && !tb[i]; i-- ); |
for ( i = w-1; i >= 0 && !tb[i]; i-- ); |
if ( i < 0 ) |
if ( i < 0 ) |
*r = 0; |
*r = 0; |
else { |
else { |
t->w = i+1; MKGF2N(t,*r); |
t->w = i+1; MKGF2N(t,*r); |
} |
} |
} |
} |