Problem 1251 - FPE crash in G4QNucleus.cc at line 2416
Summary: FPE crash in G4QNucleus.cc at line 2416
Status: RESOLVED FIXED
Alias: None
Product: Geant4
Classification: Unclassified
Component: processes/hadronic/models/chiral_inv_phase_space (show other problems)
Version: 9.4
Hardware: All All
: P5 trivial
Assignee: Mikhail.Kossov
URL:
Depends on:
Blocks:
 
Reported: 2011-10-04 11:10 CEST by Andrea Dotti
Modified: 2011-10-05 12:10 CEST (History)
1 user (show)

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 Andrea Dotti 2011-10-04 11:10:05 CEST 
An FPE crash was issued at line 2416 of G4QNucleus in ref-08.
This was due to an unprotected sqrt of a possible negative number.

The code was:
if ( evalph && aFlag && mAlph+GSResNa<totMass) {
  G4double ken = totMass - mAlph - GSResNa;
  if (barf) ken -= ABarr;
  .... sqrt(ken) ....
}

The code has been modified (following the example of lline 2398:
if ( evalph && aFlag && mAlph+GSResNa+ABarr<totMass)

Modification being tested.

Andrea
Comment 1 Andrea Dotti 2011-10-04 11:19:05 CEST 
Tests are positive.
Tag hadr-chips-body-V09-04-02 proposed for testing.