=================================================================== RCS file: /home/cvs/OpenXM_contrib2/asir2000/engine/lmi.c,v retrieving revision 1.1 retrieving revision 1.7 diff -u -p -r1.1 -r1.7 --- OpenXM_contrib2/asir2000/engine/lmi.c 1999/12/03 07:39:08 1.1 +++ OpenXM_contrib2/asir2000/engine/lmi.c 2020/10/04 03:14:08 1.7 @@ -1,317 +1,354 @@ -/* $OpenXM: OpenXM/src/asir99/engine/lmi.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */ +/* + * Copyright (c) 1994-2000 FUJITSU LABORATORIES LIMITED + * All rights reserved. + * + * FUJITSU LABORATORIES LIMITED ("FLL") hereby grants you a limited, + * non-exclusive and royalty-free license to use, copy, modify and + * redistribute, solely for non-commercial and non-profit purposes, the + * computer program, "Risa/Asir" ("SOFTWARE"), subject to the terms and + * conditions of this Agreement. For the avoidance of doubt, you acquire + * only a limited right to use the SOFTWARE hereunder, and FLL or any + * third party developer retains all rights, including but not limited to + * copyrights, in and to the SOFTWARE. + * + * (1) FLL does not grant you a license in any way for commercial + * purposes. You may use the SOFTWARE only for non-commercial and + * non-profit purposes only, such as academic, research and internal + * business use. + * (2) The SOFTWARE is protected by the Copyright Law of Japan and + * international copyright treaties. If you make copies of the SOFTWARE, + * with or without modification, as permitted hereunder, you shall affix + * to all such copies of the SOFTWARE the above copyright notice. + * (3) An explicit reference to this SOFTWARE and its copyright owner + * shall be made on your publication or presentation in any form of the + * results obtained by use of the SOFTWARE. + * (4) In the event that you modify the SOFTWARE, you shall notify FLL by + * e-mail at risa-admin@sec.flab.fujitsu.co.jp of the detailed specification + * for such modification or the source code of the modified part of the + * SOFTWARE. + * + * THE SOFTWARE IS PROVIDED AS IS WITHOUT ANY WARRANTY OF ANY KIND. FLL + * MAKES ABSOLUTELY NO WARRANTIES, EXPRESSED, IMPLIED OR STATUTORY, AND + * EXPRESSLY DISCLAIMS ANY IMPLIED WARRANTY OF MERCHANTABILITY, FITNESS + * FOR A PARTICULAR PURPOSE OR NONINFRINGEMENT OF THIRD PARTIES' + * RIGHTS. NO FLL DEALER, AGENT, EMPLOYEES IS AUTHORIZED TO MAKE ANY + * MODIFICATIONS, EXTENSIONS, OR ADDITIONS TO THIS WARRANTY. + * UNDER NO CIRCUMSTANCES AND UNDER NO LEGAL THEORY, TORT, CONTRACT, + * OR OTHERWISE, SHALL FLL BE LIABLE TO YOU OR ANY OTHER PERSON FOR ANY + * DIRECT, INDIRECT, SPECIAL, INCIDENTAL, PUNITIVE OR CONSEQUENTIAL + * DAMAGES OF ANY CHARACTER, INCLUDING, WITHOUT LIMITATION, DAMAGES + * ARISING OUT OF OR RELATING TO THE SOFTWARE OR THIS AGREEMENT, DAMAGES + * FOR LOSS OF GOODWILL, WORK STOPPAGE, OR LOSS OF DATA, OR FOR ANY + * DAMAGES, EVEN IF FLL SHALL HAVE BEEN INFORMED OF THE POSSIBILITY OF + * SUCH DAMAGES, OR FOR ANY CLAIM BY ANY OTHER PARTY. EVEN IF A PART + * OF THE SOFTWARE HAS BEEN DEVELOPED BY A THIRD PARTY, THE THIRD PARTY + * DEVELOPER SHALL HAVE NO LIABILITY IN CONNECTION WITH THE USE, + * PERFORMANCE OR NON-PERFORMANCE OF THE SOFTWARE. + * + * $OpenXM: OpenXM_contrib2/asir2000/engine/lmi.c,v 1.6 2018/03/29 01:32:52 noro Exp $ +*/ #include "ca.h" #include "base.h" #include "inline.h" typedef struct oGEN_LM { - int id; - N dense; - N sparse; - int bit; - unsigned int lower; + int id; + N dense; + N sparse; + int bit; + unsigned int lower; } *GEN_LM; void pwrlm0(N,N,N *); void gen_simpn(N,N*); void gen_simpn_force(N,N*); +void setmod_lf(N p); int lm_lazy; static GEN_LM current_mod_lm; -void random_lm(r) -LM *r; +void random_lm(LM *r) { - N n,t; + N n,t; - if ( !current_mod_lm ) - error("random_lm : current_mod_lm is not set"); - randomn(current_mod_lm->bit+1,&n); - gen_simpn_force(n,&t); - MKLM(t,*r); + if ( !current_mod_lm ) + error("random_lm : current_mod_lm is not set"); + randomn(current_mod_lm->bit+1,&n); + gen_simpn_force(n,&t); + MKLM(t,*r); } -void ntosparsen(p,bits) -N p; -N *bits; +void ntosparsen(N p,N *bits) { - int l,i,j,nz; - N r; - unsigned int *pb,*rb; - if ( !p ) - *bits = 0; - else { - l = n_bits(p); - pb = BD(p); - for ( i = l-1, nz = 0; i >= 0; i-- ) - if ( pb[i>>5]&(1<= 0; i-- ) - if ( pb[i>>5]&(1<= 0; i-- ) + if ( pb[i>>5]&(1<= 0; i-- ) + if ( pb[i>>5]&(1<dense = p; - ntosparsen(p,¤t_mod_lm->sparse); - current_mod_lm->bit = n_bits(p); - current_mod_lm->id = UP2_DENSE; /* use ordinary rep. by default */ - if ( !(current_mod_lm->bit%32) ) { - for ( i = PL(p)-1; i >= 1; i-- ) - if ( BD(p)[i] != 0xffffffff ) - break; - if ( !i ) { - current_mod_lm->lower = (unsigned int)((((UL)1)<<32)-((UL)BD(p)[0])); - current_mod_lm->id = UP2_SPARSE; /* XXX */ - } - } + if ( !current_mod_lm ) + current_mod_lm = (GEN_LM)MALLOC(sizeof(struct oGEN_LM)); + bzero((char *)current_mod_lm,sizeof(struct oGEN_LM)); + current_mod_lm->dense = p; + ntosparsen(p,¤t_mod_lm->sparse); + current_mod_lm->bit = n_bits(p); + current_mod_lm->id = UP2_DENSE; /* use ordinary rep. by default */ + if ( !(current_mod_lm->bit%32) ) { + for ( i = PL(p)-1; i >= 1; i-- ) + if ( BD(p)[i] != 0xffffffff ) + break; + if ( !i ) { + current_mod_lm->lower = (unsigned int)((((UL)1)<<32)-((UL)BD(p)[0])); + current_mod_lm->id = UP2_SPARSE; /* XXX */ + } + } + setmod_lf(p); } -void getmod_lm(p) -N *p; +void getmod_lm(N *p) { - if ( !current_mod_lm ) - *p = 0; - else - *p = current_mod_lm->dense; + if ( !current_mod_lm ) + *p = 0; + else + *p = current_mod_lm->dense; } -void simplm(n,r) -LM n; -LM *r; +void simplm(LM n,LM *r) { - N rem; + N rem; - if ( !n ) - *r = 0; - else if ( NID(n) != N_LM ) - *r = n; - else { - gen_simpn(n->body,&rem); - MKLM(rem,*r); - } + if ( !n ) + *r = 0; + else if ( NID(n) != N_LM ) + *r = n; + else { + gen_simpn(n->body,&rem); + MKLM(rem,*r); + } } -void qtolm(q,l) -Q q; -LM *l; +void simplm_force(LM n,LM *r) { - N rn; - LM a,b,c; + N rem; - if ( !q || (OID(q)==O_N && ((NID(q) == N_LM) || (NID(q) == N_GFPN))) ) { /* XXX */ - *l = (LM)q; - } else if ( OID(q) == O_N && NID(q) == N_Q ) { - gen_simpn(NM(q),&rn); MKLM(rn,a); - if ( SGN(q) < 0 ) { - chsgnlm(a,&c); a = c; - } - if ( DN(q) ) { - gen_simpn_force(DN(q),&rn); MKLM(rn,b); - if ( !b ) - error("qtolm : invalid denominator"); - divlm(a,b,l); - } else - *l = a; - } else - error("qtolm : invalid argument"); + if ( !n ) + *r = 0; + else if ( NID(n) != N_LM ) + *r = n; + else { + gen_simpn_force(n->body,&rem); + MKLM(rem,*r); + } } + +void qtolm(Q q,LM *l) +{ + N rn; + LM a,b,c; + + if ( !q || (OID(q)==O_N && ((NID(q) == N_LM) || (NID(q) == N_GFPN))) ) { /* XXX */ + *l = (LM)q; + } else if ( OID(q) == O_N && NID(q) == N_Q ) { + gen_simpn(NM(q),&rn); MKLM(rn,a); + if ( SGN(q) < 0 ) { + chsgnlm(a,&c); a = c; + } + if ( DN(q) ) { + gen_simpn_force(DN(q),&rn); MKLM(rn,b); + if ( !b ) + error("qtolm : invalid denominator"); + divlm(a,b,l); + } else + *l = a; + } else + error("qtolm : invalid argument"); +} + #define NZLM(a) ((a)&&(NID(a)==N_LM)) -void addlm(a,b,c) -LM a,b; -LM *c; +void addlm(LM a,LM b,LM *c) { - int s; - N t,t1; - LM z; + N t,t1; + LM z; - qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; - if ( !a ) - *c = b; - else if ( !b ) - *c = a; - else { - addn(a->body,b->body,&t); - gen_simpn(t,&t1); - MKLM(t1,*c); - } + qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; + if ( !a ) + *c = b; + else if ( !b ) + *c = a; + else { + addn(a->body,b->body,&t); + gen_simpn(t,&t1); + MKLM(t1,*c); + } } -void sublm(a,b,c) -LM a,b; -LM *c; +void sublm(LM a,LM b,LM *c) { - int s; - N t,t1,mod; - LM z; + int s; + N t; + LM z; - qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; - if ( !b ) - *c = a; - else if ( !a ) - chsgnlm(b,c); - else { - s = subn(a->body,b->body,&t); - if ( !s ) - *c = 0; - else { - MKLM(t,z); - if ( s < 0 ) - chsgnlm(z,c); - else - *c = z; - } - } + qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; + if ( !b ) + *c = a; + else if ( !a ) + chsgnlm(b,c); + else { + s = subn(a->body,b->body,&t); + if ( !s ) + *c = 0; + else { + MKLM(t,z); + if ( s < 0 ) + chsgnlm(z,c); + else + *c = z; + } + } } -void mullm(a,b,c) -LM a,b; -LM *c; +void mullm(LM a,LM b,LM *c) { - N t,q,r; - LM z; + N t,r; + LM z; - qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; - if ( !a || !b ) - *c = 0; - else { - kmuln(a->body,b->body,&t); - gen_simpn(t,&r); - MKLM(r,*c); - } + qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; + if ( !a || !b ) + *c = 0; + else { + kmuln(a->body,b->body,&t); + gen_simpn(t,&r); + MKLM(r,*c); + } } -void divlm(a,b,c) -LM a,b; -LM *c; +void divlm(LM a,LM b,LM *c) { - LM r; - Q t,m,i; - LM z; + LM r; + Q t,m,i; + LM z; - qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; - if ( !b ) - error("divlm: division by 0"); - else if ( !a ) - *c = 0; - else { - NTOQ(b->body,1,t); NTOQ(current_mod_lm->dense,1,m); - invl(t,m,&i); - MKLM(NM(i),r); - mullm(a,r,c); - } + qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; + if ( !b ) + error("divlm: division by 0"); + else if ( !a ) + *c = 0; + else { + NTOQ(b->body,1,t); NTOQ(current_mod_lm->dense,1,m); + invl(t,m,&i); + MKLM(NM(i),r); + mullm(a,r,c); + } } -void chsgnlm(a,c) -LM a,*c; +void chsgnlm(LM a,LM *c) { - LM t; - N s,u; + LM t; + N s,u; - if ( !a ) - *c = a; - else { - qtolm((Q)a,&t); a = t; - gen_simpn_force(a->body,&s); - subn(current_mod_lm->dense,s,&u); - MKLM(u,*c); - } + if ( !a ) + *c = a; + else { + qtolm((Q)a,&t); a = t; + gen_simpn_force(a->body,&s); + subn(current_mod_lm->dense,s,&u); + MKLM(u,*c); + } } -void pwrlm(a,b,c) -LM a; -Q b; -LM *c; +void pwrlm(LM a,Q b,LM *c) { - LM t; - N s; + LM t; + N s; - if ( !b ) { - MKLM(ONEN,*c); - } else if ( !a ) - *c = 0; - else { - qtolm((Q)a,&t); a = t; - pwrlm0(a->body,NM(b),&s); - MKLM(s,*c); - } + if ( !b ) { + MKLM(ONEN,*c); + } else if ( !a ) + *c = 0; + else { + qtolm((Q)a,&t); a = t; + pwrlm0(a->body,NM(b),&s); + MKLM(s,*c); + } } -void pwrlm0(a,n,c) -N a,n; -N *c; +void pwrlm0(N a,N n,N *c) { - N n1,t,t1,t2,r1; - int r; + N n1,t,t1,t2,r1; + int r; - if ( !n ) - *c = ONEN; - else if ( UNIN(n) ) - *c = a; - else { - r = divin(n,2,&n1); - pwrlm0(a,n1,&t); - kmuln(t,t,&t1); gen_simpn(t1,&r1); - if ( r ) { - kmuln(r1,a,&t2); gen_simpn(t2,&r1); - } - *c = r1; - } + if ( !n ) + *c = ONEN; + else if ( UNIN(n) ) + *c = a; + else { + r = divin(n,2,&n1); + pwrlm0(a,n1,&t); + kmuln(t,t,&t1); gen_simpn(t1,&r1); + if ( r ) { + kmuln(r1,a,&t2); gen_simpn(t2,&r1); + } + *c = r1; + } } -int cmplm(a,b) -LM a,b; +int cmplm(LM a,LM b) { - LM z; + LM z; - qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; - if ( !a ) - if ( !b ) - return 0; - else - return -1; - else if ( !b ) - return 1; - else - return cmpn(a->body,b->body); + qtolm((Q)a,&z); a = z; qtolm((Q)b,&z); b = z; + if ( !a ) + if ( !b ) + return 0; + else + return -1; + else if ( !b ) + return 1; + else + return cmpn(a->body,b->body); } void remn_special(N,N,int,unsigned int ,N *); -void gen_simpn(a,b) -N a; -N *b; +void gen_simpn(N a,N *b) { - if ( !current_mod_lm ) - error("gen_simpn: current_mod_lm is not set"); - if ( lm_lazy ) - *b = a; - else if ( current_mod_lm->id == UP2_SPARSE ) - remn_special(a,current_mod_lm->dense,current_mod_lm->bit,current_mod_lm->lower,b); - else - remn(a,current_mod_lm->dense,b); + if ( !current_mod_lm ) + error("gen_simpn: current_mod_lm is not set"); + if ( lm_lazy ) + *b = a; + else if ( current_mod_lm->id == UP2_SPARSE ) + remn_special(a,current_mod_lm->dense,current_mod_lm->bit,current_mod_lm->lower,b); + else + remn(a,current_mod_lm->dense,b); } -void gen_simpn_force(a,b) -N a; -N *b; +void gen_simpn_force(N a,N *b) { - if ( !current_mod_lm ) - error("gen_simpn_force: current_mod_lm is not set"); - else if ( current_mod_lm->id == UP2_SPARSE ) - remn_special(a,current_mod_lm->dense,current_mod_lm->bit,current_mod_lm->lower,b); - else - remn(a,current_mod_lm->dense,b); + if ( !current_mod_lm ) + error("gen_simpn_force: current_mod_lm is not set"); + else if ( current_mod_lm->id == UP2_SPARSE ) + remn_special(a,current_mod_lm->dense,current_mod_lm->bit,current_mod_lm->lower,b); + else + remn(a,current_mod_lm->dense,b); }