Page MenuHomec4science

atom_vec_gle.h
No OneTemporary

File Metadata

Created
Fri, Sep 13, 00:17

atom_vec_gle.h

/* ----------------------------------------------------------------------
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.
------------------------------------------------------------------------- */
/* ----------------------------------------------------------------------
* A class for an atom kind which carries around a variable number of
* additional degrees of freedom, to integrate a Generalized Langevin
* Equation dynamics.
*
* TODO: Probably this is better implemented as a hybrid type.
* ---------------------------------------------------------------------- */
#ifdef ATOM_CLASS
AtomStyle(gle,AtomVecGLE)
#else
#ifndef LMP_ATOM_VEC_GLE_H
#define LMP_ATOM_VEC_GLE_H
#include "atom_vec.h"
namespace LAMMPS_NS {
class AtomVecGLE : public AtomVec {
public:
AtomVecGLE(class LAMMPS *); //, int, char **); AtomVec classes used to accept arguments
virtual ~AtomVecGLE() {}
void grow(int);
void grow_reset();
void copy(int, int, int);
virtual int pack_comm(int, int *, double *, int, int *);
virtual int pack_comm_vel(int, int *, double *, int, int *);
virtual void unpack_comm(int, int, double *);
virtual void unpack_comm_vel(int, int, double *);
int pack_reverse(int, int, double *);
void unpack_reverse(int, int *, double *);
virtual int pack_border(int, int *, double *, int, int *);
virtual int pack_border_vel(int, int *, double *, int, int *);
int pack_border_hybrid(int, int *, double *);
virtual void unpack_border(int, int, double *);
virtual void unpack_border_vel(int, int, double *);
int unpack_border_hybrid(int, int, double *);
virtual int pack_exchange(int, double *);
virtual int unpack_exchange(double *);
int size_restart();
int pack_restart(int, double *);
int unpack_restart(double *);
void create_atom(int, double *);
void data_atom(double *, imageint, char **);
int data_atom_hybrid(int, char **);
void pack_data(double **);
int pack_data_hybrid(int, double *);
void write_data(FILE *, int, double **);
int write_data_hybrid(FILE *, double *);
bigint memory_usage();
protected:
int ns;
tagint *tag;
int *type,*mask;
imageint *image;
// s[i][j] i is the atom index, j contains 3*(ns+1) items
// corresponding to x,y,z coordinates of different
// additional DOFs, i.e. s[0] will be [ x0 y0 z0 x1 y1 z1 ... ]
double **x, **v, **f, **s;
double *q;
tagint *molecule;
int **nspecial;
tagint **special;
int *num_bond;
int **bond_type;
tagint **bond_atom;
int *num_angle;
int **angle_type;
tagint **angle_atom1,**angle_atom2,**angle_atom3;
int *num_dihedral;
int **dihedral_type;
tagint **dihedral_atom1,**dihedral_atom2,**dihedral_atom3,**dihedral_atom4;
int *num_improper;
int **improper_type;
tagint **improper_atom1,**improper_atom2,**improper_atom3,**improper_atom4;
};
}
#endif
#endif

Event Timeline