diff --git a/src/MOLECULE/angle_hybrid.cpp b/src/MOLECULE/angle_hybrid.cpp index 869b2bf41a..6cd641e73b 100644 --- a/src/MOLECULE/angle_hybrid.cpp +++ b/src/MOLECULE/angle_hybrid.cpp @@ -282,6 +282,16 @@ void AngleHybrid::coeff(int narg, char **arg) } } +/* ---------------------------------------------------------------------- + run angle style specific initialization +------------------------------------------------------------------------- */ + +void AngleHybrid::init_style() +{ + for (int m = 0; m < nstyles; m++) + if (styles[m]) styles[m]->init_style(); +} + /* ---------------------------------------------------------------------- return an equilbrium angle length ------------------------------------------------------------------------- */ diff --git a/src/MOLECULE/angle_hybrid.h b/src/MOLECULE/angle_hybrid.h index b8c858224f..8f12bcb3db 100644 --- a/src/MOLECULE/angle_hybrid.h +++ b/src/MOLECULE/angle_hybrid.h @@ -36,6 +36,7 @@ class AngleHybrid : public Angle { void compute(int, int); void settings(int, char **); void coeff(int, char **); + void init_style(); double equilibrium_angle(int); void write_restart(FILE *); void read_restart(FILE *); diff --git a/src/angle.cpp b/src/angle.cpp index 54a3a5d051..709ace1826 100644 --- a/src/angle.cpp +++ b/src/angle.cpp @@ -53,6 +53,8 @@ void Angle::init() if (!allocated) error->all(FLERR,"Angle coeffs are not set"); for (int i = 1; i <= atom->nangletypes; i++) if (setflag[i] == 0) error->all(FLERR,"All angle coeffs are not set"); + + init_style(); } /* ----------------------------------------------------------------------