bool newton_step(GroupMap &mole_map, const valarray< double > &oldmols, valarray< double > &newmols, realnum *eqerror, realnum *error, const long n, double *rlimit, double *rmax, valarray< double > &escale, void(*jacobn)(GroupMap &mole_map, const valarray< double > &b2vec, double *const ervals, double *const amat, const bool lgJac, bool *lgConserved))