version 1.5, 1999/11/04 19:33:17 |
version 1.12, 2000/01/05 06:09:11 |
|
|
/* -*- mode: C; coding: euc-japan -*- */ |
/* -*- mode: C; coding: euc-japan -*- */ |
/* $OpenXM: OpenXM/src/ox_math/serv2.c,v 1.4 1999/11/04 03:05:51 ohara Exp $ */ |
/* $OpenXM: OpenXM/src/ox_math/serv2.c,v 1.11 1999/12/14 09:31:56 ohara Exp $ */ |
|
|
/* Open Mathematica サーバ */ |
/* Open Mathematica サーバ */ |
/* ファイルディスクリプタ 3, 4 は open されていると仮定して動作する. */ |
/* ファイルディスクリプタ 3, 4 は open されていると仮定して動作する. */ |
|
|
/* MathLink との通信部分 */ |
|
|
|
#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.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 非依存部分 */ |
MLINK lp = NULL; |
#define INIT_S_SIZE 2048 |
|
#define EXT_S_SIZE 2048 |
|
|
|
static int stack_size = 0; |
|
static int stack_pointer = 0; |
|
static cmo **stack = NULL; |
|
|
typedef cmo mlo; |
int initialize_stack() |
typedef cmo_string mlo_string; |
|
typedef cmo_zz mlo_zz; |
|
|
|
/* cmo_list の派生クラス*/ |
|
typedef struct { |
|
int tag; |
|
int length; |
|
cell head[1]; |
|
char *function; |
|
} mlo_function; |
|
|
|
|
|
mlo *receive_mlo_zz() |
|
{ |
{ |
char *s; |
stack_pointer = 0; |
mlo *m; |
stack_size = INIT_S_SIZE; |
|
stack = malloc(stack_size*sizeof(cmo*)); |
fprintf(stderr, "--debug: MLO == MLTKINT.\n"); |
|
MLGetString(lp, &s); |
|
fprintf(stderr, "--debug: zz = %s.\n", s); |
|
m = (mlo *)new_cmo_zz_set_string(s); |
|
MLDisownString(lp, s); |
|
return m; |
|
} |
} |
|
|
mlo *receive_mlo_string() |
static int extend_stack() |
{ |
{ |
char *s; |
int size2 = stack_size + EXT_S_SIZE; |
mlo *m; |
cmo **stack2 = malloc(size2*sizeof(cmo*)); |
fprintf(stderr, "--debug: MLO == MLTKSTR.\n"); |
memcpy(stack2, stack, stack_size*sizeof(cmo *)); |
MLGetString(lp, &s); |
free(stack); |
fprintf(stderr, "--debug: string = \"%s\".\n", s); |
stack = stack2; |
m = (cmo *)new_cmo_string(s); |
stack_size = size2; |
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(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(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; |
|
} |
|
|
|
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 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 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 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; |
|
|
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 { |
fprintf(stderr, "ox_math:: a cmo(%d) was pushed.\n", m->tag); |
symp = lookup_by_tag(m->tag); |
} |
fprintf(stderr, "ox_math:: a %s was pushed.\n", symp->key); |
|
} |
#endif |
#endif |
Operand_Stack[Stack_Pointer] = m; |
stack[stack_pointer] = m; |
Stack_Pointer++; |
stack_pointer++; |
if (Stack_Pointer >= SIZE_OPERAND_STACK) { |
if (stack_pointer >= stack_size) { |
fprintf(stderr, "stack over flow.\n"); |
extend_stack(); |
exit(1); /* 手抜き */ |
|
} |
} |
} |
} |
|
|
/* スタックが空のときは, (CMO_NULL) をかえす. */ |
/* スタックが空のときは, (CMO_NULL) をかえす. */ |
cmo* pop() |
cmo* pop() |
{ |
{ |
if (Stack_Pointer > 0) { |
if (stack_pointer > 0) { |
Stack_Pointer--; |
stack_pointer--; |
return Operand_Stack[Stack_Pointer]; |
return stack[stack_pointer]; |
} |
} |
return new_cmo_null(); |
return new_cmo_null(); |
} |
} |
|
|
void pops(int n) |
void pops(int n) |
{ |
{ |
Stack_Pointer -= n; |
stack_pointer -= n; |
if (Stack_Pointer < 0) { |
if (stack_pointer < 0) { |
Stack_Pointer = 0; |
stack_pointer = 0; |
} |
} |
} |
} |
|
|
Line 293 void pops(int n) |
|
|
int sm_popCMO(int fd_write) |
int sm_popCMO(int fd_write) |
{ |
{ |
cmo* m = pop(); |
cmo* m = pop(); |
|
#ifdef DEBUG |
|
symbol *symp = lookup_by_tag(m->tag); |
|
fprintf(stderr, "ox_math:: opecode = SM_popCMO. (%s)\n", symp->key); |
|
#endif |
|
|
fprintf(stderr, "ox_math:: opecode = SM_popCMO. (tag = %d)\n", m->tag); |
|
if (m != NULL) { |
if (m != NULL) { |
send_ox_cmo(fd_write, m); |
send_ox_cmo(fd_write, m); |
return 0; |
return 0; |
Line 309 int sm_pops(int fd_write) |
|
Line 101 int sm_pops(int fd_write) |
|
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 依存部分 */ |
int sm_popString(int fd_write) |
int sm_popString(int fd_write) |
{ |
{ |
char* s; |
char *s; |
cmo* m; |
cmo *err; |
|
cmo *m; |
|
|
#ifdef DEBUG |
#ifdef DEBUG |
fprintf(stderr, "ox_math:: opecode = SM_popString.\n"); |
fprintf(stderr, "ox_math:: opecode = SM_popString.\n"); |
#endif |
#endif |
|
|
if ((m = pop()) != NULL && (s = convert_cmo_to_string(m)) != NULL) { |
m = pop(); |
|
if (m->tag == CMO_STRING) { |
|
send_ox_cmo(fd_write, m); |
|
}else if ((s = new_string_set_cmo(m)) != NULL) { |
send_ox_cmo(fd_write, (cmo *)new_cmo_string(s)); |
send_ox_cmo(fd_write, (cmo *)new_cmo_string(s)); |
return 0; |
}else { |
|
err = make_error_object(SM_popString, m); |
|
send_ox_cmo(fd_write, err); |
} |
} |
return SM_popString; |
return 0; |
} |
} |
|
|
|
int local_execute(char *s) |
|
{ |
|
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; |
|
} |
|
|
/* この関数はサーバに依存する. */ |
/* この関数はサーバに依存する. */ |
int sm_executeStringByLocalParser(int fd_write) |
int sm_executeStringByLocalParser(int fd_write) |
{ |
{ |
cmo* m = NULL; |
symbol *symp; |
|
cmo* m = pop(); |
|
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 = pop()) != NULL && m->tag == CMO_STRING) { |
|
/* for mathematica */ |
if (m->tag == CMO_STRING |
/* mathematica に文字列を送って評価させる */ |
&& strlen(s = ((cmo_string *)m)->s) != 0) { |
MATH_evaluateStringByLocalParser(((cmo_string *)m)->s); |
if (s[0] == ':') { |
push(MATH_getObject2()); |
local_execute(++s); |
|
}else { |
|
/* for mathematica */ |
|
/* mathematica に文字列を送って評価させる */ |
|
ml_evaluateStringByLocalParser(s); |
|
ml_select(); |
|
push(receive_mlo()); |
|
} |
return 0; |
return 0; |
} |
} |
fprintf(stderr, "cannot execute: top of stack is not string!(%p, %d)\n", m, m->tag); |
#ifdef DEBUG |
|
symp = lookup_by_tag(m->tag); |
|
fprintf(stderr, "ox_math:: error. the top of stack is %s.\n", symp->key); |
|
#endif |
return SM_executeStringByLocalParser; |
return SM_executeStringByLocalParser; |
} |
} |
|
|
Line 362 int sm_executeFunction(int fd_write) |
|
Line 188 int sm_executeFunction(int fd_write) |
|
if ((m = pop()) == NULL || m->tag != CMO_INT32) { |
if ((m = pop()) == NULL || m->tag != CMO_INT32) { |
return SM_executeFunction; |
return SM_executeFunction; |
} |
} |
|
|
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++) { |
if ((argv[i] = pop()) == NULL) { |
argv[i] = pop(); |
return SM_executeFunction; |
|
} |
|
} |
} |
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日 */ |
#define VERSION 0x11121400 |
#define VERSION 0x11102700 |
#define ID_STRING "ox_math server 1999/12/14 15:25:00" |
#define ID_STRING "ox_math server 1999/10/28 17:29:25" |
|
|
|
int sm_mathcap(int fd_write) |
int sm_mathcap(int fd_write) |
{ |
{ |
cmo* c = make_mathcap_object(VERSION, ID_STRING); |
push(make_mathcap_object(VERSION, ID_STRING)); |
push(c); |
|
return 0; |
return 0; |
} |
} |
|
|
Line 393 int receive_sm_command(int fd_read) |
|
Line 217 int receive_sm_command(int fd_read) |
|
int execute_sm_command(int fd_write, int code) |
int execute_sm_command(int fd_write, int code) |
{ |
{ |
int err = 0; |
int err = 0; |
|
#ifdef DEBUG |
|
symbol *sp = lookup_by_tag(code); |
|
fprintf(stderr, "ox_math:: %s received.\n", sp->key); |
|
#endif |
|
|
switch(code) { |
switch(code) { |
case SM_popCMO: |
case SM_popCMO: |
Line 408 int execute_sm_command(int fd_write, int code) |
|
Line 236 int execute_sm_command(int fd_write, int code) |
|
err = sm_pops(fd_write); |
err = sm_pops(fd_write); |
break; |
break; |
case SM_executeStringByLocalParser: |
case SM_executeStringByLocalParser: |
|
case SM_executeStringByLocalParserInBatchMode: |
err = sm_executeStringByLocalParser(fd_write); |
err = sm_executeStringByLocalParser(fd_write); |
break; |
break; |
case SM_executeFunction: |
case SM_executeFunction: |
err = sm_executeFunction(fd_write); |
err = sm_executeFunction(fd_write); |
break; |
break; |
|
case SM_shutdown: |
|
shutdown(); |
|
break; |
case SM_setMathCap: |
case SM_setMathCap: |
pop(); /* 無視する */ |
pop(); /* 無視する */ |
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())); |
} |
} |
} |
} |