Re: 1D PDE solver for Octave

[Carlo De Falco]

On 6 Jul 2016, at 13:42, Bill Greene <address@hidden> wrote:

> lsode wouldn't handle the DAE that results when some of the PDE are elliptic.
> 
> I believe that a daspk implementation would be much slower than the current 
> one
> because it would use a full-matrix jacobian whereas Sundials ida is exploiting
> the narrow bandwidth of the matrix.

If you had a complete implementation based on daspk done, switching to ode15s 
at a later time would be trivial.

If you cannot wait for ode15s to be finished, you could help with its 
development.

c.