Re: GSoC ode15{i,s} - type conversion and approximation problems

[Carlo De Falco]

On 16 Jun 2016, at 17:25, Francesco Faccio <address@hidden> wrote:

> Using sparse methods, in order to use Sundials' sparse solvers, I must 
> provide both DF/DY and DF/DYP in CSC format (it cannot be done automatically 
> with different quotient by Sundials), so I need to approximate one of them if 
> the user passes only one matrix. Again, should I write my own approximation 
> or use an external library (and which could be better to use for this 
> problem)?

I think it is OK for now to use a sparse method only if the user provides a 
sparse jacobian.
matlab also has the option to specify the pattern of a sparse jacobian and to 
approximate its entries
by finite differences but we can implement this feature at a later time.

c.