PG-PuReMD: fix out-of-bounds mmeory accesses and unitialized data usage in energy and force tabulation routines for van der Waals and Coulomb interactions.

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: corrections to van der Waals and Coulomb energies calculations (far neighbor list in full format, not half). Fix QEq charge matrix initialization (full format). Other refactoring and 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 issue with portions of interaction lists not being initialized for Cuda_Init_Forces. Utilize more targeted debugging prints (DEBUG -> DEBUG_FOCUS). Other general code cleanup.

PG-PuReMD: fix initial list offsets for creating lists. Fix far neighbor and hydrogen bond reallocation checks. Begin refactoring function defitions and declarations to allow compiler to enforce intentions of pointer usage.

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: add deallocation routines at simulation termination for aiding with validating memory management.

PG-PuReMD: fix charge matrix allocation issue. Update CUDA code to reflect data structure changes. Begin removing global variables. Seperate struct declarations and definitions to avoid circular data dependencies when refactoring data structures. Small I/O updates.

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.

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.