bugGNU Octave - Bugs: bug #49871, ode15s, MStateDependence option...

 
 

bug #49871: ode15s, MStateDependence option incorrect

Submitted by:  Bill Greene <billgreene>
Submitted on:  Sat 17 Dec 2016 06:58:17 PM UTC  
 
Category: Octave FunctionSeverity: 3 - Normal
Priority: 5 - NormalItem Group: Incorrect Result
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.

 

Sat 17 Dec 2016 07:52:16 PM UTC, comment #2:

I don't know of any way that Sundials/IDA could use the
information from the MStateDependence option. It definitely
shouldn't affect the Jacobian (dF/dy) calculation. My inclination
would be to simply ignore that option.

Bill Greene <billgreene>
Sat 17 Dec 2016 07:15:03 PM UTC, comment #1:

Bill,

You seem to be experienced with sundials,
do you know the ida option to which MStateDependence
should be mapped?

Carlo de Falco <cdf>
Project Member
Sat 17 Dec 2016 06:58:17 PM UTC, original submission:

I have just started doing some testing of the new ode15s function.

I have discovered incorrect behavior when the Jacobian option is
specified. I get the following warning message:

warning: with MStateDependence != 'none' an internal approximation of Jacobian Matrix will be used. Set MStateDependence equal to 'none' if you want to provide a constant or time-dependent Jacobian
warning: called from
ode15s at line 258 column 7
test_ode15s_jac at line 23 column 13

I believe that ode15s.m is, for some reason, using the setting of
MStateDependence to ignore my Jacobian function. This is the
wrong interpretation of MStateDependence. My understanding is
that MStateDependence in MATLAB just serves as a "hint" to the code that it may be able to adjust the algorithm slightly to improve performance.

I am attaching a small script that reproduces this behavior.

Bill Greene <billgreene>

 

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

Attach File(s):
   
   
Comment:
   

Attached Files
file #39260:  test_ode15s_jac.m added by billgreene (1KiB - text/plain)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -unavailable- added by cdf (Posted a comment)
  • -unavailable- added by cdf
  • -unavailable- added by cdf
  • -unavailable- added by billgreene (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 3 latest changes.

    Date Changed By Updated Field Previous Value => Replaced By
    Sat 17 Dec 2016 07:13:08 PM UTCcdfCarbon-Copy-=>Added -unavailable-
      Carbon-Copy-=>Added -unavailable-
    Sat 17 Dec 2016 06:58:17 PM UTCbillgreeneAttached File-=>Added test_ode15s_jac.m, #39260

    Back to the top


    Powered by Savane 3.1-cleanup1