version 1.1, 2018/03/29 11:52:18 |
version 1.2, 2018/03/30 04:43:16 |
|
|
|
/* $OpenXM$ |
|
*/ |
//#include <gsl/gsl_types.h> |
//#include <gsl/gsl_types.h> |
//#include <gsl/gsl_sys.h> |
//#include <gsl/gsl_sys.h> |
|
#include <unistd.h> |
#include <gsl/gsl_sf_result.h> |
#include <gsl/gsl_sf_result.h> |
|
#include <gsl/gsl_errno.h> |
#include <gsl/gsl_sf_gamma.h> |
#include <gsl/gsl_sf_gamma.h> |
#include "ox_gsl.h" |
#include "ox_gsl.h" |
extern int Debug; |
extern int Debug; |
|
// local prototype declarations |
|
|
void call_gsl_sf_lngamma_complex_e() { |
void call_gsl_sf_lngamma_complex_e() { |
cmo *c; |
cmo *c; |
double zr; |
double zr; |
Line 13 void call_gsl_sf_lngamma_complex_e() { |
|
Line 19 void call_gsl_sf_lngamma_complex_e() { |
|
int status; |
int status; |
cmo *r[3]; |
cmo *r[3]; |
cmo *ans; |
cmo *ans; |
gsl_set_error_handler_off(); |
// gsl_set_error_handler_off(); |
// Todo, gsl_set_error_handler(my_handler); |
gsl_set_error_handler((gsl_error_handler_t *)myhandler); |
c = pop(); // number of args |
c = pop(); // number of args |
zr=get_double(); |
zr=get_double(); |
zi=get_double(); |
zi=get_double(); |
Line 23 void call_gsl_sf_lngamma_complex_e() { |
|
Line 29 void call_gsl_sf_lngamma_complex_e() { |
|
r[0] = (cmo *)new_cmo_double(lnr.val); |
r[0] = (cmo *)new_cmo_double(lnr.val); |
r[1] = (cmo *)new_cmo_double(arg.val); |
r[1] = (cmo *)new_cmo_double(arg.val); |
r[2] = (cmo *)new_cmo_int32(status); |
r[2] = (cmo *)new_cmo_int32(status); |
ans = (cmo *)new_cmo_list_array((void **)r,3); |
ans = (cmo *)new_cmo_list_array((void *)r,3); |
push(ans); |
push(ans); |
} |
} |