[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[Octave-patch-tracker] [patch #9063] Add function decic
From: |
Bill Greene |
Subject: |
[Octave-patch-tracker] [patch #9063] Add function decic |
Date: |
Mon, 19 Sep 2016 16:57:19 +0000 (UTC) |
User-agent: |
Mozilla/5.0 (Windows NT 6.1; WOW64) AppleWebKit/537.36 (KHTML, like Gecko) Chrome/53.0.2785.116 Safari/537.36 |
Follow-up Comment #5, patch #9063 (project octave):
I do not know for certain that the paper describes exactly
what is used in the MATLAB decic.
But, as Shampine points out in the paper, it does have one
very attractive feature: if the problem is not a DAE
(i.e. dF/dyp has full rank), the algorithm simply solves
for yp, leaving the user's initial conditions in y completely
unchanged.
I assume you could do the same thing with the algorithm here by fixing all the
y-components and freeing all the yp components. But the user would have to
figure out that is what they need to do.
In the general DAE case, the Shampine algorithm does not require the user to
identify which equations are odes and which are algebraic.
_______________________________________________________
Reply to this item at:
<http://savannah.gnu.org/patch/?9063>
_______________________________________________
Message sent via/by Savannah
http://savannah.gnu.org/