bugGNU Octave - Bugs: bug #47350, odepkg: ode23s fails to use option...

 
 

bug #47350: odepkg: ode23s fails to use option "JPattern" correctly

Submitter:  None
Submitted:  Sun 06 Mar 2016 02:02:00 AM UTC
   
 
Category:  Octave Package Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Incorrect Result
Status:  Duplicate Assigned to:  None
Originator Name:  Rob Carter Originator Email:  -email is unavailable-
Open/Closed:  * Closed Release:  * 4.0.0
Operating System:  * Microsoft Windows Fixed Release:  None
Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Wed 22 Jun 2016 05:37:24 AM UTC, comment #1: 

Thanks for your bug report, and sorry that no one has responded to it yet. A developer has noticed this missing feature independently and submitted a patch on another bug report, so I'm marking this as a duplicate and closing. Please test the patch on bug #48289 if you can to see if it resolves your problem, and comment there if necessary.

Mike Miller <mtmiller>
Group Member
Sun 06 Mar 2016 02:02:00 AM UTC, original submission:  

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

Anonymous

 

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

Attach Files:
   
   
Comment:
   

No files currently attached

 

Digest:
   bug dependencies.

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by mtmiller (Updated the item)
  • -email is unavailable- added by None (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.

     

    Follow 5 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2016-06-22 mtmiller StatusNone Duplicate
        Open/ClosedOpen Closed
        Dependencies- Depends on bugs #48289
    2016-03-22 mtmiller Item GroupSegfault, Bus Error, etc. Incorrect Result
        Summaryinstalled package odepkg, fails to use option &quot;JPattern&quot; correctly odepkg: ode23s fails to use option "JPattern" correctly

    Back to the top

    Powered by Savane 3.13-3230.
    Corresponding source code