version 1.4, 2018/03/30 08:48:23 |
version 1.8, 2018/04/17 00:56:38 |
|
|
/* $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.7 2018/04/06 01:56:49 takayama Exp $ |
*/ |
*/ |
|
|
#include <stdio.h> |
#include <stdio.h> |
Line 40 void gc_free(void *p,size_t size) |
|
Line 40 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 111 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 293 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 307 int sm_executeFunction() |
|
Line 323 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 { |
}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; |