version 1.159, 2009/01/05 00:52:20 |
version 1.162, 2009/01/05 02:47:39 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.158 2009/01/04 10:02:00 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.161 2009/01/05 02:08:18 noro Exp $ */ |
|
|
#include "nd.h" |
#include "nd.h" |
|
|
Line 46 static int nd_worb_len; |
|
Line 46 static int nd_worb_len; |
|
static int nd_found,nd_create,nd_notfirst; |
static int nd_found,nd_create,nd_notfirst; |
static int nmv_adv; |
static int nmv_adv; |
static int nd_demand; |
static int nd_demand; |
static int nd_module,nd_istop,nd_mpos; |
static int nd_module,nd_ispot,nd_mpos; |
|
|
NumberField get_numberfield(); |
NumberField get_numberfield(); |
UINT *nd_det_compute_bound(NDV **dm,int n,int j); |
UINT *nd_det_compute_bound(NDV **dm,int n,int j); |
Line 399 int ndl_weight(UINT *d) |
|
Line 399 int ndl_weight(UINT *d) |
|
for ( j = 0; j < nd_epw; j++, u>>=nd_bpe ) |
for ( j = 0; j < nd_epw; j++, u>>=nd_bpe ) |
t += (u&nd_mask0); |
t += (u&nd_mask0); |
} |
} |
|
if ( nd_module && current_module_weight_vector && MPOS(d) ) |
|
t += current_module_weight_vector[MPOS(d)]; |
return t; |
return t; |
} |
} |
|
|
Line 568 int ndl_module_grlex_compare(UINT *d1,UINT *d2) |
|
Line 570 int ndl_module_grlex_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i; |
int i; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 577 int ndl_module_grlex_compare(UINT *d1,UINT *d2) |
|
Line 579 int ndl_module_grlex_compare(UINT *d1,UINT *d2) |
|
for ( i = nd_nvar-1; i >= 0; i-- ) |
for ( i = nd_nvar-1; i >= 0; i-- ) |
if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return 1; |
if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return 1; |
else if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return -1; |
else if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return -1; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 588 int ndl_module_glex_compare(UINT *d1,UINT *d2) |
|
Line 590 int ndl_module_glex_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i; |
int i; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 597 int ndl_module_glex_compare(UINT *d1,UINT *d2) |
|
Line 599 int ndl_module_glex_compare(UINT *d1,UINT *d2) |
|
for ( i = 0; i < nd_nvar; i++ ) |
for ( i = 0; i < nd_nvar; i++ ) |
if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return 1; |
if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return 1; |
else if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return -1; |
else if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return -1; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 608 int ndl_module_lex_compare(UINT *d1,UINT *d2) |
|
Line 610 int ndl_module_lex_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i; |
int i; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
for ( i = 0; i < nd_nvar; i++ ) |
for ( i = 0; i < nd_nvar; i++ ) |
if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return 1; |
if ( GET_EXP(d1,i) > GET_EXP(d2,i) ) return 1; |
else if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return -1; |
else if ( GET_EXP(d1,i) < GET_EXP(d2,i) ) return -1; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 626 int ndl_module_block_compare(UINT *d1,UINT *d2) |
|
Line 628 int ndl_module_block_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i,c; |
int i,c; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
if ( c = ndl_block_compare(d1,d2) ) return c; |
if ( c = ndl_block_compare(d1,d2) ) return c; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 642 int ndl_module_matrix_compare(UINT *d1,UINT *d2) |
|
Line 644 int ndl_module_matrix_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i,c; |
int i,c; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
if ( c = ndl_matrix_compare(d1,d2) ) return c; |
if ( c = ndl_matrix_compare(d1,d2) ) return c; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
if ( MPOS(d1) < MPOS(d2) ) return 1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
else if ( MPOS(d1) > MPOS(d2) ) return -1; |
} |
} |
Line 658 int ndl_module_composite_compare(UINT *d1,UINT *d2) |
|
Line 660 int ndl_module_composite_compare(UINT *d1,UINT *d2) |
|
{ |
{ |
int i,c; |
int i,c; |
|
|
if ( nd_istop ) { |
if ( nd_ispot ) { |
if ( MPOS(d1) > MPOS(d2) ) return 1; |
if ( MPOS(d1) > MPOS(d2) ) return 1; |
else if ( MPOS(d1) < MPOS(d2) ) return -1; |
else if ( MPOS(d1) < MPOS(d2) ) return -1; |
} |
} |
if ( c = ndl_composite_compare(d1,d2) ) return c; |
if ( c = ndl_composite_compare(d1,d2) ) return c; |
if ( !nd_istop ) { |
if ( !nd_ispot ) { |
if ( MPOS(d1) > MPOS(d2) ) return 1; |
if ( MPOS(d1) > MPOS(d2) ) return 1; |
else if ( MPOS(d1) < MPOS(d2) ) return -1; |
else if ( MPOS(d1) < MPOS(d2) ) return -1; |
} |
} |
Line 725 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
Line 727 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
{ |
{ |
int i; |
int i; |
|
|
|
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_add : invalid operation"); |
|
} |
#if 1 |
#if 1 |
switch ( nd_wpd ) { |
switch ( nd_wpd ) { |
case 2: |
case 2: |
Line 743 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
Line 749 INLINE void ndl_add(UINT *d1,UINT *d2,UINT *d) |
|
#else |
#else |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
for ( i = 0; i < nd_wpd; i++ ) d[i] = d1[i]+d2[i]; |
#endif |
#endif |
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_add : invalid operation"); |
|
MPOS(d) = MPOS(d1); |
|
} |
|
} |
} |
|
|
/* d1 += d2 */ |
/* d1 += d2 */ |
Line 755 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
Line 756 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
{ |
{ |
int i; |
int i; |
|
|
|
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_addto : invalid operation"); |
|
} |
#if 1 |
#if 1 |
switch ( nd_wpd ) { |
switch ( nd_wpd ) { |
case 2: |
case 2: |
Line 773 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
Line 778 INLINE void ndl_addto(UINT *d1,UINT *d2) |
|
#else |
#else |
for ( i = 0; i < nd_wpd; i++ ) d1[i] += d2[i]; |
for ( i = 0; i < nd_wpd; i++ ) d1[i] += d2[i]; |
#endif |
#endif |
if ( nd_module ) { |
|
if ( MPOS(d1) && MPOS(d2) && (MPOS(d1) != MPOS(d2)) ) |
|
error("ndl_addto : invalid operation"); |
|
MPOS(d1) = MPOS(d2); |
|
} |
|
} |
} |
|
|
INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d) |
INLINE void ndl_sub(UINT *d1,UINT *d2,UINT *d) |
Line 4453 void nd_init_ord(struct order_spec *ord) |
|
Line 4453 void nd_init_ord(struct order_spec *ord) |
|
|
|
/* module order */ |
/* module order */ |
case 256: |
case 256: |
nd_istop = ord->istop; |
nd_ispot = ord->ispot; |
nd_dcomp = -1; |
nd_dcomp = -1; |
nd_isrlex = 0; |
nd_isrlex = 0; |
switch ( ord->ord.simple ) { |
switch ( ord->ord.simple ) { |