[Octave-bug-tracker] [bug #47756] [ode suit] integrate_adaptive fails with vertical tspan

[Juan Pablo Carbajal]

URL:
  <http://savannah.gnu.org/bugs/?47756>

                 Summary: [ode suit] integrate_adaptive fails with vertical
tspan
                 Project: GNU Octave
            Submitted by: juanpi
            Submitted on: Thu 21 Apr 2016 07:12:39 PM GMT
                Category: Octave Function
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Regression
                  Status: None
             Assigned to: None
         Originator Name: juanpi
        Originator Email: 
             Open/Closed: Open
         Discussion Lock: Any
                 Release: dev
        Operating System: Any

    _______________________________________________________

Details:

When tspan is an array of vertical timestamps integrate_adaptive.m fails when
attempting hermite interpolation, due to the size of the argument
tspan(t_caught).

Example

pf   = [-0, -50, 0.25, 0.5, 0.65];
func =@(t,x) pf(1) * x(:) + pf(2) * (x(:) - pf(3)) .* (x(:) - pf(4)) .* (x(:)
- pf(5));
opt  = odeset ('RelTol', 1e-6, ...
                'AbsTol', 1e-6, ...
                'InitialStep', 1e-3, ...
                'MaxStep', t(end) / 10);

t      = linspace (0,1,100).';
[t_ y] = ode45 (func, t, (-1:0.5:1).', opt);


The patch attached should fix the problem





    _______________________________________________________

Reply to this item at:

  <http://savannah.gnu.org/bugs/?47756>

_______________________________________________
  Message sent via/by Savannah
  http://savannah.gnu.org/