version 1.85, 2003/10/28 03:23:40 |
version 1.90, 2004/03/13 06:49:15 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.84 2003/10/18 01:39:42 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.89 2004/02/03 23:31:57 noro Exp $ */ |
|
|
#include "ca.h" |
#include "ca.h" |
#include "parse.h" |
#include "parse.h" |
#include "ox.h" |
#include "ox.h" |
|
#include "base.h" |
#include "inline.h" |
#include "inline.h" |
#include <time.h> |
#include <time.h> |
|
|
Line 117 typedef struct oIndArray |
|
Line 118 typedef struct oIndArray |
|
int (*ndl_compare_function)(UINT *a1,UINT *a2); |
int (*ndl_compare_function)(UINT *a1,UINT *a2); |
|
|
static int ndv_alloc; |
static int ndv_alloc; |
|
#if 1 |
static int nd_f4_nsp=0x7fffffff; |
static int nd_f4_nsp=0x7fffffff; |
|
#else |
|
static int nd_f4_nsp=50; |
|
#endif |
static double nd_scale=2; |
static double nd_scale=2; |
static UINT **nd_bound; |
static UINT **nd_bound; |
static struct order_spec *nd_ord; |
static struct order_spec *nd_ord; |
Line 141 static int nd_psn,nd_pslen; |
|
Line 146 static int nd_psn,nd_pslen; |
|
static RHist *nd_red; |
static RHist *nd_red; |
|
|
static int nd_found,nd_create,nd_notfirst; |
static int nd_found,nd_create,nd_notfirst; |
static int nm_adv; |
|
static int nmv_adv; |
static int nmv_adv; |
static int nd_dcomp; |
static int nd_dcomp; |
static int nd_demand; |
static int nd_demand; |
|
|
extern struct order_spec dp_current_spec; |
extern struct order_spec *dp_current_spec; |
extern char *Demand; |
extern char *Demand; |
extern VL CO; |
extern VL CO; |
extern int Top,Reverse,DP_Print,dp_nelim,do_weyl,NoSugar; |
extern int Top,Reverse,DP_Print,dp_nelim,do_weyl,NoSugar; |
Line 174 extern int *current_weyl_weight_vector; |
|
Line 178 extern int *current_weyl_weight_vector; |
|
#define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0) |
#define GET_EXP_OLD(d,a) (((d)[oepos[a].i]>>oepos[a].s)&omask0) |
#define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s)) |
#define PUT_EXP_OLD(r,a,e) ((r)[oepos[a].i] |= ((e)<<oepos[a].s)) |
|
|
|
#define ROUND_FOR_ALIGN(s) ((((s)+sizeof(void *)-1)/sizeof(void *))*sizeof(void *)) |
|
|
/* macros for term comparison */ |
/* macros for term comparison */ |
#define TD_DL_COMPARE(d1,d2)\ |
#define TD_DL_COMPARE(d1,d2)\ |
(TD(d1)>TD(d2)?1:(TD(d1)<TD(d2)?-1:ndl_lex_compare(d1,d2))) |
(TD(d1)>TD(d2)?1:(TD(d1)<TD(d2)?-1:ndl_lex_compare(d1,d2))) |
Line 2401 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
Line 2407 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
ndv_alloc = 0; |
ndv_alloc = 0; |
get_vars((Obj)f,&fv); pltovl(v,&vv); |
get_vars((Obj)f,&fv); pltovl(v,&vv); |
for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ ); |
for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ ); |
|
switch ( ord->id ) { |
|
case 1: |
|
if ( ord->nv != nvar ) |
|
error("nd_{gr,f4} : invalid order specification"); |
|
break; |
|
default: |
|
break; |
|
} |
nd_init_ord(ord); |
nd_init_ord(ord); |
for ( t = BDY(f), max = 0; t; t = NEXT(t) ) |
for ( t = BDY(f), max = 0; t; t = NEXT(t) ) |
for ( tv = vv; tv; tv = NEXT(tv) ) { |
for ( tv = vv; tv; tv = NEXT(tv) ) { |
Line 2433 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
Line 2447 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
|
|
void nd_gr_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp) |
void nd_gr_trace(LIST f,LIST v,int trace,int homo,struct order_spec *ord,LIST *rp) |
{ |
{ |
struct order_spec ord1; |
struct order_spec *ord1; |
VL tv,fv,vv,vc; |
VL tv,fv,vv,vc; |
NODE fd,fd0,in0,in,r,r0,t,s,cand; |
NODE fd,fd0,in0,in,r,r0,t,s,cand; |
int m,nocheck,nvar,mindex,e,max; |
int m,nocheck,nvar,mindex,e,max; |
Line 2445 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
Line 2459 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
|
|
get_vars((Obj)f,&fv); pltovl(v,&vv); |
get_vars((Obj)f,&fv); pltovl(v,&vv); |
for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ ); |
for ( nvar = 0, tv = vv; tv; tv = NEXT(tv), nvar++ ); |
|
switch ( ord->id ) { |
|
case 1: |
|
if ( ord->nv != nvar ) |
|
error("nd_gr_trace : invalid order specification"); |
|
break; |
|
default: |
|
break; |
|
} |
nocheck = 0; |
nocheck = 0; |
mindex = 0; |
mindex = 0; |
|
|
Line 2484 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
Line 2506 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
wmax = MAX(TD(DL(a)),wmax); |
wmax = MAX(TD(DL(a)),wmax); |
} |
} |
homogenize_order(ord,nvar,&ord1); |
homogenize_order(ord,nvar,&ord1); |
nd_init_ord(&ord1); |
nd_init_ord(ord1); |
nd_setup_parameters(nvar+1,wmax); |
nd_setup_parameters(nvar+1,wmax); |
for ( t = fd0; t; t = NEXT(t) ) |
for ( t = fd0; t; t = NEXT(t) ) |
ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos); |
ndv_homogenize((NDV)BDY(t),obpe,oadv,oepos); |
Line 2523 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
Line 2545 void nd_gr_trace(LIST f,LIST v,int trace,int homo,stru |
|
m = get_lprime(++mindex); |
m = get_lprime(++mindex); |
/* reset the parameters */ |
/* reset the parameters */ |
if ( !ishomo && homo ) { |
if ( !ishomo && homo ) { |
nd_init_ord(&ord1); |
nd_init_ord(ord1); |
nd_setup_parameters(nvar+1,wmax); |
nd_setup_parameters(nvar+1,wmax); |
} else { |
} else { |
nd_init_ord(ord); |
nd_init_ord(ord); |
Line 2754 void ndv_dehomogenize(NDV p,struct order_spec *ord) |
|
Line 2776 void ndv_dehomogenize(NDV p,struct order_spec *ord) |
|
for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) |
for ( m = BDY(p), i = 0; i < len; NMV_ADV(m), i++ ) |
ndl_dehomogenize(DL(m)); |
ndl_dehomogenize(DL(m)); |
if ( newwpd != nd_wpd ) { |
if ( newwpd != nd_wpd ) { |
newadv = sizeof(struct oNMV)+(newwpd-1)*sizeof(UINT); |
newadv = ROUND_FOR_ALIGN(sizeof(struct oNMV)+(newwpd-1)*sizeof(UINT)); |
for ( m = r = BDY(p), i = 0; i < len; NMV_ADV(m), NDV_NADV(r), i++ ) { |
for ( m = r = BDY(p), i = 0; i < len; NMV_ADV(m), NDV_NADV(r), i++ ) { |
CQ(r) = CQ(m); |
CQ(r) = CQ(m); |
for ( j = 0; j < newexporigin; j++ ) DL(r)[j] = DL(m)[j]; |
for ( j = 0; j < newexporigin; j++ ) DL(r)[j] = DL(m)[j]; |
Line 2950 void nd_setup_parameters(int nvar,int max) { |
|
Line 2972 void nd_setup_parameters(int nvar,int max) { |
|
nd_mask[nd_epw-i-1] = (nd_mask0<<(i*nd_bpe)); |
nd_mask[nd_epw-i-1] = (nd_mask0<<(i*nd_bpe)); |
nd_mask1 |= (1<<(nd_bpe-1))<<(i*nd_bpe); |
nd_mask1 |= (1<<(nd_bpe-1))<<(i*nd_bpe); |
} |
} |
nm_adv = sizeof(struct oNM)+(nd_wpd-1)*sizeof(UINT); |
nmv_adv = ROUND_FOR_ALIGN(sizeof(struct oNMV)+(nd_wpd-1)*sizeof(UINT)); |
nmv_adv = sizeof(struct oNMV)+(nd_wpd-1)*sizeof(UINT); |
|
nd_epos = nd_create_epos(nd_ord); |
nd_epos = nd_create_epos(nd_ord); |
nd_blockmask = nd_create_blockmask(nd_ord); |
nd_blockmask = nd_create_blockmask(nd_ord); |
} |
} |
Line 4183 NODE nd_f4(int m) |
|
Line 4204 NODE nd_f4(int m) |
|
d = nd_reconstruct(m,0,d); |
d = nd_reconstruct(m,0,d); |
continue; |
continue; |
} |
} |
|
if ( bucket->m < 0 ) continue; |
col = nd_symbolic_preproc(bucket,&s0vect,&rp0); |
col = nd_symbolic_preproc(bucket,&s0vect,&rp0); |
if ( !col ) { |
if ( !col ) { |
for ( t = l; NEXT(t); t = NEXT(t) ); |
for ( t = l; NEXT(t); t = NEXT(t) ); |
Line 4411 NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,in |
|
Line 4433 NODE nd_f4_red_dist(int m,ND_pairs sp0,UINT *s0vect,in |
|
nd_send_int(nd_wpd); |
nd_send_int(nd_wpd); |
nd_send_int(nmv_adv); |
nd_send_int(nmv_adv); |
|
|
saveobj(nd_write,dp_current_spec.obj); fflush(nd_write); |
saveobj(nd_write,dp_current_spec->obj); fflush(nd_write); |
|
|
nd_send_int(nd_psn); |
nd_send_int(nd_psn); |
for ( i = 0; i < nd_psn; i++ ) nd_send_ndv(nd_ps[i]); |
for ( i = 0; i < nd_psn; i++ ) nd_send_ndv(nd_ps[i]); |
Line 4457 void nd_exec_f4_red_dist() |
|
Line 4479 void nd_exec_f4_red_dist() |
|
ND_pairs *sp0; |
ND_pairs *sp0; |
int *colstat; |
int *colstat; |
int a,sprow,rank; |
int a,sprow,rank; |
struct order_spec ord; |
struct order_spec *ord; |
Obj ordspec; |
Obj ordspec; |
ND spol; |
ND spol; |
int maxrs; |
int maxrs; |
Line 4472 void nd_exec_f4_red_dist() |
|
Line 4494 void nd_exec_f4_red_dist() |
|
nmv_adv = nd_recv_int(); |
nmv_adv = nd_recv_int(); |
|
|
loadobj(nd_read,&ordspec); |
loadobj(nd_read,&ordspec); |
create_order_spec(ordspec,&ord); |
create_order_spec(0,ordspec,&ord); |
nd_init_ord(&ord); |
nd_init_ord(ord); |
nd_setup_parameters(nd_nvar,0); |
nd_setup_parameters(nd_nvar,0); |
|
|
nd_psn = nd_recv_int(); |
nd_psn = nd_recv_int(); |