#ifdef PAIR_CLASS
PairStyle(sph/taitwater/morris,PairSPHTaitwaterMorris)
#else
#ifndef LMP_PAIR_TAITWATER_MORRIS_H
#define LMP_PAIR_TAITWATER_MORRIS_H
#include "pair.h"
namespace LAMMPS_NS {
class PairSPHTaitwaterMorris : public Pair {
public:
PairSPHTaitwaterMorris(class LAMMPS *);
virtual ~PairSPHTaitwaterMorris();
virtual void compute(int, int);
void settings(int, char **);
void coeff(int, char **);
virtual double init_one(int, int);
protected:
double *rho0, *soundspeed, *B;
double **cut,**viscosity;
int first;
void allocate();
};
}
#endif
#endif