Merge pull request ComputationalRadiationPhysics#891 from psychocoderHPC/topic-cleanDSNotTypeSaveImpl

clean field solver DS float usage

Author: ax3l

src/picongpu/include/fields/MaxwellSolver/DirSplitting/DirSplitting.kernel

@@ -42,17 +42,17 @@ struct DirSplittingKernel
     template<typename CursorE, typename CursorB >
     DINLINE void propagate(CursorE cursorE, CursorB cursorB) const
     {
-        float a_plus = (*cursorB(-1, 0, 0)).z() + (*cursorE(-1, 0, 0)).y();
-        float a_minus = (*cursorB(1, 0, 0)).z() - (*cursorE(1, 0, 0)).y();
-        float a_prime_plus = (*cursorB(-1, 0, 0)).y() - (*cursorE(-1, 0, 0)).z();
-        float a_prime_minus = (*cursorB(1, 0, 0)).y() + (*cursorE(1, 0, 0)).z();
+        float_X a_plus = (*cursorB(-1, 0, 0)).z() + (*cursorE(-1, 0, 0)).y();
+        float_X a_minus = (*cursorB(1, 0, 0)).z() - (*cursorE(1, 0, 0)).y();
+        float_X a_prime_plus = (*cursorB(-1, 0, 0)).y() - (*cursorE(-1, 0, 0)).z();
+        float_X a_prime_minus = (*cursorB(1, 0, 0)).y() + (*cursorE(1, 0, 0)).z();
 
         __syncthreads();
 
-        (*cursorB).z() = 0.5f * (a_plus + a_minus);
-        (*cursorE).y() = 0.5f * (a_plus - a_minus);
-        (*cursorB).y() = 0.5f * (a_prime_plus + a_prime_minus);
-        (*cursorE).z() = 0.5f * (a_prime_minus - a_prime_plus);
+        (*cursorB).z() = float_X(0.5) * (a_plus + a_minus);
+        (*cursorE).y() = float_X(0.5) * (a_plus - a_minus);
+        (*cursorB).y() = float_X(0.5) * (a_prime_plus + a_prime_minus);
+        (*cursorE).z() = float_X(0.5) * (a_prime_minus - a_prime_plus);
 
         __syncthreads();
     }