Hi,
I am trying to solve a model that consists of rather stiff ODEs in R.
I use the package ODEsolve (lsoda) to solve these ODEs.
To speed up the integration, the jacobian is also specified.
Basically, the model is a one-dimensional advection-diffusion problem,
and thus the jacobian is a tridiagonal matrix.
The size of this jacobian is 100*100.
In the original package LSODA it is possible to specify that the
jacobian is banded, which makes its inversion very efficient.
However, this feature seems to have been removed in the R version.
Is there a way to overcome this limitation?
Thanks
dr. Karline Soetaert
NIOO - CEME
PO box 140
4400 AC Yerseke
the Netherlands
Phone: ++ 31 113 577487
fax: ++ 31 113 573616
e-mail: [EMAIL PROTECTED]
[[alternative HTML version deleted]]
______________________________________________
R-help@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
