version 1.86, 2003/10/28 08:47:01 |
version 1.89, 2004/02/03 23:31:57 |
|
|
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.85 2003/10/28 03:23:40 noro Exp $ */ |
/* $OpenXM: OpenXM_contrib2/asir2000/engine/nd.c,v 1.88 2003/12/24 02:20:19 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 146 static int nmv_adv; |
|
Line 151 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 2401 void nd_gr(LIST f,LIST v,int m,int f4,struct order_spe |
|
Line 2406 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 2446 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 2458 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 2505 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 2544 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 4412 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 4458 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 4473 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(); |