Odeint check routine

Asked by Shivam Gola

Hi

I want to use odeint check routine for education perpouse, but it seems it is not working. I am using latest version of maddm and madgraph.

Question information

Language:
English Edit question
Status:
Answered
For:
MadDM Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Can you create a bug report and attach the debug file?
(or copy-paste it in this question)

Cheers,

Olivier

Revision history for this message
Shivam Gola (shivay12) said :
#2

Hi,

Many thanks for such quick response.

 I just add "call odeint_check(1)" routine in the maddm.f file from directory ../MG5_aMC_v2_6_7/alp_subroutins/src/.

Then i compile using make command and run it. There is no error message as well as no new file in output folder.

Thanks,
Shivam

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Ok,

I'm actually not an expert of that part of the code,
However the issue is likely that the information is not available at the time you put such lines.
And that you should have (in one of the log file) the following printed:

c checks to make sure that the ODE information has been saved
      if (dxsav.eq.0.d0) then
        write(*,*) '---------------------------------------------------------'
        write(*,*) 'odeint_check:'
        write(*,*) ' ODE information not saved!!!'
        write(*,*) ' dxsav, kount, and kmax need to be set in relicdensity.f'
        write(*,*) '---------------------------------------------------------'
        return
      endif

Cheers,

Olivier

Revision history for this message
Shivam Gola (shivay12) said :
#4

Sorry, I am bit ignorant. The above command line structure is written in 'tests.f' file in the odeint subroutine. The file 'tests.f' is in same 'src' folder.

I really need bit more direction in this regard.

Thanks,
Shivam

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

Hi,

Not sure what I can do to help you to be honest, the expert on that parrt of the code left physics a while ago.
I guess that you know that part of the code as well as me.
And I'm not sure what you try to do/achieve on top of that.

Sorry,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Shivam Gola for more information if necessary.

To post a message you must log in.