Hi everybody.
I have been working with a model where the equations of motion for the background can't be solved numerically in CAMB, but they can be solved using mathematica or python. Briefly the situation is the following: I have a system of three ODE y'(1), y'(2), y'(3) whose initial conditions are set by theoretical considerations. Both y(2) and y(3) undergo an initial sudden drop (of eight and 23 orders of magnitude, respectively) and then they evolve smoothly. Whereas this situation does not pose any problem in mathematica or python, I have found that dverk is unable to handle this drop and for that reason the program crashes.
So my question is if there are other ODE solvers that can be incorporated to CAMB, besides dverk. If so, any hint or a link where somebody has done something like that would be appreciated.
Regards
a dverk performance issue
-
- Posts: 29
- Joined: June 10 2012
- Affiliation: University of Mexico
-
- Posts: 29
- Joined: June 10 2012
- Affiliation: University of Mexico
a dverk performance issue
Hello everbody
Just in case somebody is interested on this topic, I could embed a python script into CAMB. The way I proceeded was to use the instruction
where arg are the input arguments of the python script. So far the program has run perfectly.
Regads
Just in case somebody is interested on this topic, I could embed a python script into CAMB. The way I proceeded was to use the instruction
Code: Select all
call EXECUTE_COMMAND_LINE ("my_python_script arg1 arg2 arg", wait=.true.)
Regads