version 1.25, 2002/06/06 01:18:05 |
version 1.27, 2003/01/06 01:16: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.c,v 1.24 2002/01/28 00:54:41 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/builtin/dp.c,v 1.26 2003/01/04 09:06:15 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "base.h" |
#include "base.h" |
#include "parse.h" |
#include "parse.h" |
|
|
extern int dp_fcoeffs; |
|
extern int dp_nelim; |
extern int dp_nelim; |
extern int dp_order_pair_length; |
extern int dp_order_pair_length; |
extern struct order_pair *dp_order_pair; |
extern struct order_pair *dp_order_pair; |
|
|
{ |
{ |
LIST f,v; |
LIST f,v; |
Num homo; |
Num homo; |
|
int m,field,t; |
struct order_spec ord; |
struct order_spec ord; |
|
NODE n; |
|
|
do_weyl = 0; |
do_weyl = 0; |
asir_assert(ARG0(arg),O_LIST,"dp_gr_f_main"); |
asir_assert(ARG0(arg),O_LIST,"dp_gr_f_main"); |
|
|
*rp = f; return; |
*rp = f; return; |
} |
} |
homo = (Num)ARG2(arg); |
homo = (Num)ARG2(arg); |
|
#if 0 |
|
asir_assert(ARG3(arg),O_N,"dp_gr_f_main"); |
|
m = QTOS((Q)ARG3(arg)); |
|
if ( m ) |
|
error("dp_gr_f_main : trace lifting is not implemented yet"); |
|
create_order_spec(ARG4(arg),&ord); |
|
#else |
|
m = 0; |
create_order_spec(ARG3(arg),&ord); |
create_order_spec(ARG3(arg),&ord); |
dp_gr_main(f,v,homo,0,1,&ord,rp); |
#endif |
|
field = 0; |
|
for ( n = BDY(f); n; n = NEXT(n) ) { |
|
t = get_field_type(BDY(n)); |
|
if ( !t ) |
|
continue; |
|
if ( t < 0 ) |
|
error("dp_gr_f_main : incosistent coefficients"); |
|
if ( !field ) |
|
field = t; |
|
else if ( t != field ) |
|
error("dp_gr_f_main : incosistent coefficients"); |
|
} |
|
dp_gr_main(f,v,homo,m?1:0,field,&ord,rp); |
} |
} |
|
|
void Pdp_f4_main(arg,rp) |
void Pdp_f4_main(arg,rp) |
Line 1572 LIST remove_zero_from_list(LIST l) |
|
Line 1594 LIST remove_zero_from_list(LIST l) |
|
NEXT(r) = 0; |
NEXT(r) = 0; |
MKLIST(rl,r0); |
MKLIST(rl,r0); |
return rl; |
return rl; |
|
} |
|
|
|
int get_field_type(P p) |
|
{ |
|
int type,t; |
|
DCP dc; |
|
|
|
if ( !p ) |
|
return 0; |
|
else if ( NUM(p) ) |
|
return NID((Num)p); |
|
else { |
|
type = 0; |
|
for ( dc = DC(p); dc; dc = NEXT(dc) ) { |
|
t = get_field_type(COEF(dc)); |
|
if ( !t ) |
|
continue; |
|
if ( t < 0 ) |
|
return t; |
|
if ( !type ) |
|
type = t; |
|
else if ( t != type ) |
|
return -1; |
|
} |
|
return type; |
|
} |
} |
} |