[BACK]Return to serv2.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / ox_math

Diff for /OpenXM/src/ox_math/Attic/serv2.c between version 1.6 and 1.18

version 1.6, 1999/11/06 21:39:37 version 1.18, 2000/12/03 15:19:23
Line 1 
Line 1 
 /* -*- mode: C; coding: euc-japan -*- */  /* -*- mode: C; coding: euc-japan -*- */
 /* $OpenXM: OpenXM/src/ox_math/serv2.c,v 1.5 1999/11/04 19:33:17 ohara Exp $ */  /* $OpenXM: OpenXM/src/ox_math/serv2.c,v 1.17 2000/11/28 20:16:03 ohara Exp $ */
   
 /* Open Mathematica サーバ */  /*
 /* ファイルディスクリプタ 3, 4 は open されていると仮定して動作する. */     Copyright (C) Katsuyoshi OHARA, 2000.
      Portions copyright 1999 Wolfram Research, Inc.
   
 /* MathLink との通信部分 */     You must see OpenXM/Copyright/Copyright.generic.
      The MathLink Library is licensed from Wolfram Research Inc..
      See OpenXM/Copyright/Copyright.mathlink for detail.
   */
   
   /*
      Remarks:
      file descripter 3 and 4 are already opened by the parent process.
   */
   
 #include <stdio.h>  #include <stdio.h>
 #include <stdlib.h>  #include <stdlib.h>
 #include <unistd.h>  #include <unistd.h>
 #include <gmp.h>  #include <gmp.h>
 #include <mathlink.h>  #include <mathlink.h>
 #include "ox.h"  #include <ox_toolkit.h>
 #include "parse.h"  
 #include "serv2.h"  #include "serv2.h"
   
 #define UNKNOWN_SM_COMMAND 50000  extern int flag_mlo_symbol;
 #define MATH_ERROR         50001  
   
 /* MLINK はポインタ型. */  /* MathLink independent */
 MLINK lp = NULL;  static cmo **stack = NULL;
   static int stack_size = 0;
   static int stack_ptr = 0;
   OXFILE *stack_oxfp = NULL;
   
 typedef cmo mlo;  #define DIFFERENCE_OF_STACK  1024
 typedef cmo_string mlo_string;  
 typedef cmo_zz mlo_zz;  
   
 /* cmo_list の派生クラス*/  static void stack_extend()
 typedef struct {  
         int tag;  
         int length;  
         cell head[1];  
         char *function;  
 } mlo_function;  
   
   
 mlo *receive_mlo_zz()  
 {  {
         char *s;      int newsize = stack_size + DIFFERENCE_OF_STACK;
         mlo  *m;      cmo **newstack = (cmo **)malloc(sizeof(cmo *)*newsize);
       if (stack != NULL) {
         fprintf(stderr, "--debug: MLO == MLTKINT.\n");          memcpy(newstack, stack, sizeof(cmo *)*stack_size);
         MLGetString(lp, &s);          free(stack);
         fprintf(stderr, "--debug: zz = %s.\n", s);  
         m = (mlo *)new_cmo_zz_set_string(s);  
         MLDisownString(lp, s);  
         return m;  
 }  
   
 mlo *receive_mlo_string()  
 {  
         char *s;  
         mlo  *m;  
         fprintf(stderr, "--debug: MLO == MLTKSTR.\n");  
         MLGetString(lp, &s);  
         fprintf(stderr, "--debug: string = \"%s\".\n", s);  
         m = (cmo *)new_cmo_string(s);  
         MLDisownString(lp, s);  
         return m;  
 }  
   
 cmo *receive_mlo_function()  
 {  
         char *s;  
         cmo *m;  
     cmo  *ob;  
     int  i,n;  
   
         fprintf(stderr, "--debug: MLO == MLTKFUNC.\n");  
         MLGetFunction(lp, &s, &n);  
         fprintf(stderr, "--debug: Function = \"%s\", # of args = %d\n", s, n);  
         m = new_cmo_list();  
         append_cmo_list((cmo_list *)m, new_cmo_string(s));  
   
         for (i=0; i<n; i++) {  
                 fprintf(stderr, "--debug: arg[%d]\n", i);  
                 fflush(stderr);  
                 ob = receive_mlo();  
                 append_cmo_list((cmo_list *)m, ob);  
         }  
   
         MLDisownString(lp, s);  
         return m;  
 }  
   
 cmo *receive_mlo_symbol()  
 {  
         cmo *ob;  
         char *s;  
   
         fprintf(stderr, "--debug: MLO == MLTKSYM.\n");  
         MLGetSymbol(lp, &s);  
         fprintf(stderr, "--debug: Symbol \"%s\".\n", s);  
   
         ob = new_cmo_indeterminate(new_cmo_string(s));  
   
         MLDisownString(lp, s);  
         return ob;  
 }  
   
 /* Mathematica を起動する. */  
 int MATH_init()  
 {  
     int argc = 2;  
     char *argv[] = {"-linkname", "math -mathlink"};  
   
     if(MLInitialize(NULL) == NULL  
            || (lp = MLOpen(argc, argv)) == NULL) {  
                 fprintf(stderr, "Mathematica Kernel not found.\n");  
                 exit(1);  
     }      }
         return 0;      stack_size = newsize;
       stack = newstack;
 }  }
   
 int MATH_exit()  
 {  
     /* quit Mathematica then close the link */  
     MLPutFunction(lp, "Exit", 0);  
     MLClose(lp);  
 }  
   
 cmo *MATH_getObject2()  
 {  
     /* skip any packets before the first ReturnPacket */  
     while (MLNextPacket(lp) != RETURNPKT) {  
         usleep(10);  
         MLNewPacket(lp);  
     }  
         return receive_mlo();  
 }  
   
 cmo *receive_mlo()  
 {  
     char *s;  
         int type;  
   
     switch(type = MLGetNext(lp)) {  
     case MLTKINT:  
                 return receive_mlo_zz();  
     case MLTKSTR:  
                 return receive_mlo_string();  
     case MLTKREAL:  
                 /* double はまだ... */  
         fprintf(stderr, "--debug: MLO == MLTKREAL.\n");  
         MLGetString(lp, &s);  
         return (cmo *)new_cmo_string(s);  
     case MLTKSYM:  
         return receive_mlo_symbol();  
     case MLTKFUNC:  
                 return receive_mlo_function();  
     case MLTKERR:  
         fprintf(stderr, "--debug: MLO == MLTKERR.\n");  
         return (cmo *)gen_error_object(MATH_ERROR);  
     default:  
         fprintf(stderr, "--debug: MLO(%d) is unknown.\n", type);  
         MLGetString(lp, &s);  
                 fprintf(stderr, "--debug: \"%s\"\n", s);  
         return (cmo *)new_cmo_string(s);  
     }  
 }  
   
 int send_mlo_int32(cmo *m)  
 {  
         MLPutInteger(lp, ((cmo_int32 *)m)->i);  
 }  
   
 int send_mlo_string(cmo *m)  
 {  
         char *s = ((cmo_string *)m)->s;  
         MLPutString(lp, s);  
         fprintf(stderr, "ox_math:: put %s.", s);  
 }  
   
 int send_mlo_zz(cmo *m)  
 {  
         char *s;  
         MLPutFunction(lp, "ToExpression", 1);  
         s = convert_cmo_to_string(m);  
         MLPutString(lp, s);  
         fprintf(stderr, "put %s.", s);  
 }  
   
 int send_mlo_list(cmo *c)  
 {  
         char *s;  
         cell *cp = ((cmo_list *)c)->head;  
         int len = length_cmo_list((cmo_list *)c);  
   
         fprintf(stderr, "ox_math:: put List with %d args.\n", len);  
         MLPutFunction(lp, "List", len);  
         while(cp->next != NULL) {  
                 send_mlo(cp->cmo);  
                 cp = cp->next;  
         }  
 }  
   
 int MATH_sendObject(cmo *m)  
 {  
         send_mlo(m);  
         MLEndPacket(lp);  
 }  
   
 int send_mlo(cmo *m)  
 {  
     char *s;  
     switch(m->tag) {  
     case CMO_INT32:  
                 send_mlo_int32(m);  
         break;  
     case CMO_STRING:  
                 send_mlo_string(m);  
         break;  
         case CMO_LIST:  
                 send_mlo_list(m);  
         break;  
     default:  
         MLPutFunction(lp, "ToExpression", 1);  
         s = convert_cmo_to_string(m);  
         MLPutString(lp, s);  
         fprintf(stderr, "put %s.", s);  
         break;  
     }  
 }  
   
 int MATH_evaluateStringByLocalParser(char *str)  
 {  
     MLPutFunction(lp, "ToExpression", 1);  
     MLPutString(lp, str);  
     MLEndPacket(lp);  
 }  
   
 int MATH_executeFunction(char *function, int argc, cmo *argv[])  
 {  
     int i;  
     MLPutFunction(lp, function, argc);  
     for (i=0; i<argc; i++) {  
         send_mlo(argv[i]);  
     }  
     MLEndPacket(lp);  
 }  
   
 /* MathLink 非依存部分 */  
   
 #define SIZE_OPERAND_STACK   2048  
   
 static cmo* Operand_Stack[SIZE_OPERAND_STACK];  
 static int Stack_Pointer = 0;  
   
 int initialize_stack()  
 {  
     Stack_Pointer = 0;  
 }  
   
 int push(cmo* m)  int push(cmo* m)
 {  {
 #if DEBUG  #if DEBUG
         symbol *symp;      symbol_t symp;
   
     if (m->tag == CMO_STRING) {      if (m->tag == CMO_STRING) {
         fprintf(stderr, "ox_math:: a cmo_string(%s) was pushed.\n", ((cmo_string *)m)->s);          fprintf(stderr, "ox_math:: a CMO_STRING(%s) was pushed.\n", ((cmo_string *)m)->s);
     }else {      }else {
                 symp = lookup_by_tag(m->tag);          symp = lookup_by_tag(m->tag);
                 fprintf(stderr, "ox_math:: a %s was pushed.\n", symp->key);          fprintf(stderr, "ox_math:: a %s was pushed.\n", symbol_get_key(symp));
         }      }
 #endif  #endif
     Operand_Stack[Stack_Pointer] = m;      if (stack_ptr >= stack_size) {
     Stack_Pointer++;          stack_extend();
     if (Stack_Pointer >= SIZE_OPERAND_STACK) {  
         fprintf(stderr, "stack over flow.\n");  
                 Stack_Pointer--;  
     }      }
       stack[stack_ptr] = m;
       stack_ptr++;
 }  }
   
 /* スタックが空のときは, (CMO_NULL) をかえす. */  /* if the stack is empty, then pop() returns (CMO_NULL). */
 cmo* pop()  cmo* pop()
 {  {
     if (Stack_Pointer > 0) {      if (stack_ptr > 0) {
         Stack_Pointer--;          return stack[--stack_ptr];
         return Operand_Stack[Stack_Pointer];  
     }      }
     return new_cmo_null();      return new_cmo_null();
 }  }
   
 void pops(int n)  void pops(int n)
 {  {
     Stack_Pointer -= n;      stack_ptr -= n;
     if (Stack_Pointer < 0) {      if (stack_ptr < 0) {
         Stack_Pointer = 0;          stack_ptr = 0;
     }      }
 }  }
   
 /* sm_XXX 関数群は、エラーのときは 0 以外の値を返し、呼び出し元で  void push_error(int errcode, cmo* pushback)
    エラーオブジェクトをセットする */  
 int sm_popCMO(int fd_write)  
 {  {
       return push((cmo *)make_error_object(errcode, pushback));
   }
   
   /*
      if error occurs, then a sm_*() function returns non-zero and
      an error obect is set by a function which calls sm_*().
   */
   int sm_popCMO(OXFILE* oxfp)
   {
     cmo* m = pop();      cmo* m = pop();
 #ifdef DEBUG  #ifdef DEBUG
         symbol *symp = lookup_by_tag(m->tag);      symbol_t symp = lookup_by_tag(m->tag);
       fprintf(stderr, "ox_math:: opecode = SM_popCMO. (%s)\n", symbol_get_key(symp));
     fprintf(stderr, "ox_math:: opecode = SM_popCMO. (%s)\n", symp->key);  
 #endif  #endif
   
     if (m != NULL) {      if (m != NULL) {
         send_ox_cmo(fd_write, m);          send_ox_cmo(oxfp, m);
         return 0;          return 0;
     }      }
     return SM_popCMO;      return SM_popCMO;
 }  }
   
 int sm_pops(int fd_write)  int sm_pops(OXFILE* oxfp)
 {  {
     cmo* m = pop();      cmo* m = pop();
     if (m != NULL && m->tag == CMO_INT32) {      if (m != NULL && m->tag == CMO_INT32) {
         pops(((cmo_int32 *)m)->i);          pops(((cmo_int32 *)m)->i);
         return 0;          return 0;
     }      }
     return UNKNOWN_SM_COMMAND;      return ERROR_ID_UNKNOWN_SM;
 }  }
   
 /* MathLink 依存部分 */  /* MathLink dependent */
 int sm_popString(int fd_write)  int sm_popString(OXFILE* oxfp)
 {  {
     char *s;      char *s;
     cmo *err;      cmo *err;
Line 327  int sm_popString(int fd_write)
Line 126  int sm_popString(int fd_write)
     fprintf(stderr, "ox_math:: opecode = SM_popString.\n");      fprintf(stderr, "ox_math:: opecode = SM_popString.\n");
 #endif  #endif
   
         m = pop();      m = pop();
         if (m->tag == CMO_STRING) {      if (m->tag == CMO_STRING) {
         send_ox_cmo(fd_write, m);          send_ox_cmo(oxfp, m);
         }else if ((s = convert_cmo_to_string(m)) != NULL) {      }else if ((s = new_string_set_cmo(m)) != NULL) {
         send_ox_cmo(fd_write, (cmo *)new_cmo_string(s));          send_ox_cmo(oxfp, (cmo *)new_cmo_string(s));
     }else {      }else {
                 err = new_cmo_error2(m);          err = make_error_object(SM_popString, m);
                 send_ox_cmo(fd_write, err);          send_ox_cmo(oxfp, err);
         }      }
         return 0;      return 0;
 }  }
   
 int local_execute(char *s)  int local_execute(char *s)
 {  {
         return 0;      if(*s == 'i') {
           switch(s[1]) {
           case '+':
               flag_mlo_symbol = FLAG_MLTKSYM_IS_STRING;
               break;
           case '-':
           case '=':
           default:
               flag_mlo_symbol = FLAG_MLTKSYM_IS_INDETERMINATE;
           }
       }
       return 0;
 }  }
   
 /* この関数はサーバに依存する. */  /* The following function is depend on an implementation of a server. */
 int sm_executeStringByLocalParser(int fd_write)  int sm_executeStringByLocalParser(OXFILE* oxfp)
 {  {
         symbol *symp;      symbol_t symp;
     cmo* m = pop();      cmo* m = pop();
         char *s = NULL;      char *s = NULL;
 #ifdef DEBUG  #ifdef DEBUG
     fprintf(stderr, "ox_math:: opecode = SM_executeStringByLocalParser.\n");      fprintf(stderr, "ox_math:: opecode = SM_executeStringByLocalParser.\n");
 #endif  #endif
   
     if (m->tag == CMO_STRING      if (m->tag == CMO_STRING
                 && strlen(s = ((cmo_string *)m)->s) != 0) {          && strlen(s = ((cmo_string *)m)->s) != 0) {
                 if (s[0] == ':') {          if (s[0] == ':') {
                         local_execute(s);              local_execute(++s);
                 }else {          }else {
                         /* for mathematica */              /* for mathematica */
                         /* mathematica に文字列を送って評価させる */              /* Sending the string `s' to mathematica for its evaluation. */
                         MATH_evaluateStringByLocalParser(s);              ml_evaluateStringByLocalParser(s);
                         push(MATH_getObject2());              ml_select();
                 }              push(receive_mlo());
                 return 0;          }
           return 0;
     }      }
 #ifdef DEBUG  #ifdef DEBUG
         if ((symp = lookup_by_tag(m->tag)) != NULL) {      symp = lookup_by_tag(m->tag);
                 fprintf(stderr, "ox_math:: error. the top of stack is %s.\n", symp->key);      fprintf(stderr, "ox_math:: error. the top of stack is %s.\n", symbol_get_key(symp));
         }else {  
                 fprintf(stderr, "ox_math:: error. the top of stack is unknown cmo. (%d)\n", m->tag);  
         }  
 #endif  #endif
     return SM_executeStringByLocalParser;      return SM_executeStringByLocalParser;
 }  }
   
 int sm_executeFunction(int fd_write)  int sm_executeFunction(OXFILE* oxfp)
 {  {
     int i, argc;      int i, argc;
     cmo **argv;      cmo **argv;
Line 393  int sm_executeFunction(int fd_write)
Line 201  int sm_executeFunction(int fd_write)
     }      }
   
     argc = ((cmo_int32 *)m)->i;      argc = ((cmo_int32 *)m)->i;
     argv = malloc(sizeof(cmo *)*argc);      argv = malloc(argc*sizeof(cmo *));
     for (i=0; i<argc; i++) {      for (i=0; i<argc; i++) {
         argv[i] = pop();          argv[i] = pop();
     }      }
     MATH_executeFunction(func, argc, argv);      ml_executeFunction(func, argc, argv);
     push(MATH_getObject2());      ml_select();
       push(receive_mlo());
     return 0;      return 0;
 }  }
   
 /* 平成11年10月13日 */  int sm_mathcap(OXFILE* oxfp)
 #define VERSION 0x11102700  
 #define ID_STRING  "ox_math server 1999/10/28 17:29:25"  
   
 int sm_mathcap(int fd_write)  
 {  {
     cmo* c = make_mathcap_object(VERSION, ID_STRING);      push((cmo *)oxf_cmo_mathcap(oxfp));
     push(c);  
     return 0;      return 0;
 }  }
   
 int receive_sm_command(int fd_read)  void sm_set_mathcap(OXFILE *oxfp)
 {  {
     return receive_int32(fd_read);      cmo_mathcap *m = (cmo_mathcap *)pop();
       if (m->tag == CMO_MATHCAP) {
           oxf_mathcap_update(oxfp, m);
       }else {
           push_error(-1, m);
           /* an error object must be pushed */
       }
 }  }
   
 int execute_sm_command(int fd_write, int code)  int receive_sm_command(OXFILE* oxfp)
 {  {
       return receive_int32(oxfp);
   }
   
   int execute_sm_command(OXFILE* oxfp, int code)
   {
     int err = 0;      int err = 0;
   #ifdef DEBUG
       symbol_t sp = lookup_by_tag(code);
       fprintf(stderr, "ox_math:: %s received.\n", symbol_get_key(sp));
   #endif
   
     switch(code) {      switch(code) {
     case SM_popCMO:      case SM_popCMO:
         err = sm_popCMO(fd_write);          err = sm_popCMO(oxfp);
         break;          break;
     case SM_popString:      case SM_popString:
         err = sm_popString(fd_write);          err = sm_popString(oxfp);
         break;          break;
     case SM_mathcap:      case SM_mathcap:
         err = sm_mathcap(fd_write);          err = sm_mathcap(oxfp);
         break;          break;
     case SM_pops:      case SM_pops:
         err = sm_pops(fd_write);          err = sm_pops(oxfp);
         break;          break;
     case SM_executeStringByLocalParser:      case SM_executeStringByLocalParser:
         err = sm_executeStringByLocalParser(fd_write);      case SM_executeStringByLocalParserInBatchMode:
           err = sm_executeStringByLocalParser(oxfp);
         break;          break;
     case SM_executeFunction:      case SM_executeFunction:
         err = sm_executeFunction(fd_write);          err = sm_executeFunction(oxfp);
         break;          break;
       case SM_shutdown:
           shutdown();
           break;
     case SM_setMathCap:      case SM_setMathCap:
         pop();  /* 無視する */          pop();  /* ignore */
         break;          break;
     default:      default:
         fprintf(stderr, "unknown command: %d.\n", code);          fprintf(stderr, "unknown command: %d.\n", code);
         err = UNKNOWN_SM_COMMAND;          err = ERROR_ID_UNKNOWN_SM;
     }      }
   
     if (err != 0) {      if (err != 0) {
         push((cmo *)gen_error_object(err));          push((cmo *)make_error_object(err, new_cmo_null()));
     }      }
 }  }
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.18
FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>