bugGNU Octave - Bugs: bug #49890, error if event function stops ode45

 
 

bug #49890: error if event function stops ode45

Submitted by:  None
Submitted on:  Tue 20 Dec 2016 04:50:30 AM UTC  
 
Category: NoneSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Regression
Status: FixedAssigned to: Carlo de Falco <cdf>
Originator Name: BP MinakerOriginator Email: -unavailable-
Open/Closed: ClosedRelease: dev
Operating System: GNU/Linux

Add a New Comment(Rich Markup)
   

You are not logged in

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

 

Mon 26 Dec 2016 04:33:12 PM UTC, comment #5:

I fixed the issued with ode23, ode45 on the stable branch here (http://hg.savannah.gnu.org/hgweb/octave/rev/2b2ac017c901).

Rik <rik5>
Project Administrator
Fri 23 Dec 2016 01:23:21 PM UTC, comment #4:

I was hoping the submitter would provide an example so hadn't commented on this bug earlier. The error message reported by the user can be reproduced by changing the ode45 call in the sample below to

[t,y,te,ye] = ode45 (@(t,y) 2, [0, 2], 1, opts)

It is valid in matlab to have less than 5 return arguments when events are included. A fix for this and several other odeIJ bugs
is in my repository here

https://bitbucket.org/wgreene71/octave

Bill Greene <billgreene>
Fri 23 Dec 2016 06:03:25 AM UTC, comment #3:

Without an example from the original reporter, we can close this bug. It works for me with a trivial event function

and the call

Rik <rik5>
Project Administrator
Tue 20 Dec 2016 10:23:26 AM UTC, comment #2:

... or can you send a minimal example? That would help us debugging it and finally to add a test.

Seb.

Sebastian <sschoeps>
Tue 20 Dec 2016 09:13:28 AM UTC, comment #1:

Thanks for reporting this.

Do you happen to also know a solution?
I'm cc-ing the odepkg package maintainers as this may also affect solvers in that package.

c.

Carlo de Falco <cdf>
Project MemberIn charge of this item.
Tue 20 Dec 2016 04:50:30 AM UTC, original submission:

If an event function triggers a stop in ode45, the result is an error: element number 1 undefined in return list
It appears that ode45 fails to return any solution arguments if it is stopped before the end of the solution interval.

Anonymous

 

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

Attach File(s):
   
   
Comment:
   

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 billgreene (Posted a comment)
  • -unavailable- added by rik5 (Posted a comment)
  • -unavailable- added by sschoeps (Posted a comment)
  • -unavailable- added by cdf (Posted a comment)
  • -unavailable- added by cdf
  • -unavailable- added by cdf
  • -unavailable- added by cdf
  • -unavailable- added by cdf
  • -unavailable- added by None (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):

     

     

    Follow 8 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Mon 26 Dec 2016 04:33:12 PM UTCrik5StatusWorks For Me=>Fixed
      Open/ClosedOpen=>Closed
    Fri 23 Dec 2016 06:03:25 AM UTCrik5StatusNone=>Works For Me
    Tue 20 Dec 2016 09:13:28 AM UTCcdfAssigned toNone=>cdf
      Carbon-Copy-=>Added -unavailable-
      Carbon-Copy-=>Added -unavailable-
      Carbon-Copy-=>Added -unavailable-
      Carbon-Copy-=>Added -unavailable-

    Back to the top


    Powered by Savane 3.1-cleanup1