25 Mar 15:00
possible integration problem with v5.1 and v5.2
Frederic Y. Bois <frederic.bois <at> ineris.fr>
2008-03-25 14:00:20 GMT
2008-03-25 14:00:20 GMT
You may encounter problems when integrating a set of differential equations: with a very dense scheduling of output or input function discontinuities, the integrator may get stuck at a stopping point and does not advance past it. This seems to happen only when the code is compiled with the optimization option on (at least with gcc), which is by default when you use the Makefiles provided. You can try running an non-optimized executable. The simplest way to get at that is, for example: cd ~/mcsim/sim mod mymodel.model gcc -o mymodel.exe *.c -lgsl -lgslcblas -lm for optimized code you would do gcc -O -o mymodel.exe *.c -lgsl -lgslcblas -lm This assumes that "mod", the model translation software is on you command path. In any case, I am trying to understand and fix the problem. It is not a sneaky bug that would produce wrong results silently. The failure is reported and the program aborts immediately, the problem does not go unnoticed. F. Bois
RSS Feed