double r,r2inv,r6inv,forcecoul,forcelj,cforce,negforce;

  double factor_coul,factor_lj;

  double grij,expm2,prefactor,t,erfc;

  double philj;

  int iH1,iH2,jH1,jH2;

  double xiM[3],xjM[3];

  double *x1,*x2;

	  if (eflag) {

	    evdwl = r6inv*(lj3[itype][jtype]*r6inv-lj4[itype][jtype]) -

	      offset[itype][jtype];

	    evdwl *= factor_lj*philj;

	    evdwl *= factor_lj;

	  }

	} else evdwl = 0.0;