Ticket #26 (closed defect - FPE: fixed)

Opened 1 year ago

Last modified 1 year ago

Arithmetic exception in H2_LevelPops()

Reported by: will Assigned to: nobody
Priority: critical Milestone:
Component: H chemistry (deprecated) Version:
Keywords: dynamics Cc: gary@pa.uky.edu, rjrwilliams@googlemail.com

Description

After several hours, the attached model falls over with an fpe like this: 

Starting program: /fs/pinga/other0/will/Cloudy/trunk/source/cloudy.exe < ../in/ flow6_n35big.in > flow6_n35big.out

Program received signal SIGFPE, Arithmetic exception.
0x0000000000460175 in H2_LevelPops () at mole_h2.cpp:2351
2351                                            SDIV( H2_Boltzmann[0][iVib][iRo t] ) );

Attachments

flow6_n35big.in (1.3 kB) - added by will on 05/11/07 11:48:53.

Change History

05/11/07 11:48:53 changed by will

  • attachment flow6_n35big.in added.

05/16/07 11:12:52 changed by will

Gary proposed a fix to the data in the file H2_coll_He_ORNL.dat. However, this does not resolve the problem. The fpe still occurs at exactly the same point.

05/26/07 09:57:03 changed by gary

debugger shows exception thrown with these particulars

crash mole_h2.cpp line 302, nColl 5 iVibHi = 8 iRotHi = 3 iVibLo = 8 iRotLo = 1

Te = 1201 collider density 7818

H2_CollRate = 1e40

H2_CollRate set line 1674 of mole_h2.cpp

05/28/07 10:25:54 changed by gary

  • status changed from new to closed.
  • resolution set to fixed.

fixed in rev 1163 - the H2-He collision interpolation equation has a pole within the temperature range needed. Phillip Stancil tested and proposed taking the absolute value of two of the fitting coefficients so that no pole is possible.