-
Notifications
You must be signed in to change notification settings - Fork 11
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Errors reading *.eqn file in version 3.0 #78
Comments
@jnavarla, thanks for writing. Can you also post any other files like *.kpp and *.spc? I will try to replicate the issue. |
in your code? |
Thanks. Yes, here you have the kpp and spc files |
Hello Rolf, |
Yes, it should be removed. However, I wonder which integrator are you using, and where have you defined it? |
I was able to get KPP to work with: #INTEGRATOR rosenbrock
#LANGUAGE C
#DRIVER general
#JACOBIAN SPARSE_LU_ROW
#HESSIAN on
#STOICMAT on
#FUNCTION AGGREGATE
#DOUBLE ON
#include aquMTv2.spc
#include aquMTv2.eqn I believe that when you specify the |
I was able to get KPP to work with: #INTEGRATOR rosenbrock
#LANGUAGE C
#DRIVER general
#JACOBIAN SPARSE_LU_ROW
#HESSIAN on
#STOICMAT on
#FUNCTION AGGREGATE
#DOUBLE ON
#include aqu.spc
#include aqu.eqn I believe that when you specify the |
@yantosca: Maybe we should mention this in the box of our documentation at: Something like: The dummy integrator "none" does not exist anymore. Thus, commands like
should be removed. |
@RolfSander, I'll work on that. Thanks! |
Thx, @yantosca ! |
I'm using Rosenbrock.m defined in $KPP_HOME/int. Now, removing #INTEGRATOR none , KPP works with my *.eqn files and generates the matlab code, but it gives errors when running. Maybe some call is written in a different way in version 3.0 and gives the error in matlab; I'm running R2021a. I will check the generated matlab code to look at. |
Thank you @yantosca for your reply. Yes, I'm using a *.def file (attached here), but the problem was on #INTEGRATOR command as @RolfSander suggested. The matlab code gererated with kpp v2.2.3 works without problems, but the code generated with kpp v3.0.1 doens't. |
Sorry but I'm still trying to understand your code... You say that you are using Rosenbrock.m. This means that somewhere you
However, I cannot find these lines anywhere in the code that you have |
Sorry @RolfSander and @yantosca . I realized that I send you a wrong *.kpp file. I'm sorry to waste your time. Instead #INTEGRATOR rosenbrock, I'm using a #DRIVER rosendrv, With these files I can have a matlab code but I have execution errors related to Aout I read something about this problem in another matlab thread here and I will check that later. |
Thanks @jnavarla. In KPP 3.0.0 we added Aout as an argument so as to get the reaction rates out for diagnostics (if so desired). Also @RolfSander and I have mostly worked with the F90 and C interfaces but not with matlab, so there may still be a few issues. We had @jenniethompson and her group help us out with the matlab version for 3.0.0. |
Thanks @RolfSander . Yes, the problem in *_Fun.m generation still persist in kpp v3.0.1. As it is explained in #56 the arguments in the function call are wrong. |
Does it work if you copy the matlab files from the old KPP version into KPP-3.0.0? As @yantosca already mentioned, we are not actively using matlab. However, if you let us know what works for you, we can update the matlab files accordingly for the next KPP release... |
Yes, if you copy the ROOT_Fun.m matlab file generated from a previous version ok kpp, for instance in my case v2.2.3, changing the file generated by v3.0.1, the matlab code works properly without execution errors. I've tested this with my model, but the small_strato.* either. This model is small and it is easy to check the rate equations by hand. Both functions generate the same individual reaction rate terms and the time derivative of the state variables of the system (Vdot), |
Great to see that it works for you if you use ROOT_Fun.m from v2.2.3. Still, it would be nice to solve this problem in a better way. Could you post the diff between both ROOT_Fun.m files? That should give us a hint... |
@jimmielin had posted this issue #56 about issues with Fun in matlab. Not sure if this is why the behavior was changed in KPP 3.0.0. |
|
@RolfSander @jnavarla @jimmielin: Wondering if we can close out this issue. I believe #99 may have addressed this. |
@yantosca: Thx for taking care of those old issues. Yes, I think we can close this one. |
Closing this issue. Thanks @RolfSander |
Hi all,
I'm using KPP to develop an Ultrasound driven chemical reaction and I generated the model in KPP v2.2.
I have updated to KPP v3.0, but I have errors reading the *.eqn file looking for a 'none.def' file which is no defined.
I have reviewed all the files, but I cannot seek any difference when I process the files with v2.2.
Is there any care/suggestion for the eqn files in the new KPP version?
I attach a copy of *.eqn file
aquMTv2.eqn.txt
Thanks in advance for any suggestion.
Javier
The text was updated successfully, but these errors were encountered: