![[BACK]](/icons/cvsweb/back.gif){kind=icon}
Return to ox.c CVS log ![[TXT]](/icons/cvsweb/text.gif){kind=icon}
![[DIR]](/icons/cvsweb/dir.gif){kind=icon}
Up to [local] / OpenXM_contrib2 / asir2000 / io|
Return to ox.c CVS log | Up to [local] / OpenXM_contrib2 / asir2000 / io |
1.1 ! noro 1: /* $OpenXM: OpenXM/src/asir99/io/ox.c,v 1.7 1999/11/19 04:15:49 noro Exp $ */
! 2: #include "ca.h"
! 3: #include "parse.h"
! 4: #include "wsio.h"
! 5: #include "ox.h"
! 6:
! 7: #define ISIZ sizeof(int)
! 8:
! 9: void ox_flush_stream(),ox_write_int(),ox_write_cmo();
! 10: void ox_read_int(),ox_read_cmo();
! 11: void mclist_to_mc();
! 12: void ox_read_local();
! 13:
! 14: extern Obj VOIDobj;
! 15:
! 16: extern int ox_need_conv;
! 17: extern int ox_usr1_sent, ox_int_received, critical_when_signal;
! 18: unsigned int ox_serial;
! 19: int ox_flushing;
! 20: int ox_batch;
! 21: int ox_check=1;
! 22: jmp_buf ox_env;
! 23:
! 24: MATHCAP my_mathcap;
! 25:
! 26: struct oxcap {
! 27: unsigned int ox;
! 28: int ncap;
! 29: int *cap;
! 30: };
! 31:
! 32: struct mathcap {
! 33: LIST mc;
! 34: unsigned int version;
! 35: char *servername;
! 36: int nsmcap;
! 37: unsigned int *smcap;
! 38: int noxcap;
! 39: struct oxcap *oxcap;
! 40: };
! 41:
! 42: struct oxcap *my_oxcap;
! 43:
! 44: static struct mathcap my_mc;
! 45: static struct mathcap *remote_mc;
! 46: static int remote_mc_len;
! 47:
! 48: void ox_resetenv(s)
! 49: char *s;
! 50: {
! 51: fprintf(stderr,"%s\n",s);
! 52: longjmp(ox_env,1);
! 53: }
! 54:
! 55: static int available_cmo[] = {
! 56: CMO_NULL, CMO_INT32, CMO_DATUM, CMO_STRING, CMO_MATHCAP,
! 57: CMO_ERROR, CMO_ERROR2, CMO_LIST, CMO_MONOMIAL32,
! 58: CMO_ZZ, CMO_QQ, CMO_ZERO,
! 59: CMO_DMS_GENERIC, CMO_DMS_OF_N_VARIABLES,
! 60: CMO_RING_BY_NAME, CMO_DISTRIBUTED_POLYNOMIAL,
! 61: CMO_RECURSIVE_POLYNOMIAL, CMO_UNIVARIATE_POLYNOMIAL,
! 62: CMO_INDETERMINATE,
! 63: 0
! 64: };
! 65:
! 66: static int available_sm[] = {
! 67: SM_dupErrors, SM_getsp, SM_popSerializedLocalObject,
! 68: SM_popCMO, SM_popString, SM_setName,
! 69: SM_evalName, SM_executeStringByLocalParser,
! 70: SM_executeStringByLocalParserInBatchMode,
! 71: SM_executeFunction, SM_shutdown, SM_pops,
! 72: SM_mathcap, SM_setMathcap, SM_nop,
! 73: SM_beginBlock, SM_endBlock,
! 74: 0
! 75: };
! 76:
! 77: /*
! 78: mathcap =
! 79: [
! 80: version list,
! 81: SMlist,
! 82: [
! 83: [OX tag,CMO tag list],
! 84: [OX tag,CMO tag list],
! 85: ...
! 86: ]
! 87: ]
! 88: */
! 89:
! 90: void create_my_mathcap(char *system)
! 91: {
! 92: NODE n,n0;
! 93: int i,k;
! 94: STRING str;
! 95: LIST sname,smlist,oxlist,cmolist,asirlist,oxtag,oxasir,r;
! 96: USINT tag,t,t1;
! 97:
! 98: if ( my_mathcap )
! 99: return;
! 100: /* version */
! 101: MKSTR(str,system);
! 102: MKUSINT(t,OX_VERSION);
! 103: n0 = mknode(2,t,str); MKLIST(sname,n0);
! 104:
! 105: /* cmo tag */
! 106: for ( n0 = 0, i = 0; k = available_sm[i]; i++ ) {
! 107: NEXTNODE(n0,n); MKUSINT(t,k); BDY(n) = (pointer)t;
! 108: }
! 109: NEXT(n) = 0; MKLIST(smlist,n0);
! 110:
! 111: /* creation of [OX_DATA,CMO list] */
! 112: /* ox tag */
! 113: MKUSINT(tag,OX_DATA);
! 114: /* cmo tag */
! 115: for ( n0 = 0, i = 0; k = available_cmo[i]; i++ ) {
! 116: NEXTNODE(n0,n); MKUSINT(t,k); BDY(n) = (pointer)t;
! 117: }
! 118: NEXT(n) = 0; MKLIST(cmolist,n0);
! 119: /* [ox tag, cmo list] */
! 120: n0 = mknode(2,tag,cmolist);
! 121: MKLIST(oxlist,n0);
! 122:
! 123: /* creation of [OX_LOCAL_OBJECT_ASIR,ASIR tag] */
! 124: /* ox tag */
! 125: MKUSINT(tag,OX_LOCAL_OBJECT_ASIR);
! 126: /* local tag */
! 127: MKUSINT(t,ASIR_VL);
! 128: MKUSINT(t1,ASIR_OBJ);
! 129: n0 = mknode(2,t,t1); MKLIST(cmolist,n0);
! 130: /* [ox tag, local list] */
! 131: n0 = mknode(2,tag,cmolist);
! 132: MKLIST(asirlist,n0);
! 133:
! 134: /* [oxlist,asirlist] */
! 135: n0 = mknode(2,oxlist,asirlist); MKLIST(oxasir,n0);
! 136:
! 137: /* [version,sm,oxasir] */
! 138: n0 = mknode(3,sname,smlist,oxasir); MKLIST(r,n0);
! 139:
! 140: MKMATHCAP(my_mathcap,r);
! 141: mclist_to_mc(r,&my_mc);
! 142: my_oxcap = my_mc.oxcap;
! 143: }
! 144:
! 145: void store_remote_mathcap(int s,MATHCAP mc)
! 146: {
! 147: if ( !remote_mc ) {
! 148: remote_mc_len = 16;
! 149: remote_mc = (struct mathcap *)
! 150: CALLOC(remote_mc_len,sizeof(struct mathcap));
! 151: }
! 152: if ( s >= remote_mc_len ) {
! 153: remote_mc_len *= 2;
! 154: remote_mc = (struct mathcap *)REALLOC(remote_mc,
! 155: remote_mc_len*sizeof(struct mathcap));
! 156: }
! 157: mclist_to_mc(BDY(mc),&remote_mc[s]);
! 158: }
! 159:
! 160: /*
! 161: mathcap =
! 162: [
! 163: version list,
! 164: SMlist,
! 165: [
! 166: [OX tag,CMO tag list],
! 167: [OX tag,CMO tag list],
! 168: ...
! 169: ]
! 170: ]
! 171:
! 172: ===>
! 173:
! 174: mathcap
! 175: | version | &servername | nsmcap | &smcap | noxcap | &oxcap |
! 176: smcap
! 177: | SM_xxx | SM_yyy | ... |
! 178: oxcap
! 179: | oxcap[0] | oxcap[1] | ... |
! 180: oxcap[i]
! 181: | ox | ncap | &cap |
! 182: cap
! 183: | CMO_xxx | CMO_yyy | ... |
! 184: */
! 185:
! 186: void mclist_to_mc(LIST mclist,struct mathcap *mc)
! 187: {
! 188: int id,l,i,j;
! 189: NODE n,t,oxcmo,ox,cap;
! 190: int *ptr;
! 191:
! 192: /*
! 193: [
! 194: [ version,servername ]
! 195: [sm1,sm2,...],
! 196: [
! 197: [o1,[n11,n12,...]],
! 198: [o2,[n21,n22,...]],
! 199: ...
! 200: ]
! 201: ]
! 202: */
! 203: n = BDY(mclist);
! 204: mc->mc = mclist;
! 205: mc->version = BDY((USINT)BDY(BDY((LIST)BDY(n))));
! 206: mc->servername = BDY((STRING)BDY(NEXT(BDY((LIST)BDY(n)))));
! 207:
! 208: /* smcap */
! 209: n = NEXT(n);
! 210: t = BDY((LIST)BDY(n));
! 211: mc->nsmcap = length(t);
! 212: mc->smcap = (int *)MALLOC_ATOMIC(mc->nsmcap*sizeof(int));
! 213: for ( j = 0, ptr = mc->smcap; j < mc->nsmcap; j++, t = NEXT(t) )
! 214: ptr[j] = BDY((USINT)BDY(t));
! 215:
! 216: n = NEXT(n);
! 217: n = BDY((LIST)BDY(n));
! 218: /* n -> BDY([[OX1,CMOlist1], [OX2,CMOlist2], ...]) */
! 219: mc->noxcap = length(n);
! 220: mc->oxcap = (struct oxcap *)MALLOC(mc->noxcap*sizeof(struct oxcap));
! 221: for ( j = 0; j < mc->noxcap; j++, n = NEXT(n) ) {
! 222: oxcmo = BDY((LIST)BDY(n));
! 223: /* oxcmo = BDY([OXj,CMOlistj]) */
! 224: mc->oxcap[j].ox = BDY((USINT)BDY(oxcmo));
! 225: cap = BDY((LIST)BDY(NEXT(oxcmo)));
! 226: /* cap ->BDY(CMOlistj) */
! 227: l = length(cap);
! 228: mc->oxcap[j].ncap = l;
! 229: mc->oxcap[j].cap = (unsigned int *)CALLOC(l+1,sizeof(unsigned int));
! 230: for ( t = cap, ptr = mc->oxcap[j].cap, i = 0; i < l; t = NEXT(t), i++ )
! 231: ptr[i] = BDY((USINT)BDY(t));
! 232: }
! 233: }
! 234:
! 235: int check_sm_by_mc(s,smtag)
! 236: int s;
! 237: unsigned int smtag;
! 238: {
! 239: struct mathcap *rmc;
! 240: int nsmcap,i;
! 241: unsigned int *smcap;
! 242:
! 243: /* XXX : return 1 if remote_mc is not available. */
! 244: if ( !remote_mc )
! 245: return 1;
! 246: rmc = &remote_mc[s];
! 247: nsmcap = rmc->nsmcap;
! 248: smcap = rmc->smcap;
! 249: if ( !smcap )
! 250: return 1;
! 251: for ( i = 0; i < nsmcap; i++ )
! 252: if ( smcap[i] == smtag )
! 253: break;
! 254: if ( i == nsmcap )
! 255: return 0;
! 256: else
! 257: return 1;
! 258: }
! 259:
! 260: int check_by_mc(s,oxtag,cmotag)
! 261: int s;
! 262: unsigned int oxtag,cmotag;
! 263: {
! 264: struct mathcap *rmc;
! 265: int noxcap,ncap,i,j;
! 266: struct oxcap *oxcap;
! 267: unsigned int *cap;
! 268:
! 269: /* XXX : return 1 if remote_mc is not available. */
! 270: if ( !remote_mc )
! 271: return 1;
! 272: rmc = &remote_mc[s];
! 273: noxcap = rmc->noxcap;
! 274: oxcap = rmc->oxcap;
! 275: if ( !oxcap )
! 276: return 1;
! 277: for ( i = 0; i < noxcap; i++ )
! 278: if ( oxcap[i].ox == oxtag )
! 279: break;
! 280: if ( i == noxcap )
! 281: return 0;
! 282: ncap = oxcap[i].ncap;
! 283: cap = oxcap[i].cap;
! 284: for ( j = 0; j < ncap; j++ )
! 285: if ( cap[j] == cmotag )
! 286: break;
! 287: if ( j == ncap )
! 288: return 0;
! 289: else
! 290: return 1;
! 291: }
! 292:
! 293: void begin_critical() {
! 294: critical_when_signal = 1;
! 295: }
! 296:
! 297: void end_critical() {
! 298: critical_when_signal = 0;
! 299: if ( ox_usr1_sent ) {
! 300: ox_usr1_sent = 0; ox_usr1_handler();
! 301: }
! 302: if ( ox_int_received ) {
! 303: ox_int_received = 0; int_handler(SIGINT);
! 304: }
! 305: }
! 306:
! 307: void ox_usr1_handler(sig)
! 308: int sig;
! 309: {
! 310: extern jmp_buf env;
! 311: unsigned int cmd;
! 312:
! 313: #if !defined(VISUAL)
! 314: signal(SIGUSR1,ox_usr1_handler);
! 315: #endif
! 316: if ( critical_when_signal ) {
! 317: fprintf(stderr,"usr1 : critical\n");
! 318: ox_usr1_sent = 1;
! 319: } else {
! 320: ox_flushing = 1;
! 321: ox_resetenv("usr1 : return to toplevel by SIGUSR1");
! 322: }
! 323: }
! 324:
! 325: void clear_readbuffer()
! 326: {
! 327: char c;
! 328: fd_set r,w,e;
! 329: struct timeval interval;
! 330: int n,sock;
! 331:
! 332: #if defined(linux)
! 333: iofp[0].in->_IO_read_ptr = iofp[0].in->_IO_read_end;
! 334: #elif defined(__FreeBSD__)
! 335: fpurge(iofp[0].in);
! 336: #endif
! 337: /*
! 338: sock = fileno(iofp[0].in);
! 339: interval.tv_sec = (int)0;
! 340: interval.tv_usec = (int)0;
! 341:
! 342: FD_ZERO(&r); FD_ZERO(&w); FD_ZERO(&e);
! 343: FD_SET(sock,&r);
! 344: while ( 1 ) {
! 345: n = select(FD_SETSIZE,&r,&w,&e,&interval);
! 346: if ( !n )
! 347: break;
! 348: read(sock,&c,1);
! 349: }
! 350: */
! 351: }
! 352:
! 353: #if MPI
! 354: int ox_data_is_available(int s)
! 355: {
! 356: return 1;
! 357: }
! 358:
! 359: void wait_for_data(int s)
! 360: {
! 361: return;
! 362: }
! 363: #else
! 364: int ox_data_is_available(int s)
! 365: {
! 366: return FP_DATA_IS_AVAILABLE(iofp[s].in);
! 367: }
! 368:
! 369: void wait_for_data(int s)
! 370: {
! 371: fd_set r;
! 372: int sock;
! 373:
! 374: if ( !FP_DATA_IS_AVAILABLE(iofp[s].in) ) {
! 375: #if defined(VISUAL)
! 376: sock = iofp[s].in->fildes;
! 377: FD_ZERO(&r);
! 378: FD_SET(sock,&r);
! 379: select(0,&r,NULL,NULL,NULL);
! 380: #else
! 381: sock = fileno(iofp[s].in);
! 382: FD_ZERO(&r);
! 383: FD_SET(sock,&r);
! 384: select(FD_SETSIZE,&r,NULL,NULL,NULL);
! 385: #endif
! 386: }
! 387: }
! 388: #endif
! 389:
! 390: void ox_send_data(int s,pointer p)
! 391: {
! 392: if ( ox_check && !ox_check_cmo(s,(Obj)p) )
! 393: error("ox_send_data : Mathcap violation");
! 394: begin_critical();
! 395: ox_write_int(s,OX_DATA);
! 396: ox_write_int(s,ox_serial++);
! 397: ox_write_cmo(s,p);
! 398: ox_flush_stream(s);
! 399: end_critical();
! 400: }
! 401:
! 402: void ox_send_cmd(int s,int id)
! 403: {
! 404: if ( ox_check && !check_sm_by_mc(s,id) )
! 405: error("ox_send_cmd : Mathcap violation");
! 406: begin_critical();
! 407: ox_write_int(s,OX_COMMAND);
! 408: ox_write_int(s,ox_serial++);
! 409: ox_write_int(s,id);
! 410: ox_flush_stream(s);
! 411: end_critical();
! 412: }
! 413:
! 414: void ox_send_sync(int s)
! 415: {
! 416: begin_critical();
! 417: ox_write_int(s,OX_SYNC_BALL);
! 418: ox_write_int(s,ox_serial++);
! 419: ox_flush_stream(s);
! 420: end_critical();
! 421: }
! 422:
! 423: void ox_send_local_data(int s,Obj p)
! 424: {
! 425: begin_critical();
! 426: ox_write_int(s,OX_LOCAL_OBJECT_ASIR);
! 427: ox_write_int(s,ox_serial++);
! 428: ox_write_int(s,ASIR_OBJ);
! 429: saveobj(iofp[s].out,p);
! 430: ox_flush_stream(s);
! 431: end_critical();
! 432: }
! 433:
! 434: void ox_send_local_ring(int s,VL vl)
! 435: {
! 436: begin_critical();
! 437: ox_write_int(s,OX_LOCAL_OBJECT_ASIR);
! 438: ox_write_int(s,ox_serial++);
! 439: ox_write_int(s,ASIR_VL);
! 440: savevl(iofp[s].out,vl);
! 441: ox_flush_stream(s);
! 442: end_critical();
! 443: }
! 444:
! 445: unsigned int ox_recv(int s, int *id, pointer *p)
! 446: {
! 447: unsigned int cmd,serial;
! 448: USINT ui;
! 449:
! 450: wait_for_data(s);
! 451: begin_critical();
! 452: ox_read_int(s,id);
! 453: ox_read_int(s,&serial);
! 454: switch ( *id ) {
! 455: case OX_COMMAND:
! 456: ox_read_int(s,&cmd);
! 457: MKUSINT(ui,cmd);
! 458: *p = (pointer)ui;
! 459: break;
! 460: case OX_DATA:
! 461: ox_read_cmo(s,p);
! 462: break;
! 463: case OX_LOCAL_OBJECT_ASIR:
! 464: ox_read_local(s,p);
! 465: break;
! 466: default:
! 467: *p = 0;
! 468: break;
! 469: }
! 470: end_critical();
! 471: return serial;
! 472: }
! 473:
! 474: void ox_get_result(s,rp)
! 475: int s;
! 476: Obj *rp;
! 477: {
! 478: int id;
! 479: Obj obj,r;
! 480: int level;
! 481:
! 482: level = 0;
! 483: r = 0;
! 484: do {
! 485: ox_recv(s,&id,(pointer *)&obj);
! 486: if ( id == OX_COMMAND ) {
! 487: switch ( ((USINT)obj)->body ) {
! 488: case SM_beginBlock:
! 489: level++;
! 490: break;
! 491: case SM_endBlock:
! 492: level--;
! 493: }
! 494: } else
! 495: r = obj;
! 496: } while ( level );
! 497: *rp = r;
! 498: }
! 499:
! 500: void ox_read_int(int s, int *n)
! 501: {
! 502: ox_need_conv = iofp[s].conv;
! 503: read_int(iofp[s].in,n);
! 504: }
! 505:
! 506: void ox_read_cmo(int s, Obj *rp)
! 507: {
! 508: ox_need_conv = iofp[s].conv;
! 509: read_cmo(iofp[s].in,rp);
! 510: }
! 511:
! 512: void ox_read_local(int s, Obj *rp)
! 513: {
! 514: int id;
! 515:
! 516: ox_need_conv = iofp[s].conv;
! 517: read_int(iofp[s].in,&id);
! 518: switch ( id ) {
! 519: case ASIR_VL:
! 520: loadvl(iofp[s].in);
! 521: *rp = VOIDobj;
! 522: break;
! 523: case ASIR_OBJ:
! 524: loadobj(iofp[s].in,rp);
! 525: break;
! 526: default:
! 527: error("ox_read_local : unsupported id");
! 528: break;
! 529: }
! 530: }
! 531:
! 532: void ox_write_int(int s, int n)
! 533: {
! 534: ox_need_conv = iofp[s].conv;
! 535: write_int(iofp[s].out,&n);
! 536: }
! 537:
! 538: void ox_write_cmo(int s, Obj obj)
! 539: {
! 540: ox_need_conv = iofp[s].conv;
! 541: write_cmo(iofp[s].out,obj);
! 542: }
! 543:
! 544: int ox_check_cmo(int s, Obj obj)
! 545: {
! 546: NODE m;
! 547:
! 548: if ( !obj )
! 549: return 1;
! 550: switch ( OID(obj) ) {
! 551: case O_MATHCAP: case O_STR: case O_ERR: case O_USINT: case O_VOID:
! 552: return 1;
! 553: case O_P:
! 554: if ( !check_by_mc(s,OX_DATA,CMO_RECURSIVE_POLYNOMIAL) )
! 555: return 0;
! 556: else
! 557: return ox_check_cmo_p(s,(P)obj);
! 558: case O_R:
! 559: if ( !check_by_mc(s,OX_DATA,CMO_RATIONAL) )
! 560: return 0;
! 561: else if ( !check_by_mc(s,OX_DATA,CMO_RECURSIVE_POLYNOMIAL) )
! 562: return 0;
! 563: else
! 564: return ox_check_cmo_p(s,NM((R)obj)) && ox_check_cmo_p(s,DN((R)obj));
! 565: case O_DP:
! 566: return ox_check_cmo_dp(s,(DP)obj);
! 567: case O_N:
! 568: if ( NID((Num)obj) == N_Q ) {
! 569: if ( INT((Q)obj) )
! 570: return check_by_mc(s,OX_DATA,CMO_ZZ);
! 571: else
! 572: return check_by_mc(s,OX_DATA,CMO_QQ);
! 573: } else
! 574: return 0;
! 575: case O_LIST:
! 576: for ( m = BDY((LIST)obj); m; m = NEXT(m) )
! 577: if ( !ox_check_cmo(s,(BDY(m))) )
! 578: return 0;
! 579: return 1;
! 580: default:
! 581: return 0;
! 582: }
! 583: }
! 584:
! 585: void ox_get_serverinfo(int s, LIST *rp)
! 586: {
! 587: if ( remote_mc )
! 588: *rp = remote_mc[s].mc;
! 589: else {
! 590: MKLIST(*rp,0);
! 591: }
! 592: }
! 593:
! 594: int ox_check_cmo_p(int s, P p)
! 595: {
! 596: DCP dc;
! 597:
! 598: if ( NUM(p) )
! 599: return ox_check_cmo(s,(Obj)p);
! 600: else {
! 601: for ( dc = DC(p); dc; dc = NEXT(dc) )
! 602: if ( !ox_check_cmo_p(s,COEF(dc)) )
! 603: return 0;
! 604: return 1;
! 605: }
! 606: }
! 607:
! 608: int ox_check_cmo_dp(int s, DP p)
! 609: {
! 610: MP m;
! 611:
! 612: for ( m = BDY(p); m; m = NEXT(m) )
! 613: if ( !ox_check_cmo(s,(Obj)m->c) )
! 614: return 0;
! 615: return 1;
! 616: }
! 617:
! 618: void ox_flush_stream(s)
! 619: int s;
! 620: {
! 621: if ( ox_batch )
! 622: return;
! 623: #if defined(VISUAL)
! 624: if ( _fileno(&iofp[s].out->fp) < 0 )
! 625: cflush(iofp[s].out);
! 626: else
! 627: #elif MPI
! 628: if ( (char)fileno(&iofp[s].out->fp) < 0 )
! 629: cflush(iofp[s].out);
! 630: else
! 631: #endif
! 632: fflush(iofp[s].out);
! 633: }
! 634:
! 635: void ox_flush_stream_force(s)
! 636: int s;
! 637: {
! 638: #if defined(VISUAL)
! 639: if ( _fileno(&iofp[s].out->fp) < 0 )
! 640: cflush(iofp[s].out);
! 641: else
! 642: #elif MPI
! 643: if ( (char)fileno(&iofp[s].out->fp) < 0 )
! 644: cflush(iofp[s].out);
! 645: else
! 646: #endif
! 647: fflush(iofp[s].out);
! 648: }