search

wgreene310 / pde1dm

1D PDE Solver for MATLAB and Octave written in m-code
GNU General Public License v3.0
7 stars 1 forks source link

Solving DGLs with complex numbers seems not to work #8

Open MQwertzy opened 1 year ago

MQwertzy commented 1 year ago

dde1dm seems to abort the calculation under octave as soon as values in my DGLs become complex with the following error: error: ode15i: "Jacobian" function must evaluate to a real square matrix error: called from ode15i at line 207 column 11 solveTransient at line 253 column 21 pde1dm at line 124 column 7 I can avoid the program abort in my case if I require s=real(s); in the pdeFunc definition. Unfortunately this does not work for my mathematical problem. Using Matlab and replacing pde1dm with pdepe my DGLs can be solved without problems. Is this intended/expected behavior?

wgreene310 commented 1 year ago

Unfortunately this is due to a limitation of the ode15i implementation in octave. I added a sentence to the pde1dm manual about this.

pde1dm in Matlab should give the same solution as pdepe to your complex equations. Please let me know if this is not the case.