sPuReMD: refactor and clean up preconditioner application code (level scheduling, graph coloring, Jacobi iteration). Conditionally change charge matrix for setting up preconditioner for ICHOLT and ILU(T).