I am trying to work with RT3D using a rnxs.dll file as a user-defined reaction function and it works fine. My systems of reactions are so stiff and I would like to define a jacobian function that works along with it. I follow the same steps of making and compiling the rxns.dll for jacrxns.dll file. But the result that I get by Isolver=2 (gear with explicit jacobian) is different from Isolver=1 (general gear solver).
I would like to know if the jacrxns.dll file should be compiled separately from rxns.dll or not?
If yes, should they be somehow linked together? Or there is no need for that?
Any comment on this topic is greatly appreciated.