version 1.39, 2005/08/25 18:59:11 |
version 1.40, 2006/12/12 11:50:37 |
|
|
* 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/builtin/dp-supp.c,v 1.38 2004/12/06 09:29:34 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp-supp.c,v 1.39 2005/08/25 18:59:11 ohara Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
Line 511 void heu_nezgcdnpz(VL vl,P *pl,int m,P *pr) |
|
Line 511 void heu_nezgcdnpz(VL vl,P *pl,int m,P *pr) |
|
int i,r; |
int i,r; |
P gcd,t,s1,s2,u; |
P gcd,t,s1,s2,u; |
Q rq; |
Q rq; |
|
DCP dc; |
|
extern int DP_Print; |
|
|
while ( 1 ) { |
while ( 1 ) { |
for ( i = 0, s1 = 0; i < m; i++ ) { |
for ( i = 0, s1 = 0; i < m; i++ ) { |
r = random(); UTOQ(r,rq); |
r = random(); UTOQ(r,rq); |
Line 522 void heu_nezgcdnpz(VL vl,P *pl,int m,P *pr) |
|
Line 524 void heu_nezgcdnpz(VL vl,P *pl,int m,P *pr) |
|
mulp(vl,pl[i],(P)rq,&t); addp(vl,s2,t,&u); s2 = u; |
mulp(vl,pl[i],(P)rq,&t); addp(vl,s2,t,&u); s2 = u; |
} |
} |
ezgcdp(vl,s1,s2,&gcd); |
ezgcdp(vl,s1,s2,&gcd); |
|
if ( DP_Print > 2 ) |
|
{ fprintf(asir_out,"(%d)",nmonop(gcd)); fflush(asir_out); } |
for ( i = 0; i < m; i++ ) { |
for ( i = 0; i < m; i++ ) { |
if ( !divtpz(vl,pl[i],gcd,&t) ) |
if ( !divtpz(vl,pl[i],gcd,&t) ) |
break; |
break; |