bugGNU Octave - Bugs: bug #46074, missing function ode23tb

 
 

bug #46074: missing function ode23tb

Submitter:  Axel <axelt83>
Submitted:  Tue 29 Sep 2015 08:23:25 AM UTC
   
 
Category:  Octave Function Severity:  3 - Normal
Priority:  5 - Normal Item Group:  Feature Request
Status:  In Progress Assigned to:  None
Originator Name:  Axel Open/Closed:  * Open
Release:  * dev Operating System:  * Any
Fixed Release:  None Planned Release:  None
* Mandatory Fields

Add a New Comment Rich Markup
   

Jump to the original submission

Fri 02 Sep 2022 11:40:33 AM UTC, comment #14: 

I'm attaching the changeset that adds the ode23tb function from Fabio Cassini's repository here as a patch (to reduce the risk of losing it completely).

Maybe someone can pick it up from there.

I haven't done anything with it apart from quickly scrolling through it. I noticed that there are still some TODO comments. But it looks pretty advanced otherwise.


(file #53652)

Markus Mützel <mmuetzel>
Group administrator
Thu 01 Sep 2022 03:05:32 PM UTC, comment #13: 

I just found out that the initial implementation of ode23tb
that was done by Fabio Cassini did not go completely lost when
BitBucket removed mercurial support.

A backup archive may be found here :

https://bitbucket-archive.softwareheritage.org/projects/fa/FabioCassini/octave.html


Anonymous
Sun 29 Aug 2021 03:46:34 AM UTC, comment #12: 

I think it is not possible without required methods in octave or by using SUNDIALS.
 
           can anyone suggest me any other possible way to work on it?

kumar <introcoding>
Sat 28 Aug 2021 06:37:47 AM UTC, comment #11: 

Don't use Matlab's code. Don't even look at it. Their license doesn't allow taking any part or "inspiration" from it.
Nevertheless, the interface should be compatible. Afaict, it is ok to read their documentation of the function or ask here or on Octave's discourse for running some tests on Matlab's function to see its output for specific input parameters.

Markus Mützel <mmuetzel>
Group administrator
Sat 28 Aug 2021 06:22:37 AM UTC, comment #10: 

So there's a matlab code on ode23tb,we may follow in that manner?

kumar <introcoding>
Tue 24 Aug 2021 03:14:16 PM UTC, comment #9: 

From their website:
https://de.mathworks.com/help/matlab/ref/ode23tb.html

> ode23tb is an implementation of TR-BDF2, an implicit Runge-Kutta formula with a trapezoidal rule step as its first stage and a backward differentiation formula of order two as its second stage. By construction, the same iteration matrix is used in evaluating both stages. Like ode23s and ode23t, this solver may be more efficient than ode15s for problems with crude tolerances [1], [2].


I guess that means that the TR-BDF2 algorithm would be at the heart of the implementation.

Preferably, the interface to the new function would be compatible to Matlab's `ode23tb`. So, this interface should be implemented additionally to the actual algorithm.

Markus Mützel <mmuetzel>
Group administrator
Tue 24 Aug 2021 09:47:13 AM UTC, comment #8: 

I think we've to implement the code only related to tr-bdf2?

kumar <introcoding>
Fri 16 Jul 2021 04:40:13 PM UTC, comment #7: 

Yeah tq,i'll try the same

kumar <introcoding>
Fri 16 Jul 2021 03:43:36 PM UTC, comment #6: 

Lawrence F. Shampine has a web page on some interesting info:

https://faculty.smu.edu/shampine/current.html

I assume he would not mind sending a copy of his paper if asked politely.

Dmitri.
--



Dmitri A. Sergatskov <dasergatskov>
Fri 16 Jul 2021 08:34:51 AM UTC, comment #5: 

According to the Matlab documentation, they are using an algorithm described in these papers:
https://de.mathworks.com/help/matlab/ref/ode23tb.html

> [1] Bank, R. E., W. C. Coughran, Jr., W. Fichtner, E. Grosse, D. Rose, and R. Smith, “Transient Simulation of Silicon Devices and Circuits,” IEEE Trans. CAD, 4 (1985), pp. 436–451.
> [2] Shampine, L. F. and M. E. Hosea, “Analysis and Implementation of TR-BDF2,” Applied Numerical Mathematics 20, 1996.


I only found those behind paywalls though.

Markus Mützel <mmuetzel>
Group administrator
Fri 16 Jul 2021 03:29:38 AM UTC, comment #4: 


ohh sorry,i thought it as same as ode15.You know any packages/software to integrate it,i'm working on it's code.

kumar <introcoding>
Wed 14 Jul 2021 07:52:09 AM UTC, comment #3: 


comment #2:

> I want to do this by using sundial ode packages.Can you assign to me?



AFAIK sundials does not offer this feature so what you propose does not sound feasible.

Carlo de Falco <cdf>
Group Member
Wed 14 Jul 2021 07:43:11 AM UTC, comment #2: 

I want to do this by using sundial ode packages.Can you assign to me?

kumar <introcoding>
Mon 14 Aug 2017 02:41:44 PM UTC, comment #1: 

This is in progress, but will be part of Octave core,
not of odepkg as ode23tb is also a core function in Matlab.

Carlo de Falco <cdf>
Group Member
Tue 29 Sep 2015 08:23:25 AM UTC, original submission:  

Dear all,

could someone implement "ode23tb" solver in octave, still missing in odepkg?

Thank you.

Axel

Axel <axelt83>

 

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

Attach Files:
   
   
Comment:
   

Attached Files
file #53652:  bug46074-ode23tb.patch added by mmuetzel (41KiB - text/plain)

 

Depends on the following items: None found

Items that depend on this one: None found

 

Carbon-Copy List
  • -email is unavailable- added by dasergatskov (Posted a comment)
  • -email is unavailable- added by mmuetzel (Posted a comment)
  • -email is unavailable- added by introcoding (Posted a comment)
  • -email is unavailable- added by cdf (Posted a comment)
  • -email is unavailable- added by cdf
  • -email is unavailable- added by cdf
  • -email is unavailable- added by cdf
  • -email is unavailable- added by jwe (Updated the item)
  • -email is unavailable- added by axelt83 (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 12 latest changes.

    Date Changed by Updated Field Previous Value => Replaced by
    2022-09-02 mmuetzel Attached File- Added bug46074-ode23tb.patch, #53652
    2017-08-14 jwe CategoryOctave Package Octave Function
        StatusNone In Progress
        Release4.0.0 dev
        Operating SystemMicrosoft Windows Any
        Summary[octave forge] (odepkg) missing function ode23tb missing function ode23tb
    2017-08-14 cdf Carbon-Copy- Added -email is unavailable-
        Carbon-Copy- Added -email is unavailable-
        Carbon-Copy- Added -email is unavailable-
    2017-08-13 jwe Summary[odepkg] ode23tb [octave forge] (odepkg) missing function ode23tb
    2015-09-30 rik5 Carbon-CopyRemoved 72865 -
    2015-09-30 rik5 CategoryOctave Function Octave Package

    Back to the top

    Powered by Savane 3.13-3230.
    Corresponding source code