[Octave-bug-tracker] [bug #47350] installed package odepkg, fails to use option "JPattern" correctly

[anonymous]

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

                 Summary: installed package odepkg, fails to use option
"JPattern" correctly
                 Project: GNU Octave
            Submitted by: None
            Submitted on: Sun 06 Mar 2016 02:02:00 AM UTC
                Category: Octave Forge Package
                Severity: 3 - Normal
                Priority: 5 - Normal
              Item Group: Crash
                  Status: None
             Assigned to: None
         Originator Name: Rob Carter
        Originator Email: address@hidden
             Open/Closed: Open
         Discussion Lock: Any
                 Release: 4.0.0
        Operating System: Microsoft Windows

    _______________________________________________________

Details:

I installed the latest version of odepkg from Ocatve forge.  I am trying to
use ode23s, but first I use odeset to set up options for ode23s.  My problem
is that my program fails due to the way ode23a is using the option 'JPattern'.
 See my code snippet below...


M = spalloc(nodes, nodes, 2*nodes);

%Populate Jacobian pattern to ease ode23s' job a bit
jac_patt = M;
for i = 3:nodes-1
  jac_patt(i, i+1) = 1;
  jac_patt(i, i) = 1;
  jac_patt(i, i-1) = 1;
  jac_patt(i, i-2) = 1;
end
jac_patt(1, 1) = 1;
jac_patt(1, 2) = 1;
jac_patt(2, 1) = 1;
jac_patt(2, 2) = 1;
jac_patt(2, 3) = 1;
jac_patt(nodes, nodes) = 1;
jac_patt(nodes, nodes-1) = 1;
jac_patt(nodes, nodes-2) = 1;

options = odeset('OutputFcn',@query_solver,
                 'JPattern',jac_patt,
                 'MaxStep',max_step_time);

%Call the ODE solver
[T, C] = ode23s(@ode_statements,time_span,initial_conds,options);



I get the following error message when trying to run...

error: element number 2 undefined in return list
error: called from
    ode23s at line 89 column 12
    pde_solution at line 68 column 6
    main_program at line 15 column 14

Please help me understand how to use this option

Thanks,
Rob Carter




    _______________________________________________________

Reply to this item at:

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

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