Return to init.c CVS log

Up to [local] / OpenXM_contrib2 / asir2000 / engine

Annotation of OpenXM_contrib2/asir2000/engine/init.c, Revision 1.1.1.1

1.1       noro        1: /* $OpenXM: OpenXM/src/asir99/engine/init.c,v 1.1.1.1 1999/11/10 08:12:26 noro Exp $ */
                      2: #include "ca.h"
                      3: #include "version.h"
                      4:
                      5: struct oN oUNIN,oTWON,oTHREEN,oFOURN,oFIVEN,oSIXN,oSEVENN,oEIGHTN;
                      6: struct oQ oUNIQ,oTWOQ,oTHREEQ,oFOURQ,oFIVEQ,oSIXQ,oSEVENQ,oEIGHTQ;
                      7:
                      8: struct oUP2 oONEUP2;
                      9:
                     10: N ONEN = &oUNIN;
                     11: Q ONE = &oUNIQ;
                     12: Q TWO = &oTWOQ;
                     13: Q THREE = &oTHREEQ;
                     14: Q FOUR = &oFOURQ;
                     15: Q FIVE = &oFIVEQ;
                     16: Q SIX = &oSIXQ;
                     17: Q SEVEN = &oSEVENQ;
                     18: Q EIGHT = &oEIGHTQ;
                     19:
                     20: LM THREE_LM,FOUR_LM,EIGHT_LM;
                     21:
                     22: struct oR oUNIR;
                     23: struct oMQ oUNIMQ;
                     24: struct oLM oUNILM;
                     25: struct oC oIU;
                     26: struct oUSINT oVOID;
                     27:
                     28: MQ ONEM = &oUNIMQ;
                     29: LM ONELM = &oUNILM;
                     30: UP2 ONEUP2 = &oONEUP2;
                     31: C IU = &oIU;
                     32: USINT VOIDobj = &oVOID;
                     33:
                     34: int bigfloat;
                     35: int nez = 1;
                     36: int current_mod = 0;
                     37:
                     38: #if PARI
                     39: #if defined(THINK_C)
                     40: int paristack = 1<<16;
                     41: #else
                     42: int paristack = 1<<20;
                     43: #endif
                     44:
                     45: void (*addnumt[])() = { addq, addreal, addalg, addbf, addcplx, addmi, addlm, addgf2n, addgfpn };
                     46: void (*subnumt[])() = { subq, subreal, subalg, subbf, subcplx, submi, sublm, subgf2n, subgfpn };
                     47: void (*mulnumt[])() = { mulq, mulreal, mulalg, mulbf, mulcplx, mulmi, mullm, mulgf2n, mulgfpn };
                     48: void (*divnumt[])() = { divq, divreal, divalg, divbf, divcplx, divmi, divlm, divgf2n, divgfpn };
                     49: void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, pwrbf, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn };
                     50: void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, chsgnbf, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn };
                     51: int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, cmpbf, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn };
                     52: #else
                     53: void (*addnumt[])() = { addq, addreal, addalg, 0, addcplx, addmi, addlm, addgf2n, addgfpn };
                     54: void (*subnumt[])() = { subq, subreal, subalg, 0, subcplx, submi, sublm, subgf2n, subgfpn };
                     55: void (*mulnumt[])() = { mulq, mulreal, mulalg, 0, mulcplx, mulmi, mullm, mulgf2n, mulgfpn };
                     56: void (*divnumt[])() = { divq, divreal, divalg, 0, divcplx, divmi, divlm, divgf2n,divgfpn };
                     57: void (*pwrnumt[])() = { pwrq, pwrreal, pwralg, 0, pwrcplx, pwrmi, pwrlm, pwrgf2n, pwrgfpn };
                     58: void (*chsgnnumt[])() = { chsgnq, chsgnreal, chsgnalg, 0, chsgncplx, chsgnmi, chsgnlm, chsgngf2n, chsgngfpn };
                     59: int (*cmpnumt[])() = { cmpq, cmpreal, cmpalg, 0, cmpcplx, cmpmi, cmplm, cmpgf2n, cmpgfpn };
                     60: #endif
                     61:
                     62: double get_current_time();
                     63:
                     64: void nglob_init() {
                     65:        PL(&oUNIN) = 1; BD(&oUNIN)[0] = 1;
                     66:        PL(&oTWON) = 1; BD(&oTWON)[0] = 2;
                     67:        PL(&oTHREEN) = 1; BD(&oTHREEN)[0] = 3;
                     68:        PL(&oFOURN) = 1; BD(&oFOURN)[0] = 4;
                     69:        PL(&oFIVEN) = 1; BD(&oFIVEN)[0] = 5;
                     70:        PL(&oSIXN) = 1; BD(&oSIXN)[0] = 6;
                     71:        PL(&oSEVENN) = 1; BD(&oSEVENN)[0] = 7;
                     72:        PL(&oEIGHTN) = 1; BD(&oEIGHTN)[0] = 8;
                     73:
                     74:        oONEUP2.w = 1; oONEUP2.b[0] = 1;
                     75:
                     76:        OID(&oUNIQ) = O_N; NID(&oUNIQ) = N_Q;
                     77:        SGN(&oUNIQ) = 1; NM(&oUNIQ) = &oUNIN; DN(&oUNIQ) = 0;
                     78:
                     79:        OID(&oTWOQ) = O_N; NID(&oTWOQ) = N_Q;
                     80:        SGN(&oTWOQ) = 1; NM(&oTWOQ) = &oTWON; DN(&oTWOQ) = 0;
                     81:
                     82:        OID(&oTHREEQ) = O_N; NID(&oTHREEQ) = N_Q;
                     83:        SGN(&oTHREEQ) = 1; NM(&oTHREEQ) = &oTHREEN; DN(&oTHREEQ) = 0;
                     84:
                     85:        OID(&oFOURQ) = O_N; NID(&oFOURQ) = N_Q;
                     86:        SGN(&oFOURQ) = 1; NM(&oFOURQ) = &oFOURN; DN(&oFOURQ) = 0;
                     87:
                     88:        OID(&oFIVEQ) = O_N; NID(&oFIVEQ) = N_Q;
                     89:        SGN(&oFIVEQ) = 1; NM(&oFIVEQ) = &oFIVEN; DN(&oFIVEQ) = 0;
                     90:
                     91:        OID(&oSIXQ) = O_N; NID(&oSIXQ) = N_Q;
                     92:        SGN(&oSIXQ) = 1; NM(&oSIXQ) = &oSIXN; DN(&oSIXQ) = 0;
                     93:
                     94:        OID(&oSEVENQ) = O_N; NID(&oSEVENQ) = N_Q;
                     95:        SGN(&oSEVENQ) = 1; NM(&oSEVENQ) = &oSEVENN; DN(&oSEVENQ) = 0;
                     96:
                     97:        OID(&oEIGHTQ) = O_N; NID(&oEIGHTQ) = N_Q;
                     98:        SGN(&oEIGHTQ) = 1; NM(&oEIGHTQ) = &oEIGHTN; DN(&oEIGHTQ) = 0;
                     99:
                    100:        OID(&oUNIR) = O_R; NM(&oUNIR) = (P)&oUNIQ; DN(&oUNIR) = (P)&oUNIQ; oUNIR.reduced = 1;
                    101:        OID(&oUNIMQ) = O_N; NID(&oUNIMQ) = N_M; CONT(&oUNIMQ) = 1;
                    102:        OID(&oUNILM) = O_N; NID(&oUNILM) = N_LM; BDY(&oUNILM) = ONEN;
                    103:        OID(&oIU) = O_N; NID(&oIU) = N_C; oIU.r = 0; oIU.i = (Num)ONE;
                    104:
                    105:        MKLM(&oTHREEN,THREE_LM);
                    106:        MKLM(&oFOURN,FOUR_LM);
                    107:        MKLM(&oEIGHTN,EIGHT_LM);
                    108:
                    109:        OID(&oVOID) = O_VOID;
                    110:
                    111: /* moved to parse/main.c */
                    112: #if 0
                    113: #if PARI
                    114:        risa_pari_init();
                    115: #endif
                    116:        srandom((int)get_current_time());
                    117: #endif
                    118:        init_up2_tab();
                    119: }
                    120:
                    121: extern double gctime;
                    122: double suspend_start;
                    123: double suspended_time=0;
                    124:
                    125: void get_eg(p)
                    126: struct oEGT *p;
                    127: {
                    128:        p->exectime = get_clock() - gctime - suspended_time; p->gctime = gctime;
                    129: }
                    130:
                    131: void init_eg(eg)
                    132: struct oEGT *eg;
                    133: {
                    134:        bzero((char *)eg,sizeof(struct oEGT));
                    135: }
                    136:
                    137: void add_eg(base,start,end)
                    138: struct oEGT *base,*start,*end;
                    139: {
                    140:        base->exectime += end->exectime - start->exectime;
                    141:        base->gctime += end->gctime - start->gctime;
                    142: }
                    143:
                    144: void print_eg(item,eg)
                    145: char *item;
                    146: struct oEGT *eg;
                    147: {
                    148:        printf("%s=(%.4g,%.4g)",item,eg->exectime,eg->gctime);
                    149: }
                    150:
                    151: void print_split_eg(start,end)
                    152: struct oEGT *start,*end;
                    153: {
                    154:        struct oEGT base;
                    155:
                    156:        init_eg(&base); add_eg(&base,start,end);
                    157:        printf("(%.4g,%.4g)",base.exectime,base.gctime);
                    158: }
                    159:
                    160: void print_split_e(start,end)
                    161: struct oEGT *start,*end;
                    162: {
                    163:        struct oEGT base;
                    164:
                    165:        init_eg(&base); add_eg(&base,start,end);
                    166:        printf("(%.4g)",base.exectime);
                    167: }
                    168:
                    169: void suspend_timer() {
                    170:        suspend_start = get_clock();
                    171: }
                    172:
                    173: void resume_timer() {
                    174:        suspended_time += get_clock()-suspend_start;
                    175: }
                    176:
                    177: extern int lm_lazy, up_lazy;
                    178: extern GC_dont_gc;
                    179:
                    180: void reset_engine() {
                    181:        lm_lazy = 0; up_lazy = 0;
                    182:        GC_dont_gc = 0;
                    183: }
                    184:
                    185: unsigned int get_asir_version() {
                    186:        return ASIR_VERSION;
                    187: }