[R] Problem defining a system of odes as a C library with lsoda

[Woodrow Setzer]

I think the problem is in odesolve (something I thought I'd already fixed).  I don't have a 64-bit system to test on, so could you try these changes and let me know (offline) if they fix the problem?  There is a type mismatch in odesolve; I want to know if that is the cause of your problem.  In any case, I'll get an updated version of odesolve on CRAN ASAP.

in mymod.c:
change the line 
mymod(void(* odeparms)(int *, double *))
to
mymod(void(* odeparms)(long int *, double *))

and the line
int N=3;
to 
long int N=3;

Woody

-----Original Message-----
From: Dylan Childs <D.Childs at sheffield.ac.uk>
Sent: Nov 6, 2005 2:39 PM
To: r-help at stat.math.ethz.ch
Subject: [R] Problem defining a system of odes as a C library with lsoda

I have been trying to make use of the odesolve library on my 
university's Linux grid - currently R version 2.0.1 is installed and 
the system runs 64-bit Scientific Linux based on Redhat.  

... [deleted]

However, the call to 
lsoda fails with the following error:

Error in lsoda(c(1, 0, 0), times, "myderivs", parms, rtol = 1e-04, atol 
= my.atol,  :
	Confusion over the length of parms

... [deleted]

Dr. Dylan Z. Childs
Department of Animal and Plant Sciences,
University of Sheffield,
Sheffield,
S10 2TN, UK.

______________________________________________
R-help at stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html


Woodrow Setzer
National Center for Computational Toxicology
US Environmental Protection Agency
Research Triangle Park, NC 27711