Page Menu
Home
c4science
Search
Configure Global Search
Log In
Files
F90987384
pair_hbond_dreiding_lj.h
No One
Temporary
Actions
Download File
Edit File
Delete File
View Transforms
Subscribe
Mute Notifications
Award Token
Subscribers
None
File Metadata
Details
File Info
Storage
Attached
Created
Wed, Nov 6, 16:44
Size
1 KB
Mime Type
text/x-c
Expires
Fri, Nov 8, 16:44 (1 d, 22 h)
Engine
blob
Format
Raw Data
Handle
22173638
Attached To
rLAMMPS lammps
pair_hbond_dreiding_lj.h
View Options
/* ----------------------------------------------------------------------
LAMMPS - Large-scale Atomic/Molecular Massively Parallel Simulator
http://lammps.sandia.gov, Sandia National Laboratories
Steve Plimpton, sjplimp@sandia.gov
Copyright (2003) Sandia Corporation. Under the terms of Contract
DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government retains
certain rights in this software. This software is distributed under
the GNU General Public License.
See the README file in the top-level LAMMPS directory.
------------------------------------------------------------------------- */
#ifdef PAIR_CLASS
PairStyle(hbond/dreiding/lj,PairHbondDreidingLJ)
#else
#ifndef LMP_PAIR_HBOND_DREIDING_LJ_H
#define LMP_PAIR_HBOND_DREIDING_LJ_H
#include "pair.h"
namespace LAMMPS_NS {
class PairHbondDreidingLJ : public Pair {
public:
PairHbondDreidingLJ(class LAMMPS *);
virtual ~PairHbondDreidingLJ();
virtual void compute(int, int);
void settings(int, char **);
virtual void coeff(int, char **);
virtual void init_style();
double init_one(int, int);
virtual double single(int, int, int, int, double, double, double, double &);
protected:
double cut_inner_global,cut_outer_global,cut_angle_global;
int ap_global;
struct Param {
double epsilon,sigma;
double lj1,lj2,lj3,lj4;
double d0,alpha,r0;
double morse1;
double denom_vdw;
double cut_inner,cut_outer,cut_innersq,cut_outersq,cut_angle,offset;
int ap;
};
Param *params; // parameter set for an I-J-K interaction
int nparams; // number of parameters read
int maxparam;
int *donor; // 1 if this type is ever a donor, else 0
int *acceptor; // 1 if this type is ever an acceptor, else 0
int ***type2param; // mapping from D,A,H to params, -1 if no map
void allocate();
};
}
#endif
#endif
Event Timeline
Log In to Comment