version 1.3, 2000/08/22 05:04:03 |
version 1.5, 2001/10/09 01:36:09 |
|
|
* 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/A.c,v 1.2 2000/08/21 08:31:24 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/A.c,v 1.4 2001/06/07 04:54:39 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
|
|
|
int get_lprime(); |
|
|
void pdiva(vl,p0,p1,p2,pr) |
void pdiva(vl,p0,p1,p2,pr) |
VL vl; |
VL vl; |
P p1,p2,p0; |
P p1,p2,p0; |
|
|
wt = W_UMALLOC(m+n); ws = W_UMALLOC(m+n); |
wt = W_UMALLOC(m+n); ws = W_UMALLOC(m+n); |
wm = W_UMALLOC(m+n); |
wm = W_UMALLOC(m+n); |
for ( q = ONE, t = 0, c = 0, index = 0; ; ) { |
for ( q = ONE, t = 0, c = 0, index = 0; ; ) { |
mod = lprime[index++]; |
mod = get_lprime(index++); |
if ( !mod ) |
|
error("pinva : lprime[] exhausted."); |
|
if ( !rem(NM((Q)LC(p)),mod) || !rem(NM((Q)LC(p0)),mod) ) |
if ( !rem(NM((Q)LC(p)),mod) || !rem(NM((Q)LC(p0)),mod) ) |
continue; |
continue; |
ptomp(mod,p,&tg); ptomp(mod,p0,&th); srchump(mod,tg,th,&tr); |
ptomp(mod,p,&tg); ptomp(mod,p0,&th); srchump(mod,tg,th,&tr); |