#ifdef PAIR_CLASS
PairStyle(tri/lj,PairTriLJ)
#else
#ifndef LMP_PAIR_TRI_LJ_H
#define LMP_PAIR_TRI_LJ_H
#include "pair.h"
namespace LAMMPS_NS {
class PairTriLJ : public Pair {
public:
PairTriLJ(class LAMMPS *);
virtual ~PairTriLJ();
virtual void compute(int, int);
void settings(int, char **);
void coeff(int, char **);
virtual void init_style();
double init_one(int, int);
protected:
double cut_global;
double **cut;
double **epsilon,**sigma;
double **lj1,**lj2,**lj3,**lj4;
class AtomVecTri *avec;
struct Discrete {
double dx,dy,dz;
double sigma;
};
Discrete *discrete; int ndiscrete; int dmax; int *dnum; int *dfirst; int nmax;
void allocate();
void discretize(int, double, double *, double *, double *);
};
}
#endif
#endif