Changeset 1340 for trunk/source/global/management/src/G4Pow.cc

Timestamp:

Nov 5, 2010, 3:45:55 PM (14 years ago)

Author:

garnier

Message:

update ti head

File:

• trunk/source/global/management/src/G4Pow.cc (modified) (4 diffs)

trunk/source/global/management/src/G4Pow.cc

@@ -24 +24 @@
 // ********************************************************************
 //
-// $Id: G4Pow.cc,v 1.2 2009/07/03 11:35:06 vnivanch Exp $
-// GEANT4 tag $Name: geant4-09-04-beta-01 $
+// $Id: G4Pow.cc,v 1.4 2010/08/24 08:12:08 gcosmo Exp $
+// GEANT4 tag $Name: global-V09-03-22 $
 //
 // -------------------------------------------------------------------
@@ -65 +65 @@
   lz.resize(maxZ,0.0);
   fact.resize(maxZ,0.0);
+  logfact.resize(maxZ,0.0);
 
   onethird = 1.0/3.0;
   G4double f = 1.0;
 
-  for(G4int i=1; i<maxZ; i++)
+  for(G4int i=1; i<maxZ; ++i)
   {
     G4double x  = G4double(i);
@@ -76 +77 @@
     f      *= x;
     fact[i] = f;
+    logfact[i] = std::log(f);
   }
   fact[0] = 1.0;
+  logfact[0] = 0.0;
 }
 
@@ -86 +89 @@
 
 // -------------------------------------------------------------------
+
+G4double G4Pow::powN(G4double x, G4int n)
+{
+  G4double res = 1.0;
+  if(n >= 0) { for(G4int i=0; i<n; ++i) { res *= x; } }
+  else if((n < 0) && (x != 0.0))
+  {
+    G4double y = 1.0/x;
+    G4int nn = -n;
+    for(G4int i=0; i<nn; ++i) { res *= y; }
+  }
+  return res;
+}