void PairReaxC::init_style( )

{

  if (!atom->q_flag) error->all("Pair reax/c requires atom attribute q");

  firstwarn = 1;

  int iqeq;

  for (iqeq = 0; iqeq < modify->nfix; iqeq++)

  if (eflag || vflag) ev_setup(eflag,vflag);

  else evflag = vflag_fdotr = eflag_global = vflag_global = 0;

  if ((eflag_atom || vflag_atom) && firstwarn) {

    firstwarn = 0;

    if (comm->me == 0) 

      error->warning("Pair reax/c cannot yet compute "

		     "per-atom energy or stress");

  }

  if (vflag_global) control->virial = 1;

  else control->virial = 0;

  class FixReaxC *fix_reax;

  double *chi,*eta,*gamma;

  int qeqflag;

  int setup_flag;

  int firstwarn;

  void allocate();

  void write_reax_atoms();