PG-PuReMD: fix bugs in the SAI preconditioner for the charge solver (memory leak, incorrect row numbers).

PG-PuReMD: port SAI preconditioner from MPI codebase. Rework charge solver preconditioner code. Other code clean-up.

PG-PuReMD: fix cross-stream race conditions on intermediary force calculation variables (CdDelta, Cdbo, etc.). Rework valency and torsion calculations to minmize the number of atomic operations performed. Split lone pair and over-/under-coordination kernels. Fix small size mismatch in allocation routines (over-allocation). Other code clean-up.

sPuReMD, PG-PuReMD: add control file keyword (include_polarization_energy) to toggle polarization energy calcuation in ReaxFF.

PG-PuReMD: fix warnings (and errors) from several uninitialized memory regions due to multiple stream change. Fix race condition with valence angle computation due to mistake with dynamic shared memory. Add control file parameter for setting the number of CUDA streams. Other code clean-up.

PG-PuReMD: fix issue with post-processing step with QEq pseudo-charges (too little storage in q array). Fix issue with device-wide CUB reductions and prefix sums reallocating device memory too frequently. Rework CUDA-aware MPI code to fix issues (still seeing out-of-memory errors on Ascent cluster at OLCF).

PG-PuReMD: ensure MPI pack/unpack CUDA kernels are run in the correct stream. Other changes to bring the CUDA-aware MPI functions for the charge solver in line with their MPI counterparts.

PG-PuReMD: utlize multiple CUDA streams to run independent kernels in the ReaxFF+QEq task dependency graph.

PG-PuReMD: fix stream issue with cudaMemset's happening in default stream (switch to cudaMemsetAsync). Rework memory management in light of this change.

WIP: finish the explicit single stream port by using cudaMemcpyAsync and cudaStreamSynchronize calls.

WIP: changes to support multiple CUDA streams. This commit contains porting of all kernels from the default stream to a single non-default stream. This commit has bugs regarding cudaMemcpy's still being in the default stream (and synchronous).

PG-PuReMD: implement MPI+CUDA versions of SDM, PIPECG, PIPECR for charge solver methods (along with missing dual QEq versions).

PG-PuReMD: use a warp of threads for addtional bonded interaction computations (bond orders, bonds, lone pair/over/under).

PG-PuReMD: use a warp of threads for apply bond derivatives to forces and simplify these code paths.

PG-PuReMD: fix incorrect summations in kinetic energy and temperature calculation functions. Replace several hand-coded device-wide reduction kernels with CUB libary reduction.

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).

sPuReMD: make zlib dependency optional (trajectory file compression). PG-PuReMD: remove erroneous zlib header file inclusion.

PG-PuReMD: handle a few missed reductions in nonbonded interactions more efficiently. Clean-up bonded interactions functions and do not export kernels via header files. Increase C++ standard to c++14 for CUB (Thrust dependency). Continue with bundled CUB version and transition to version included with CUDA SDK in the future.

PG-PuReMD: use atomics for both partial energy and force accumulation (small reworks to reduce number of atomic operations). Remove several unneccesary cudaMemset calls. Other code clean-up.

Test: fix Python wrapper for QM/MM (interfaces out-of-date). Fix Fortran tester bugs. Add test artifacts.