Resolve "mixed-precision"

Closes #146 (closed)

src/Solver/FFTPoissonSolver.h

@@ -33,7 +33,7 @@ namespace ippl {
         // types for LHS and RHS
         using lhs_type = typename Solver<Tlhs, Trhs, Dim, Mesh, Centering>::lhs_type;
         using rhs_type = typename Solver<Tlhs, Trhs, Dim, Mesh, Centering>::rhs_type;
-
+	using Tg = typename Tlhs::value_type;
         // type of output
         using Base = Electrostatics<Tlhs, Trhs, Dim, Mesh, Centering>;
 
@@ -41,8 +41,10 @@ namespace ippl {
         // define a type of Field with integers to be used for the helper Green's function
         // also define a type for the Fourier transformed complex valued fields
         typedef Field<Trhs, Dim, Mesh, Centering> Field_t;
+        typedef Field<Tg, Dim, Mesh, Centering> Field_gt;
         typedef Field<int, Dim, Mesh, Centering> IField_t;
         typedef Field<Kokkos::complex<Trhs>, Dim, Mesh, Centering> CxField_t;
+        typedef Field<Kokkos::complex<Tg>, Dim, Mesh, Centering> CxField_gt;
         typedef Vector<Trhs, Dim> Vector_t;
 
         // define type for field layout
@@ -75,7 +77,7 @@ namespace ippl {
         void initializeFields();
 
         // communication used for multi-rank Vico-Greengard's Green's function
-        void communicateVico(Vector<int, Dim> size, typename CxField_t::view_type view_g,
+        void communicateVico(Vector<int, Dim> size, typename CxField_gt::view_type view_g,
                              const ippl::NDIndex<Dim> ldom_g, const int nghost_g,
                              typename Field_t::view_type view, const ippl::NDIndex<Dim> ldom,
                              const int nghost);
@@ -131,9 +133,9 @@ namespace ippl {
         std::string alg_m;
 
         // members for Vico-Greengard
-        CxField_t grnL_m;
+        CxField_gt grnL_m;
 
-        std::unique_ptr<FFT<CCTransform, Dim, double, Mesh, Centering>> fft4n_m;
+        std::unique_ptr<FFT<CCTransform, Dim, Tg, Mesh, Centering>> fft4n_m;
 
         std::unique_ptr<Mesh> mesh4_m;
         std::unique_ptr<FieldLayout_t> layout4_m;
@@ -145,7 +147,6 @@ namespace ippl {
 
         // buffer for communication
         detail::FieldBufferData<Trhs> fd_m;
-
     protected:
         virtual void setDefaultParameters() override {
             using heffteBackend       = typename FFT_t::heffteBackend;