bugGNU Octave - Bugs: bug #49427, structure returned by ode solvers...


bug #49427: structure returned by ode solvers invoked with one output is incomplete

Submitted by:  Carlo de Falco <cdf>
Submitted on:  Mon 24 Oct 2016 10:41:06 AM UTC  
Category: Octave FunctionSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Matlab Compatibility
Status: NoneAssigned to: None
Originator Name: Open/Closed: Open
Release: devOperating System: Any

Add a New Comment(Rich Markup)

You are not logged in

Please log in, so followups can be emailed to you.


Thu 13 Jul 2017 06:47:21 AM UTC, comment #3:

I'm interested in giving it a try. I've looked a little at it, and have almost made the jump to making my own version of it. Any help on where to start would be great.

Marshall <marsian>
Wed 12 Jul 2017 02:23:29 PM UTC, comment #2:

sure, it shouldn't take too much effort to implement deval at the moment, are you willing to help with this?

If so, I can send you precise pointers to the current code
where to start from.

Carlo de Falco <cdf>
Project Member
Wed 12 Jul 2017 01:51:40 PM UTC, comment #1:

I'm interested in having a deval function. Since odeXX is native to Octave now and does have the output form of sol = odeXX(...), it seems like what's there should be enough to extract deval from the existing functions: It already interpolates internally.

Marshall <marsian>
Mon 24 Oct 2016 10:41:06 AM UTC, original submission:

The ode solvers in Matlab can be called with two main different
forms for the output


the second form does not merely aggregate outputs in one
structure but triggers a completely different behaviour of the

The idea is that with the second form the solver advances until
the end of the timespan (or passed it) and does not care of
evaluation of the state at intermediate timesteps.

The output structure "out" contains all data that can be fed
to "deval" to evaluate the solution at any time point between
tspan(1) and tspan(end) to the maximum possible order of
accuracy consistent with the solution method used.

We don't have deval in Octave at the moment and I cannot find
documentation of the actual fields of the output structure in
matlab, so I think we don't need to copy those fields exactly.

On the other hand here is the list of data we would need to do
interpolation as a post-processing step (at least for explicit runge kutta):

  • times selected by the adaptive selection algorithm and state evaluated there
  • order of the method (to select interpolation algorithm)
  • name of the method (to select location of intermediate "stages")
  • derivatives at intermadiate "stages"
Carlo de Falco <cdf>
Project Member


(Note: upload size limit is set to 16384 kB, after insertion of the required escape characters.)

Attach File(s):

No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -unavailable- added by marsian (Posted a comment)
  • -unavailable- added by cdf (Submitted the item)

    Do you think this task is very important?
    If so, you can click here to add your encouragement to it.
    This task has 0 encouragements so far.

    Only project members can vote.


    Please enter the title of George Orwell's famous dystopian book (it's a date):



    No Changes Have Been Made to This Item

    Back to the top

    Powered by Savane 3.1-cleanup1