version 1.12, 2001/09/03 07:01:06 |
version 1.27, 2007/01/25 16:19:41 |
|
|
* 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/init.c,v 1.11 2001/06/07 04:54:40 noro Exp $ |
* $OpenXM: OpenXM_contrib2/asir2000/engine/init.c,v 1.26 2006/02/25 06:33:31 noro Exp $ |
*/ |
*/ |
#include "ca.h" |
#include "ca.h" |
#include "version.h" |
#include "version.h" |
Line 83 USINT VOIDobj = &oVOID; |
|
Line 83 USINT VOIDobj = &oVOID; |
|
int bigfloat; |
int bigfloat; |
int nez = 1; |
int nez = 1; |
int current_mod = 0; |
int current_mod = 0; |
|
int StrassenSize = 0; |
|
int outputstyle = 0; |
|
int Itvplot = 33; |
|
|
static int *lprime; |
static int *lprime; |
int lprime_size; |
static int lprime_size; |
|
|
#if PARI |
#if defined(PARI) |
int paristack = 1<<16; |
int paristack = 1<<16; |
#define ADDBF addbf |
#define ADDBF addbf |
#define SUBBF subbf |
#define SUBBF subbf |
Line 108 int paristack = 1<<16; |
|
Line 111 int paristack = 1<<16; |
|
|
|
#if defined(INTERVAL) |
#if defined(INTERVAL) |
int zerorewrite = 0; |
int zerorewrite = 0; |
void (*addnumt[])() = { addq, addreal, addalg, ADDBF, additvp, additvd, 0, additvf, addcplx, addmi, addlm, addgf2n, addgfpn, addgfs, addgfsn }; |
void (*addnumt[])() = { addq, addreal, addalg, ADDBF, 0, additvp, additvd, 0, additvf, addcplx, addmi, addlm, addgf2n, addgfpn, addgfs, addgfsn, adddalg }; |
void (*subnumt[])() = { subq, subreal, subalg, SUBBF, subitvp, subitvd, 0, subitvf, subcplx, submi, sublm, subgf2n, subgfpn, subgfs, subfspn }; |
void (*subnumt[])() = { subq, subreal, subalg, SUBBF, 0, subitvp, subitvd, 0, subitvf, subcplx, submi, sublm, subgf2n, subgfpn, subgfs, subgfsn, subdalg }; |
void (*mulnumt[])() = { mulq, mulreal, mulalg, MULBF, mulitvp, mulitvd, 0, mulitvf, mulcplx, mulmi, mullm, mulgf2n, mulgfpn, mulgfs, mulgfsn }; |
void (*mulnumt[])() = { mulq, mulreal, mulalg, MULBF, 0, mulitvp, mulitvd, 0, mulitvf, mulcplx, mulmi, mullm, mulgf2n, mulgfpn, mulgfs, mulgfsn, muldalg }; |
void (*divnumt[])() = { divq, divreal, divalg, DIVBF, divitvp, divitvd, 0, divitvf, divcplx, divmi, divlm, divgf2n, divgfpn, divgfs, divgfsn }; |
void (*divnumt[])() = { divq, divreal, divalg, DIVBF, 0, divitvp, divitvd, 0, divitvf, divcplx, divmi, divlm, divgf2n, divgfpn, divgfs, divgfsn, divdalg }; |
void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, PWRBF, pwritvp, pwritvd, 0, pwritvf, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn, pwrgfs, pwrgfsn }; |
void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, PWRBF, 0, pwritvp, pwritvd, 0, pwritvf, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn, pwrgfs, pwrgfsn, pwrdalg }; |
void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, CHSGNBF, chsgnitvp, chsgnitvd, 0, chsgnitvf, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn, chsgngfs , chsgngfsn}; |
void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, CHSGNBF, 0, chsgnitvp, chsgnitvd, 0, chsgnitvf, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn, chsgngfs , chsgngfsn, chsgndalg}; |
int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, CMPBF, cmpitvp, cmpitvd, 0, cmpitvf, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn, cmpgfs, cmpgfsn }; |
int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, CMPBF, 0, cmpitvp, cmpitvd, 0, cmpitvf, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn, cmpgfs, cmpgfsn, cmpdalg }; |
#else |
#else |
void (*addnumt[])() = { addq, addreal, addalg, ADDBF, addcplx, addmi, addlm, addgf2n, addgfpn, addgfs, addgfsn }; |
void (*addnumt[])() = { addq, addreal, addalg, ADDBF, addcplx, addmi, addlm, addgf2n, addgfpn, addgfs, addgfsn, adddalg }; |
void (*subnumt[])() = { subq, subreal, subalg, SUBBF, subcplx, submi, sublm, subgf2n, subgfpn, subgfs, subgfsn }; |
void (*subnumt[])() = { subq, subreal, subalg, SUBBF, subcplx, submi, sublm, subgf2n, subgfpn, subgfs, subgfsn, subdalg }; |
void (*mulnumt[])() = { mulq, mulreal, mulalg, MULBF, mulcplx, mulmi, mullm, mulgf2n, mulgfpn, mulgfs, mulgfsn }; |
void (*mulnumt[])() = { mulq, mulreal, mulalg, MULBF, mulcplx, mulmi, mullm, mulgf2n, mulgfpn, mulgfs, mulgfsn, muldalg }; |
void (*divnumt[])() = { divq, divreal, divalg, DIVBF, divcplx, divmi, divlm, divgf2n, divgfpn, divgfs, divgfsn }; |
void (*divnumt[])() = { divq, divreal, divalg, DIVBF, divcplx, divmi, divlm, divgf2n, divgfpn, divgfs, divgfsn, divdalg }; |
void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, PWRBF, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn, pwrgfs, pwrgfsn }; |
void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, PWRBF, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn, pwrgfs, pwrgfsn, pwrdalg }; |
void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, CHSGNBF, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn, chsgngfs, chsgngfsn }; |
void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, CHSGNBF, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn, chsgngfs, chsgngfsn, chsgndalg }; |
int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, CMPBF, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn, cmpgfs, cmpgfsn }; |
int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, CMPBF, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn, cmpgfs, cmpgfsn, cmpdalg }; |
#endif |
#endif |
|
|
double get_current_time(); |
double get_current_time(); |
Line 177 void nglob_init() { |
|
Line 180 void nglob_init() { |
|
|
|
/* moved to parse/main.c */ |
/* moved to parse/main.c */ |
#if 0 |
#if 0 |
#if PARI |
#if defined(PARI) |
risa_pari_init(); |
risa_pari_init(); |
#endif |
#endif |
srandom((int)get_current_time()); |
srandom((int)get_current_time()); |
Line 191 extern double gctime; |
|
Line 194 extern double gctime; |
|
double suspend_start; |
double suspend_start; |
double suspended_time=0; |
double suspended_time=0; |
|
|
void get_eg(p) |
void get_eg(struct oEGT *p) |
struct oEGT *p; |
|
{ |
{ |
p->exectime = get_clock() - gctime - suspended_time; p->gctime = gctime; |
p->exectime = get_clock() - gctime - suspended_time; p->gctime = gctime; |
} |
} |
|
|
void init_eg(eg) |
void init_eg(struct oEGT *eg) |
struct oEGT *eg; |
|
{ |
{ |
bzero((char *)eg,sizeof(struct oEGT)); |
bzero((char *)eg,sizeof(struct oEGT)); |
} |
} |
|
|
void add_eg(base,start,end) |
void add_eg(struct oEGT *base,struct oEGT *start,struct oEGT *end) |
struct oEGT *base,*start,*end; |
|
{ |
{ |
base->exectime += end->exectime - start->exectime; |
base->exectime += end->exectime - start->exectime; |
base->gctime += end->gctime - start->gctime; |
base->gctime += end->gctime - start->gctime; |
} |
} |
|
|
void print_eg(item,eg) |
void print_eg(char *item,struct oEGT *eg) |
char *item; |
|
struct oEGT *eg; |
|
{ |
{ |
printf("%s=(%.4g,%.4g)",item,eg->exectime,eg->gctime); |
printf("%s=(%.4g,%.4g)",item,eg->exectime,eg->gctime); |
} |
} |
|
|
void print_split_eg(start,end) |
void print_split_eg(struct oEGT *start,struct oEGT *end) |
struct oEGT *start,*end; |
|
{ |
{ |
struct oEGT base; |
struct oEGT base; |
|
|
Line 226 struct oEGT *start,*end; |
|
Line 223 struct oEGT *start,*end; |
|
printf("(%.4g,%.4g)",base.exectime,base.gctime); |
printf("(%.4g,%.4g)",base.exectime,base.gctime); |
} |
} |
|
|
void print_split_e(start,end) |
void print_split_e(struct oEGT *start,struct oEGT *end) |
struct oEGT *start,*end; |
|
{ |
{ |
struct oEGT base; |
struct oEGT base; |
|
|
Line 265 char *get_asir_distribution() { |
|
Line 261 char *get_asir_distribution() { |
|
} |
} |
|
|
|
|
void create_error(ERR *err,unsigned int serial,char *msg) |
void create_error(ERR *err,unsigned int serial,char *msg,LIST trace) |
{ |
{ |
int len; |
int len; |
USINT ui; |
USINT ui,notsupp; |
NODE n,n1; |
NODE n,n1; |
LIST list; |
LIST list; |
char *msg1; |
char *msg1; |
STRING errmsg; |
STRING errmsg; |
|
|
MKUSINT(ui,serial); |
MKUSINT(ui,serial); |
|
MKUSINT(notsupp,-1); |
len = strlen(msg); |
len = strlen(msg); |
msg1 = (char *)MALLOC(len+1); |
msg1 = (char *)MALLOC(len+1); |
strcpy(msg1,msg); |
strcpy(msg1,msg); |
MKSTR(errmsg,msg1); |
MKSTR(errmsg,msg1); |
MKNODE(n1,errmsg,0); MKNODE(n,ui,n1); MKLIST(list,n); |
if ( !trace ) |
|
MKLIST(trace,0); |
|
n = mknode(4,ui,notsupp,errmsg,trace); MKLIST(list,n); |
MKERR(*err,list); |
MKERR(*err,list); |
} |
} |
|
|
Line 303 int get_lprime(index) |
|
Line 302 int get_lprime(index) |
|
return lprime[index]; |
return lprime[index]; |
} |
} |
|
|
void create_new_lprimes(index) |
void create_new_lprimes(int index) |
int index; |
|
{ |
{ |
int count,p,i,j,d; |
int count,p,i,j,d; |
|
|