IDA_NO_RECOVERY

[Jolandi Labuschagne]

Hi,

I am using IDA solver and get the error copied below. From what I've seen this might be due to the DAE index not being equal to 1. Is there an easy way to check the DAE index for a large system of equations? I feel confident in my initial conditions so I don't believe they might be causing this problem.

-------------------------------------------------------------------------------

psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:852: Linear solve failed psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:852: Linear solve failed psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:852: Linear solve failed psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:852: Linear solve failed psetup failed: .../casadi/interfaces/sundials/idas_interface.cpp:852: Linear solve failed The residual routine or the linear setup or solve routine had a recoverable error, but IDACalcIC was unable to recover.

--------------------------------------------------------------------------- RuntimeError ​Traceback (most recent call last) Input In [260], in <module> ----> 1 ​r = F(x0=x0, z0=z0, p=p0) File ~/miniconda3/envs/research.flotationmaterialmodel/lib/python3.8/site-packages/casadi/casadi.py:13453, in Function.__call__(self, *args, **kwargs) ​13450 return tuple(ret) ​13451 else: ​13452 # Named inputs -> return dictionary > 13453 ​return self.call(kwargs) File ~/miniconda3/envs/research.flotationmaterialmodel/lib/python3.8/site-packages/casadi/casadi.py:12324, in Function.call(self, *args) ​12256 def call(self, *args) -> "void": ​12257 """ ​12258 Generate a Jacobian function of output oind with respect to input iind. ​12259 ​(...) ​12322 ​12323 """ > 12324 ​return _casadi.Function_call(self, *args) RuntimeError: .../casadi/interfaces/sundials/idas_interface.cpp:591: IDACalcIC returned