High Index DAE Solution

APMonitor and Gekko solve high index Differential and Algebraic Equations (DAEs) without rearrangement or differentiation. A common method for high index DAEs is to use Pantelides algorithm to reduce the index order. The index of a DAE is the number of times the equations must be differentiated to return to index-0 (ODE) form. Most ODE and DAE solvers have several challenges:
1. Consistent initial condition is required
2. Required index-1 or index-2 (Hessenberg form)
3. Numerical drift with repeated differentiation
APMonitor and Gekko overcome these challenges by solving high index DAEs without differentiation with orthogonal collocation on finite elements and efficient sparse solvers.
See apmonitor.com/wiki/index.php/... for source code.

[@HuyNguyen-bw4sv]

Thank you!

[@YouR.Dumbass]

Dear Prof. Hedengren,
I have the following error while running Index-1, whereas other files runs just fine.
Exception Traceback (most recent call last)
Cell In[2], line 26
24 m.options.IMODE=4
25 m.options.NODES=3
---> 26 m.solve(disp=False)
28 plt.figure(figsize=(10,5))
29 plt.subplot(3,1,1)
File c:\Users\human\anaconda3\Lib\site-packages\gekko\gekko.py:2185, in GEKKO.solve(self, disp, debug, GUI, **kwargs)
2183 #print APM error message and die
2184 if (debug >= 1) and ('@error' in response):
-> 2185 raise Exception(response)
2187 #load results
2188 def byte2str(byte):
Exception: @error: Solution Not Found

[@apm]

Thanks for letting me know about that error. Could you create a question on StackOverflow with the code you are using?