 
 // Constant memory for gridpoint weighting
 #define CONSTANT_SIZE 2048
+//! The developer has chosen not to document this variable
 __device__ __constant__ float GPU_rho_coeff[CONSTANT_SIZE];
-
-
-// global variables for thermodynamic output
-int g_Nx, g_Ny, g_Nz, g_block_size;
-cufftComplex *g_rho_real_space;
-float3 *g_vg;
-float *g_green_hat;
-float2 *o_data, *idat;
 
 /*! \file PPPMForceGPU.cu
   \brief Defines GPU kernel code for calculating the Fourier space forces for the Coulomb interaction. Used by PPPMForceComputeGPU.


}


float float_reduce(float* i_data, float* o_data, int n) {

    float gpu_result = 0.0;

    return cudaSuccess;
        }

//! The developer has chosen not to document this function
__global__ void calculate_thermo_quantities_kernel(cufftComplex* rho, 
                                                   float* green_function, 
                                                   float* energy_sum,
