bugGNU Octave - Bugs: bug #50193, ode15s gives wrong result for...

 
 

bug #50193: ode15s gives wrong result for complex input

Submitter:  Marco Caliari <caliari>
Submitted:  Tue 31 Jan 2017 10:18:37 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  None Assigned to:  None
Originator Name:  Open/Closed:  * Open
Release:  * dev Operating System:  * GNU/Linux
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Thu 16 Mar 2017 04:06:02 PM UTC, comment #2: 

The computation should be interrupted whenever fun is discovered to be complex, for instance in check_default_input.

Marco Caliari <caliari>
Group Member
Thu 09 Mar 2017 05:02:37 PM UTC, comment #1: 

This is a known problem:

sundials/IDA can only deal with float or double
and either of the two must be chosen when compiling
the library.


What to do for complex input is unclear to me at the moment.

Carlo de Falco <cdf>
Group Member
Tue 31 Jan 2017 10:18:37 AM UTC, original submission:  

The following code


f=@(t,y) 1i*y;sol=ode15s(f,[0,1],1)


gives a completely wrong result (which should be exp(1i*sol.x).

Marco Caliari <caliari>
Group Member

 

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

Attach Files:
   
   
Comment:
   

No files currently attached

 

Depends on the following items: None found

Items that depend on this one: None found

 

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

    There are 0 votes so far. Votes easily highlight which items people would like to see resolved in priority, independently of the priority of the item set by tracker managers.

    Only group members can vote.

     

    No changes have been made to this item

    Back to the top

    Powered by Savane 3.13-3230.
    Corresponding source code