sPuReMD: fix issue with far neighbor list being incorrectly generated on MD steps where the list was not being computed (issue with atomic distances being incorrectly computed for peroidic neighbors). Disable several debugging print statements (DEBUG -> DEBUG_FOCUS). Other code clean-up.

PuReMD: code refactoring to remove global variables (similar to sPuReMD). sPuReMD: minor refactoring to match naming conventions of PuReMD.

sPuReMD: fix issue with SMALL_BOX_SUPPORT incorrectly computing pairwise atomic distances for periodic self-images. Fix issue with integer vector sum function (used in bond order and torsion interactions).

sPuReMD: remove global variable (array of function pointers for printing bonded interactions during debugging).

PuReMD: clean-up of heuristic for dynamically determining when to recompute the preconditioner for the charge solver.

sPuReMD: improvements for building in debug mode (add assertion support for run-time checks). Fix issue with Verlet list recomputation logic being incorrect. Fixes to be string comparisons.

sPuReMD: fix issue with certain preconditioners not being initialized for runs using restart files. Fix issue with part of simulation box struct not being initialized for runs using restart files (box->min). Other code clean-up from debugging issues with SMALL_BOX_SUPPORT.

sPuReMD: add safety check for BGF input files (CRYSTX). Remove compiler warnings when enabling DEBUG_FOCUS.

sPuReMD: fix issues with using sizeof with dynamically allocated pointers. Change OpenMP schedule policies for linear algebra code. PuReMD: correct discrepancy with charge matrix construction. Enable dual solver for charges by default.

sPuReMD: refactoring to better align with PuReMD code. PuReMD: backport sPuReMD changes from validation efforts against Fortran ReaxFF code.

sPuReMD: fixes to SAI preconditioner to correctly apply threshold to negative entries and to handle EE/ACKS2 cases when a threshold value of 1.0 is chosen.

sPuReMD: fix bug with diagonal scaling during FG_ICHOL/FG_ILUT producing incorrect factorizations when diagonal entries were negative.

sPuReMD: complete fine-grained Cholesky implementation. TODO: fix numeric breakdowns with water 6540 atom system.

sPuReMD: workarounds for memory issues with ILUT factors. Implement ILU(0) and call when selecting ILUT with a droptol of 0.

sPuReMD: fix issue with hbonds list indices not being reset correctly. Small algebraic rearrangements of valence angle calculations to match Fortran ReaxFF code (no numeric changes). Fix issue with rvec_isZero being incorrect. Fix issue with debugging arrays for TEST_FORCES not having indicies initialized. Other clean-up work.

sPuReMD: revert force tabulation changes. Temporarily work around out-of-bound memory accesses in cubic spline interpolation routes (for force tabulation) by expanding array sizes. Revisit at a later date.