version 1.4, 2018/03/30 08:48:23 |
version 1.12, 2018/06/07 01:53:33 |
|
|
/* $OpenXM: OpenXM/src/ox_gsl/ox_gsl.c,v 1.3 2018/03/30 04:43:16 takayama Exp $ |
/* $OpenXM: OpenXM/src/ox_gsl/ox_gsl.c,v 1.11 2018/06/06 07:40:32 takayama Exp $ |
*/ |
*/ |
|
|
#include <stdio.h> |
#include <stdio.h> |
|
|
#include <setjmp.h> |
#include <setjmp.h> |
#include <string.h> |
#include <string.h> |
#include <unistd.h> |
#include <unistd.h> |
|
#include <signal.h> |
#include <math.h> |
#include <math.h> |
#include "ox_gsl.h" |
#include "ox_gsl.h" |
#include "call_gsl.h" // need only when you bind call_gsl functions. |
#include "call_gsl.h" // need only when you bind call_gsl functions. |
Line 40 void gc_free(void *p,size_t size) |
|
Line 41 void gc_free(void *p,size_t size) |
|
void init_gc() |
void init_gc() |
{ GC_INIT(); |
{ GC_INIT(); |
mp_set_memory_functions(GC_malloc,gc_realloc,gc_free); |
mp_set_memory_functions(GC_malloc,gc_realloc,gc_free); |
|
init_dic(); // initialize ox_eval.c |
} |
} |
|
|
void initialize_stack() |
void initialize_stack() |
Line 110 int sm_mathcap() |
|
Line 112 int sm_mathcap() |
|
// CMO_DISTRIBUTED_POLYNOMIAL, |
// CMO_DISTRIBUTED_POLYNOMIAL, |
// CMO_RECURSIVE_POLYNOMIAL, |
// CMO_RECURSIVE_POLYNOMIAL, |
// CMO_POLYNOMIAL_IN_ONE_VARIABLE, |
// CMO_POLYNOMIAL_IN_ONE_VARIABLE, |
|
CMO_TREE, |
CMO_ERROR2, |
CMO_ERROR2, |
0}; |
0}; |
int available_sm_command[]={ |
int available_sm_command[]={ |
Line 291 char *get_string() { |
|
Line 294 char *get_string() { |
|
return(NULL); |
return(NULL); |
} |
} |
|
|
|
void test_ox_eval() { |
|
cmo *c; |
|
double d=0; |
|
pop(); |
|
c=pop(); |
|
if (Debug) { |
|
ox_printf("cmo *c="); print_cmo(c); ox_printf("\n"); |
|
} |
|
init_dic(); |
|
register_entry("x",1.25); |
|
if (eval_cmo(c,&d) == 0) myhandler("eval_cmo failed",NULL,0,-1); |
|
push((cmo *)new_cmo_double(d)); |
|
} |
|
|
int sm_executeFunction() |
int sm_executeFunction() |
{ |
{ |
cmo_string *func = (cmo_string *)pop(); |
cmo_string *func = (cmo_string *)pop(); |
Line 298 int sm_executeFunction() |
|
Line 315 int sm_executeFunction() |
|
push(make_error2("sm_executeFunction, not CMO_STRING",NULL,0,-1)); |
push(make_error2("sm_executeFunction, not CMO_STRING",NULL,0,-1)); |
return -1; |
return -1; |
} |
} |
|
init_dic(); |
// Test functions |
// Test functions |
if (strcmp(func->s, "add_int32") == 0) { |
if (strcmp(func->s, "add_int32") == 0) { |
my_add_int32(); |
my_add_int32(); |
Line 307 int sm_executeFunction() |
|
Line 325 int sm_executeFunction() |
|
show_double_list(); |
show_double_list(); |
}else if (strcmp(func->s,"restart")==0) { |
}else if (strcmp(func->s,"restart")==0) { |
pop(); restart(); |
pop(); restart(); |
|
}else if (strcmp(func->s,"test_ox_eval")==0) { |
|
test_ox_eval(); |
// The following functions are defined in call_gsl.c |
// The following functions are defined in call_gsl.c |
}else if (strcmp(func->s,"gsl_sf_lngamma_complex_e")==0) { |
}else if (strcmp(func->s,"gsl_sf_lngamma_complex_e")==0) { |
call_gsl_sf_lngamma_complex_e(); |
call_gsl_sf_lngamma_complex_e(); |
|
}else if (strcmp(func->s,"gsl_integration_qags")==0) { |
|
call_gsl_integration_qags(); |
|
}else if (strcmp(func->s,"gsl_monte_plain_integrate")==0) { |
|
call_gsl_monte_plain_miser_vegas_integrate(0); |
|
}else if (strcmp(func->s,"gsl_monte_miser_integrate")==0) { |
|
call_gsl_monte_plain_miser_vegas_integrate(1); |
|
}else if (strcmp(func->s,"gsl_monte_vegas_integrate")==0) { |
|
call_gsl_monte_plain_miser_vegas_integrate(2); |
}else { |
}else { |
push(make_error2("sm_executeFunction, unknown function",NULL,0,-1)); |
push(make_error2("sm_executeFunction, unknown function",NULL,0,-1)); |
return -1; |
return -1; |
|
|
fd_rw = oxf_open(3); |
fd_rw = oxf_open(3); |
oxf_determine_byteorder_server(fd_rw); |
oxf_determine_byteorder_server(fd_rw); |
} |
} |
|
#if defined(__CYGWIN__) |
|
void *mysignal(int sig,void (*handler)(int m)); |
|
mysignal(SIGUSR1,usr1_handler); |
|
#else |
signal(SIGUSR1,usr1_handler); |
signal(SIGUSR1,usr1_handler); |
|
#endif |
|
|
while(1) { |
while(1) { |
receive(); |
receive(); |