Return to call_gsl.c CVS log

Up to [local] / OpenXM / src / ox_gsl

Annotation of OpenXM/src/ox_gsl/call_gsl.c, Revision 1.1

1.1     ! takayama    1: //#include <gsl/gsl_types.h>
        !             2: //#include <gsl/gsl_sys.h>
        !             3: #include <gsl/gsl_sf_result.h>
        !             4: #include <gsl/gsl_sf_gamma.h>
        !             5: #include "ox_gsl.h"
        !             6: extern int Debug;
        !             7: void  call_gsl_sf_lngamma_complex_e() {
        !             8:   cmo *c;
        !             9:   double zr;
        !            10:   double zi;
        !            11:   gsl_sf_result lnr;
        !            12:   gsl_sf_result arg;
        !            13:   int status;
        !            14:   cmo *r[3];
        !            15:   cmo *ans;
        !            16:   gsl_set_error_handler_off();
        !            17:   // Todo, gsl_set_error_handler(my_handler);
        !            18:   c = pop(); // number of args
        !            19:   zr=get_double();
        !            20:   zi=get_double();
        !            21:   if (Debug) printf("gsl_sf_lngamma_complex_e(zr=%lg,zi=%lg)\n",zr,zi);
        !            22:   status = gsl_sf_lngamma_complex_e(zr,zi,&lnr,&arg);
        !            23:   r[0] = (cmo *)new_cmo_double(lnr.val);
        !            24:   r[1] = (cmo *)new_cmo_double(arg.val);
        !            25:   r[2] = (cmo *)new_cmo_int32(status);
        !            26:   ans = (cmo *)new_cmo_list_array((void **)r,3);
        !            27:   push(ans);
        !            28: }