Resolve "mixed-precision"

Closes #146 (closed)

src/Solver/P3MSolver.hpp

@@ -250,7 +250,7 @@ namespace ippl {
         // for the P3M collision modelling method, it indicates
         // the splitting between Particle-Particle interactions
         // and the Particle-Mesh computations).
-        double alpha = 1e6;
+        Trhs alpha = 1e6;
 
         // calculate square of the mesh spacing for each dimension
         Vector_t hrsq(hr_m * hr_m);
@@ -264,7 +264,7 @@ namespace ippl {
         const int nghost                 = grn_m.getNghost();
         const auto& ldom                 = layout_mp->getLocalNDIndex();
 
-        constexpr double ke = 2.532638e8;
+        constexpr Trhs ke = 2.532638e8;
 
         // Kokkos parallel for loop to find (0,0,0) point and regularize
         Kokkos::parallel_for(
@@ -277,7 +277,7 @@ namespace ippl {
 
                 const bool isOrig = (ig == 0 && jg == 0 && kg == 0);
 
-                double r      = Kokkos::real(Kokkos::sqrt(view(i, j, k)));
+                Trhs r      = Kokkos::real(Kokkos::sqrt(view(i, j, k)));
                 view(i, j, k) = (!isOrig) * ke * (Kokkos::erf(alpha * r) / r);
             });