octave-maintainers
[Top][All Lists]
Advanced

[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: Review ode csets


From: c.
Subject: Re: Review ode csets
Date: Fri, 21 Oct 2016 19:02:06 +0200

On 21 Oct 2016, at 12:39, Carlo de Falco <address@hidden> wrote:

> Rik,
> This issue does not look like it is bad enough to block the release.

Actually, this probably is worse than I thought and deserves a fix.
It seems refinement steps are not returned whether integration is stopped or 
not:

Matlab 2016a:
------------------------------------------------------------------------------------------------------
>> [t, y] = ode45 (myfun, [0 1], 1, odeset ('maxstep', 1, 'initialstep', .1, 
>> 'refine', 3)); numel (t)

ans =

    10

>> [t, y] = ode45 (myfun, [0 1], 1, odeset ('maxstep', 1, 'initialstep', .1, 
>> 'refine', 1)); numel (t)

ans =

     4
------------------------------------------------------------------------------------------------------

Octave:
------------------------------------------------------------------------------------------------------
>> [t, y] = ode45 (myfun, [0 1], 1, odeset ('maxstep', 1, 'initialstep', .1, 
>> 'refine', 3)); numel (t)
ans =  6
>> [t, y] = ode45 (myfun, [0 1], 1, odeset ('maxstep', 1, 'initialstep', .1, 
>> 'refine', 1)); numel (t)
ans =  6
>> 
------------------------------------------------------------------------------------------------------

what do you think?
I'll post a bug report about this later.
c.






reply via email to

[Prev in Thread] Current Thread [Next in Thread]