Displaying 2 results from an estimated 2 matches for "jacfunc".
Did you mean:
cfunc
2004 Oct 06
1
Foreign code problem
Hello,
I wanted to test the odesolve package and tried to use compiled C-code.
But when I do:
erg <- lsoda(y, times, "mond", parms, rtol, atol, tcrit=NULL, jacfunc=NULL,
verbose=FALSE, dllname="mond", hmin=0, hmax=Inf)
I get the error message:
Error in lsoda(y, times, "mond", parms, rtol, atol, tcrit = NULL,
jacfunc = NULL, :
Unable to find mond in mond
The C code is:
#include <math.h...
2005 Nov 06
1
Problem defining a system of odes as a C library with lsoda
...ms <- c(k1 = 0.04, k2 = 1e4, k3=3e7)
my.atol <- c(1e-6, 1e-10, 1e-6)
times <- c(0, 4*10^(-1:10))
print(system.time(
out1 <- lsoda(c(1.0,0.0,0.0),times,"myderivs",
parms,
rtol=1e-4,atol=my.atol,jacfunc="myjac",
dllname="mymod")
)
)
### END
When I run this mymod.c compiles without any errors or warnings and the
mymod.so library seems to load fine via dyn.load (the call
is.loaded("mymod") returns TRUE for ex...