Changeset 1315 for trunk/source/global/management/include/G4Pow.hh

Timestamp:

Jun 18, 2010, 11:42:07 AM (14 years ago)

Author:

garnier

Message:

update geant4-09-04-beta-cand-01 interfaces-V09-03-09 vis-V09-03-08

File:

• trunk/source/global/management/include/G4Pow.hh (modified) (4 diffs)

trunk/source/global/management/include/G4Pow.hh

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4Pow.hh,v 1.2 2009/07/03 11:35:07 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-03 $
+// $Id: G4Pow.hh,v 1.5 2010/05/28 08:18:03 vnivanch Exp $
+// GEANT4 tag $Name: geant4-09-04-beta-cand-01 $
 //
 //
@@ -35 +35 @@
 //
 // Utility singleton class for the fast computation of log and pow
-// functions 
+// functions. Integer argument should in the interval 0-255, no
+// check is performed inside these methods for performance reasons.
+// Computations with double arguments are fast for the interval
+// 0.5-255.5, standard library is used in the opposite case
 
 // Author: Vladimir Ivanchenko
@@ -109 +112 @@
 inline G4double G4Pow::A13(G4double A)
 {
-  const G4int maxZ = 256; 
+  const G4double minA = 0.5000001; 
+  const G4double maxA = 255.5; 
 
   G4double res;
-  G4int i = G4int(A + 0.5);
-  if((i >= 1) && (i < maxZ))
+  if((A >= minA) && (A <= maxA))
   {
+    G4int i = G4int(A + 0.5);
     G4double x = (1.0 - A/G4double(i))*onethird;
     res = pz13[i]*(1.0 + x - x*x*(1.0 - 1.66666666*x));
@@ -144 +148 @@
 inline G4double G4Pow::logA(G4double A)
 {
-  const G4int maxZ = 256;
+  const G4double minA = 0.5000001; 
+  const G4double maxA = 255.5; 
 
   G4double res;
-  G4int i = G4int(A + 0.5);
-  if((i >= 1) && (i < maxZ))
+  if((A >= minA) && (A <= maxA))
   {
+    G4int i = G4int(A + 0.5);
     G4double x = 1.0 - A/G4double(i);
     res = lz[i] + x*(1.0 - 0.5*x + onethird*x*x);