PG-PuReMD: split bond and hbond initialization kernels for future optimizations. Cleanup CUDA code related to far neighbor list being in full format (remove unnecessary conditionals and arithmetic).

PG-PuReMD: fix performance logging code around GPU code. Rearrange header files to allow some preprocessor definitions to be defined via options in configure script. Rework Autoconf and Automake code to allow passing flags directly to nvcc compiler wrapper (for CUDA code). Enable C11 and C++11 standard targets during compilation.

PG-PuReMD: small fix for initializtion routine performance logging (timers not reset after outputting to log file).

PG-PuReMD: fix issue with charge solver preconditioner refactoring rate causes issues with reneighboring actions (preconditionering rate was previously coupled with reneighoring rate for SAI but this causes issues for Jacobi, etc.). Be more greedy with memory allocation sizes to decrease reallocation frequency (MPI buffers, etc.). More GPU code clean-up.

PG-PuReMD: fix MPI buffer allocations sizes. Ensure that nonblocking MPI messages have completed for each dimension before continuing. Rework reallocation checks in integration routines. Temporarily disable CUDA-aware MPI code paths (need to perform packing/unpacking first on device before handing off pointers). Other code clean-up.

PG-PuReMD: fix issue with divergent MPI_Reduce calls by multiple MPI processes (timing logging code). Add run-time MPI routine error checking. Fix issue with upper limit of hydrogen atoms allowed being hard-coded in GPU code (use dynamic memory allocation instead). Other code clean-up.

PG-PuReMD: clean-up MPI custom datatype initialization and error reporting. Avoid truncations in buffer size calculations. Other code clean-up.

PG-PuReMD: merge changes from sPuReMD and PuReMD for validation work for simulation accuracy. Small correction to lone pair energy calculation. PuReMD: code formatting changes to better align with sPuReMD. sPuReMD: fix issue with valence angle calculations not setting list indices under certain sitations (taking previous continue statement missed index assignment).

PG-PuReMD: fix charge extrapolation. Improve memory management of local portion of charge matrix and reset linear system size for each time step.

PG-PuReMD: begin merging preconditioner code from sPuReMD. Make data structures between sPuReMD and PG-PuReMD consistent. Update control file parsing.

PG-PuReMD: merge memory management and step retry logic from MPI+GPU code into MPI-only code. Fix and backport numerous bugs (void pointer arithmetic and dereferences, etc.). Update code to match new data structures. Begin merging changes from OpenMP code (sPuReMD) for preconditioning work.

PG-PuReMD: only perform device-side reductions on energy at requested update energy frequency. Fix trajectory file output. Add num. retries of main simulation loop to log file. Other refactoring.

PG-PuReMD: begin enhancements. Enable sim. step retries. Enable far neighbors list max per atom. Other general code clean-ups.

The MPI-Not-GPU version is able to run through all benchmarks with decent agreement with the other versions now. The test harness has been modified to allow testing between the three functional versions.
Also added a README file in the test_harness directory about how to use the harness. Finally, I have saved one copy of the output from the previous test run to sim_output_082516.txt, for future comparison.