diff scripts/ode/ode15i.m @ 22902:284bbb0328f2

Update odei5{i,s} code for Octave 4.3.0+ and Sundials 2.7.0. * libinterp/dldfunctions/__ode15__.cc : Specify CSC format for sparse matrices, use new API function name for clearing sparse matrix. * scripts/ode/ode15{i,s}.m : use correct calling sequence for odemergeoptions.
author Carlo de Falco <carlo.defalco@polimi.it>
date Tue, 01 Nov 2016 06:41:29 +0100
parents 4c56f3ffec04
children 23847979b91e
line wrap: on
line diff
--- a/scripts/ode/ode15i.m	Tue Aug 23 03:19:11 2016 +0200
+++ b/scripts/ode/ode15i.m	Tue Nov 01 06:41:29 2016 +0100
@@ -151,27 +151,24 @@
     endif
   endif
 
-  persistent defaults   = [];
-  persistent classes    = [];
-  persistent attributes = [];
-
-  [defaults, classes, attributes] = odedefaults (n, trange(1),
-                                                 trange(end));
+  [defaults, classes, attributes] = ...
+  odedefaults (n, trange(1), trange(end));
 
-  defaults   = rmfield (defaults,   {"NonNegative", "Mass", ...
-                                     "MStateDependence", "MvPattern", ...
-                                     "MassSingular", "InitialSlope", "BDF"});
-  classes    = rmfield (classes,    {"NonNegative", "Mass", ...
-                                     "MStateDependence", "MvPattern", ...
-                                     "MassSingular", "InitialSlope", "BDF"});
-  attributes = rmfield (attributes, {"NonNegative", "Mass", ...
-                                     "MStateDependence", "MvPattern", ...
-                                     "MassSingular", "InitialSlope", "BDF"});
+  persistent ignorefields = {"NonNegative", "Mass", ...
+                             "MStateDependence", "MvPattern", ...
+                             "MassSingular", "InitialSlope", "BDF"};
+  
+  defaults   = rmfield (defaults, ignorefields);
+  classes    = rmfield (classes, ignorefields);
+  attributes = rmfield (attributes, ignorefields);
 
   classes        = odeset (classes, 'Vectorized', {});
-  attributes     = odeset (attributes, 'Jacobian', {}, 'Vectorized', {});
+  attributes     = ...
+  odeset (attributes, 'Jacobian', {}, 'Vectorized', {});
 
-  options = odemergeopts (options, defaults, classes, attributes, solver);
+  options = ...
+  odemergeopts ("ode15i", options, defaults,
+                classes, attributes, solver);
 
   ## Jacobian
   options.havejac       = false;