Make ionization work after merge of ComputationalRadiationPhysics#594,C…

…omputationalRadiationPhysics#596,ComputationalRadiationPhysics#607

After the merge of the pull requests ComputationalRadiationPhysics#594, ComputationalRadiationPhysics#596 and ComputationalRadiationPhysics#607 functionality of the
ionization feature has been re-established. The example `LWFA_with_Ionization`
was reduced to its minimum divergence from the default parameter set.

examples/LWFA_with_Ionization/include/simulation_defines/param/componentsConfig.param

@@ -1,6 +1,6 @@
 /**
- * Copyright 2013-2014 Axel Huebl, Anton Helm, Rene Widera, Felix Schmitt, 
- *                      Marco Garten
+ * Copyright 2013-2014 Axel Huebl, Heiko Burau, Anton Helm, 
+ *                     Rene Widera, Richard Pausch, Marco Garten
  *
  * This file is part of PIConGPU.
  *
@@ -35,52 +35,43 @@ namespace simulation_starter = defaultPIConGPU;
  *  - laserGaussianBeam     : Gaussian beam (focusing)
  *  - laserPulseFrontTilt   : Gaussian beam with a tilted pulse envelope
  *                            in 'x' direction
- *  - laserWavepacket       : wavepacket (gaussian in time and space, not focusing)
+ *  - laserWavepacket       : wavepacket (Gaussian in time and space, not focusing)
  *  - laserPlaneWave        : a plane wave
  */
 namespace laserProfile = laserGaussianBeam;
 
 /*! Field Configuration --------------------------------------------------
  *  - fieldSolverYee : standard Yee solver
- *  - fieldSolverLehe: Num. Cherenkov free field solver in a choosen direction
+ *  - fieldSolverLehe: Num. Cherenkov free field solver in a chosen direction
  *  - fieldSolverDirSplitting: Sentoku's Directional Splitting Method
  *  - fieldSolverNone: disable the vacuum update of E and B
  *
  * * For development purposes: ---------------------------------------------
  *  - fieldSolverYeeNative : generic version of fieldSolverYee
- *    (need more shared memory per gpu and is slow)
+ *    (need more shared memory per GPU and is slow)
  */
-#ifndef PARAM_FIELDSOLVER
-#define PARAM_FIELDSOLVER fieldSolverYee
-#endif
-namespace fieldSolver = PARAM_FIELDSOLVER;
+namespace fieldSolver = fieldSolverYee;
 
-/*! Gas Configuration------------------------------------------------------
+/*! Gas Configuration-----------------------------------------------------
  * Select a normed gas density profile
- *  - gasGaussian : super - gaussian density profile
- *  - gasGaussianCloud : super - gaussian density profile with a center position
+ *  - gasGaussian : super - Gaussian density profile
+ *  - gasGaussianCloud : super - Gaussian density profile with a center position
  *  - gasLinExp : piecewise linear-exponential profile (linear first)
- *  - gasHomogenous : constant gas density with a certain length
- *  - gasSphereFlanks : constant sphere droplet with exponential decr envelope
+ *  - gasHomogeneous : constant gas density with a certain length
+ *  - gasSphereFlanks : constant sphere droplet with exponentially decreasing envelope
  *  - gasFreeFormula: use a custom formula (slower)
  *  - gasFromHdf5: load gas density from HDF5 file
  *  - gasNone : just stay with a vacuum
  */
-#ifndef PARAM_GASPROFILE
-#define PARAM_GASPROFILE gasSphereFlanks
-#endif
-namespace gasProfile = PARAM_GASPROFILE;
+namespace gasProfile = gasSphereFlanks;
 
 /*! Particle Initialization ----------------------------------------------
  *  - particleInitRandomPos  : distribute your particles uniformly random
  *  - particleInitQuietStart : distribute your particles lattice-like
  */
 namespace particleInit = particleInitRandomPos;
 
-
 /*enable (1) or disable (0) current calculation*/
-#ifndef ENABLE_CURRENT
 #define ENABLE_CURRENT 1
-#endif
 
 }

examples/LWFA_with_Ionization/include/simulation_defines/param/gasConfig.param

@@ -1,5 +1,6 @@
 /**
- * Copyright 2013-2014 Axel Huebl, Rene Widera, Felix Schmitt, Marco Garten
+ * Copyright 2013-2014 Axel Huebl, Heiko Burau, Rene Widera, Felix Schmitt
+ *                     Richard Pausch, Marco Garten
  *
  * This file is part of PIConGPU. 
  * 
@@ -19,7 +20,6 @@
  */ 
 
 
-
 #pragma once
 
 namespace picongpu
@@ -47,7 +47,7 @@ namespace picongpu
 
     namespace gasGaussian
     {
-        /** Gas Formular:
+        /** Gas Formula:
                     const float_X exponent = fabs((y - GAS_CENTER) / GAS_SIGMA);
                     const float_X density = __expf(GAS_FACTOR*__powf(exponent, GAS_POWER));
 
@@ -101,8 +101,8 @@ namespace picongpu
         namespace SI
         {
             /** Parameters for the exponential slope
-             *  For Y > GAS_Y_MAX_SI:
-             *    let Y' = Y - GAS_Y_MAX_SI
+             *  For Y > GAS_MAX_SI:
+             *    let Y' = Y - GAS_MAX_SI
              *    \rho = exp[ - Y' * D ]
              *         = element [0.0; 1.0]
              *  unit: 1/m
@@ -112,7 +112,7 @@ namespace picongpu
 
     }
 
-    namespace gasHomogenous
+    namespace gasHomogeneous
     {
         namespace SI
         {
@@ -135,7 +135,7 @@ namespace picongpu
          *       ..,x,..            x high (constant)
          * 
          *  Imagine a ball with a "hard core" and "soft flanks".
-         *                 (+1 for using antithetics in a description)
+         *                 (+1 for using an antithesis in a description)
          */
         namespace SI
         {
@@ -164,7 +164,7 @@ namespace picongpu
 
     namespace gasGaussianCloud
     {
-        /** Gas Formular:
+        /** Gas Formula:
                     const float_X exponent = fabs((y - GAS_CENTER) / GAS_SIGMA);
                     const float_X density = __expf(GAS_FACTOR*__powf(exponent, GAS_POWER));