From rough989 at gmail.com  Fri Jan 26 11:30:36 2007
From: rough989 at gmail.com (Stefano Mazzuco)
Date: Fri, 26 Jan 2007 11:30:36 +0100
Subject: [Rd] double integral with C code
Message-ID: <45B9D84C.5080904@gmail.com>

Hi all,
This is more a C querie rather than a R one:
I'm writing a C code passing a function F to adapt fortran subroutine. I
need to integrate over two variables of F, call them x1 and x2. Then I
call the C code in R to optimize the integrated F function.
for example F could be defined as
---------------------------------------
static double marg_like(const double *param,
            double x1,
            double x2){...........}
---------------------------------------
Then I integrate over x1 and x2 with
-------------------------------------------------
F77_CALL(adapt)(2,(-5,-5),(5,5),100,1700,F,0.01,10000)
-------------------------------------------------

So here my question: I should I define x1 and x2? For the time being I
defined them as static variables, i.e.
static double x1;
static double x2;

but I'm pretty sure this is wrong

Any hint?
thanks in advance
Stefano

P.S. Sorry for having posted this in the R-help section, was an unintended mistake.



[Rd] double integral with C code
Stefano Mazzuco