Problem 1531 - precison problem in G4PhysicsVector
Summary: precison problem in G4PhysicsVector
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: global/management (show other problems)
Version: 9.6
Hardware: PC Windows
: P5 normal
Assignee: Vladimir.Ivantchenko
URL:
Depends on:
Blocks:
 
Reported: 2013-10-28 10:37 CET by Pedro.Arce
Modified: 2013-12-13 10:14 CET (History)
0 users

See Also:

Attachments
Add an attachment (proposed patch, testcase, etc.)

Note You need to log in before you can comment on or make changes to this problem.
Description Pedro.Arce 2013-10-28 10:37:09 CET 
I found a precision problem when running geant4.9.6.p02 in windows (Visual C++10). I have an energy 1 MeV that indeed is 0.9999999999999989. Then at G4PhysicsVector::Interpolate

  } else if( theEnergy >= edgeMax ) {

says it is smaller, but at 

FindBinLocation it is located in the last bin, and when an interpolation is done between this bin and bin+1 it crashes because it is out of range.
Comment 1 Vladimir.Ivantchenko 2013-12-12 17:24:17 CET 
The problem has been fixed in the new 10.0 release. The fix will be back-ported to 9.6 in the time of the new patch.