What's the reason for this change?  this->n_reactions() calls _reaction_set.n_reactions(), and _reaction_set just got initialized from reaction_set... so there should be no difference, right?