Bug #1492

LogParabola flux method returns nan

Added by Buehler Rolf almost 9 years ago. Updated over 8 years ago.

Status:ClosedStart date:06/30/2015
Priority:NormalDue date:
Assigned To:-% Done:

0%

Category:-
Target version:-
Duration:

Description

Running on gammalib v0.10.0 on Ubuntu 14.04 the following lines:

import gammalib as gl
crab = gl.GModels("/home/buehler/Code/gammalib/test/data/model_point_logparabola.xml”)
crab[“Crab"].spectral().flux(gl.GEnergy(1,“TeV"),gl.GEnergy(10,“TeV”))

results in:

+++ WARNING in GIntegral::romberg(1e+06, 1e+07, 5): Integration uncertainty nan exceeds absolute tolerance of nan after 21 iterations. Result -nan is inaccurate.

This was very likely the cause of the problem I reported in ctools Issue #1447 https://cta-redmine.irap.omp.eu/issues/1447


Recurrence

No recurrence.

History

#1 Updated by Knödlseder Jürgen over 8 years ago

Can you check whether the problem still exists on the actual code?

Michael has reported a similar problem (#1556) which was traced back to the problem of invalid references. Michael’s problem has been fixed.

#2 Updated by Buehler Rolf over 8 years ago

Running on Version: 1.1.0.dev (22 November 2015) the problem is gone (no warning appears and the returned value looks reasonable)

#3 Updated by Knödlseder Jürgen over 8 years ago

  • Status changed from New to Closed

Thus close the issue.

Also available in: Atom PDF