[Octave-patch-tracker] [patch #8102] wrappers for stiff ode solvers

[Pooja Rao]

Follow-up Comment #1, patch #8102 (project octave):

This is an initial changeset for function ode15s (from jwe's implementation
patch #8102) to the ode-pkg. I have added two options ('InitialStep' and
'MaxStep') to this function along with changing the way it outputs its
solution. Earlier, it would take [tstart tfinal] as the time input vector and
produce output only at tstart and tfinal (because it uses LSODE /DASRT). To
make it produce solution values at points between tstart and tfinal, I used
linspace in the ode15s file. This is a temporary workaround.  I also added a
test and a couple of demos at the end.

The following is a brief list of TO-DOs for ode15s:

1) Add options to make it more Matlab compatible
2) Replace linspace with a more adaptive output
3) Replace LSODE with DASPK/DASRT  so it solves DAEs.
In its current form it uses DASRT but only when the 'Events'  option is
non-empty.

Thanks!

(file #31026)
    _______________________________________________________

Additional Item Attachment:

File name: added_ode15s.diff              Size:6 KB


    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/patch/?8102>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/