|
From: | Bill Greene |
Subject: | Re: 1D PDE solver for Octave |
Date: | Wed, 6 Jul 2016 12:30:05 -0400 |
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.
[Prev in Thread] | Current Thread | [Next in Thread] |