[BACK]Return to call_gsl_sf.c CVS log [TXT][DIR] Up to [local] / OpenXM / src / ox_gsl

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

1.1     ! takayama    1: /* $OpenXM$
        !             2: */
        !             3: //#include <gsl/gsl_types.h>
        !             4: //#include <gsl/gsl_sys.h>
        !             5: #include <unistd.h>
        !             6: #include <gsl/gsl_sf_result.h>
        !             7: #include <gsl/gsl_errno.h>
        !             8: #include <gsl/gsl_sf_gamma.h>
        !             9: #include "ox_gsl.h"
        !            10: extern int Debug;
        !            11: // local prototype declarations
        !            12: #include "call_gsl_sf.h"
        !            13:
        !            14: void  call_gsl_sf_gamma_inc() {
        !            15:   int argc;
        !            16:   double a;
        !            17:   double x;
        !            18:   double val;
        !            19:   cmo *ans;
        !            20:   //  gsl_set_error_handler_off();
        !            21:   gsl_set_error_handler((gsl_error_handler_t *)myhandler);
        !            22:   argc = get_i(); // number of args
        !            23:   if (argc != 2) {
        !            24:     pops(argc);
        !            25:     push(make_error2("The argc must be 2 for gsl_sf_gamma_inc.",NULL,0,-1));
        !            26:     return;
        !            27:   }
        !            28:   a=get_double();
        !            29:   x=get_double();
        !            30:   if (Debug) printf("gsl_sf_gamma_inc(a=%lg,x=%lg)\n",a,x);
        !            31:   val = gsl_sf_gamma_inc(a,x);
        !            32:   ans = (cmo *)new_cmo_double(val);
        !            33:   push(ans);
        !            34: }

FreeBSD-CVSweb <freebsd-cvsweb@FreeBSD.org>