#ifdef PAIR_CLASS
PairStyle(edpd,PairEDPD)
#else
#ifndef LMP_PAIR_EDPD_H
#define LMP_PAIR_EDPD_H
#include "pair.h"
namespace LAMMPS_NS {
class PairEDPD : public Pair {
public:
PairEDPD(class LAMMPS *);
virtual ~PairEDPD();
virtual void compute(int, int);
virtual void settings(int, char **);
virtual void coeff(int, char **);
void init_style();
double init_one(int, int);
virtual void write_restart(FILE *);
virtual void read_restart(FILE *);
virtual void write_restart_settings(FILE *);
virtual void read_restart_settings(FILE *);
double single(int, int, int, int, double, double, double, double &);
protected:
double cut_global;
int seed;
double **cut,**cutT;
double **a0,**gamma;
double **power;
double **slope;
double **kappa;
double **powerT;
int power_flag, kappa_flag;
double ***sc,***kc;
class RanMars *random;
class RanMars *randomT;
void allocate();
};
}
#endif
#endif