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: k.Soetaert@nioo.knaw.nl

 

 


	[[alternative HTML version deleted]]


