| version 1.3, 2000/08/22 05:04:07 |
version 1.4, 2001/10/09 01:36:14 |
|
|
| * 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/up_gf2n.c,v 1.2 2000/08/21 08:31:28 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/up_gf2n.c,v 1.3 2000/08/22 05:04:07 noro Exp $ |
| */ |
*/ |
| #include "ca.h" |
#include "ca.h" |
| #include <math.h> |
#include <math.h> |
| Line 54 extern int debug_up; |
|
| Line 54 extern int debug_up; |
|
| extern int up_lazy; |
extern int up_lazy; |
| extern GEN_UP2 current_mod_gf2n; |
extern GEN_UP2 current_mod_gf2n; |
| |
|
| void squarep_gf2n(vl,n1,nr) |
void squarep_gf2n(VL vl,P n1,P *nr) |
| VL vl; |
|
| P n1; |
|
| P *nr; |
|
| { |
{ |
| UP b1,br; |
UP b1,br; |
| |
|
|
|
| } |
} |
| } |
} |
| |
|
| void squareup_gf2n(n1,nr) |
void squareup_gf2n(UP n1,UP *nr) |
| UP n1; |
|
| UP *nr; |
|
| { |
{ |
| UP r; |
UP r; |
| GF2N *c1,*c; |
GF2N *c1,*c; |
|
|
| |
|
| /* x^(2^n) mod f */ |
/* x^(2^n) mod f */ |
| |
|
| void powermodup_gf2n(f,xp) |
void powermodup_gf2n(UP f,UP *xp) |
| UP f; |
|
| UP *xp; |
|
| { |
{ |
| UP x,t,invf; |
UP x,t,invf; |
| int k,n; |
int k,n; |
| GF2N lm; |
GF2N lm; |
| struct oEGT eg_sq,eg_rem,eg_mul,eg_inv,eg0,eg1,eg2; |
|
| |
|
| n = degup2(current_mod_gf2n->dense); |
n = degup2(current_mod_gf2n->dense); |
| MKGF2N(ONEUP2,lm); |
MKGF2N(ONEUP2,lm); |
|
|
| |
|
| /* g^d mod f */ |
/* g^d mod f */ |
| |
|
| void generic_powermodup_gf2n(g,f,d,xp) |
void generic_powermodup_gf2n(UP g,UP f,Q d,UP *xp) |
| UP g,f; |
|
| Q d; |
|
| UP *xp; |
|
| { |
{ |
| N e; |
N e; |
| UP x,y,t,invf,s; |
UP x,y,t,invf,s; |
|
|
| |
|
| /* g+g^2+...+g^(2^(nd-1)) mod f; where e = deg(mod) */ |
/* g+g^2+...+g^(2^(nd-1)) mod f; where e = deg(mod) */ |
| |
|
| void tracemodup_gf2n(g,f,d,xp) |
void tracemodup_gf2n(UP g,UP f,Q d,UP *xp) |
| UP g,f; |
|
| Q d; |
|
| UP *xp; |
|
| { |
{ |
| UP x,t,s,u,invf; |
UP x,t,s,u,invf; |
| int en,i; |
int en,i; |
|
|
| *xp = s; |
*xp = s; |
| } |
} |
| |
|
| void tracemodup_gf2n_slow(g,f,d,xp) |
void tracemodup_gf2n_slow(UP g,UP f,Q d,UP *xp) |
| UP g,f; |
|
| Q d; |
|
| UP *xp; |
|
| { |
{ |
| UP x,t,s,u; |
UP x,t,s,u; |
| int en,i; |
int en,i; |
|
|
| *xp = s; |
*xp = s; |
| } |
} |
| |
|
| static struct oEGT eg_trace_tab,eg_trace_mul; |
void tracemodup_gf2n_tab(UP g,UP f,Q d,UP *xp) |
| |
|
| void tracemodup_gf2n_tab(g,f,d,xp) |
|
| UP g,f; |
|
| Q d; |
|
| UP *xp; |
|
| { |
{ |
| UP x0,x2,t,s,u; |
UP x0,x2,t,s,u; |
| int en,i; |
int en,i; |
| UP *tab; |
UP *tab; |
| GF2N one; |
GF2N one; |
| struct oEGT eg1,eg2; |
|
| |
|
| en = QTOS(d)*degup2(current_mod_gf2n->dense); |
en = QTOS(d)*degup2(current_mod_gf2n->dense); |
| remup(g,f,&t); g = t; |
remup(g,f,&t); g = t; |
|
|
| *xp = s; |
*xp = s; |
| } |
} |
| |
|
| void square_rem_tab_up_gf2n(f,tab,rp) |
void square_rem_tab_up_gf2n(UP f,UP *tab,UP *rp) |
| UP f; |
|
| UP *tab; |
|
| UP *rp; |
|
| { |
{ |
| UP s,t,u,n; |
UP s,t,u,n; |
| Num *c; |
Num *c; |
|
|
| } |
} |
| } |
} |
| |
|
| void powertabup_gf2n(f,xp,tab) |
void powertabup_gf2n(UP f,UP xp,UP *tab) |
| UP f; |
|
| UP xp; |
|
| UP *tab; |
|
| { |
{ |
| UP y,t,invf; |
UP y,t,invf; |
| int i,d; |
int i,d; |
|
|
| } |
} |
| } |
} |
| |
|
| void find_root_gf2n(f,r) |
void find_root_gf2n(UP f,GF2N *r) |
| UP f; |
|
| GF2N *r; |
|
| { |
{ |
| UP g,ut,c,t,h,rem; |
UP g,ut,c,t,h,rem; |
| int n; |
int n; |
| GF2N rn; |
GF2N rn; |
| struct oEGT eg0,eg1,eg_trace; |
|
| |
|
| n = degup2(current_mod_gf2n->dense); |
n = degup2(current_mod_gf2n->dense); |
| g = f; |
g = f; |
![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to up_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