ICODE-ADC/lammps
4
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

add unique identifier to the history fix so the gran/*/history styles can be used multiple times with pair style hybrid

Browse Source
This commit is contained in:
Axel Kohlmeyer
2021-08-06 15:34:43 -04:00
parent e9e6cdca1d
commit 03f9ef7de9
1 changed files with 7 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
src/GRANULAR/pair_gran_hooke_history.cpp
View File

@ -68,7 +68,7 @@ PairGranHookeHistory::PairGranHookeHistory(LAMMPS *lmp) : Pair(lmp)
// this is so final order of Modify:fix will conform to input script
fix_history = nullptr;
modify->add_fix("NEIGH_HISTORY_HH_DUMMY all DUMMY");
modify->add_fix("NEIGH_HISTORY_HH_DUMMY"+std::to_string(instance_me)+" all DUMMY");
fix_dummy = (FixDummy *) modify->fix[modify->nfix-1];
}
@ -80,8 +80,8 @@ PairGranHookeHistory::~PairGranHookeHistory()
delete [] svector;
if (!fix_history) modify->delete_fix("NEIGH_HISTORY_HH_DUMMY");
else modify->delete_fix("NEIGH_HISTORY_HH");
if (!fix_history) modify->delete_fix("NEIGH_HISTORY_HH_DUMMY"+std::to_string(instance_me));
else modify->delete_fix("NEIGH_HISTORY_HH"+std::to_string(instance_me));
if (allocated) {
memory->destroy(setflag);
@ -436,10 +436,9 @@ void PairGranHookeHistory::init_style()
// this is so its order in the fix list is preserved
if (history && (fix_history == nullptr)) {
auto cmd = fmt::format("NEIGH_HISTORY_HH all NEIGH_HISTORY {}",
size_history);
modify->replace_fix("NEIGH_HISTORY_HH_DUMMY",cmd,1);
int ifix = modify->find_fix("NEIGH_HISTORY_HH");
auto cmd = fmt::format("NEIGH_HISTORY_HH{} all NEIGH_HISTORY {}", instance_me, size_history);
modify->replace_fix("NEIGH_HISTORY_HH_DUMMY"+std::to_string(instance_me),cmd,1);
int ifix = modify->find_fix("NEIGH_HISTORY_HH"+std::to_string(instance_me));
fix_history = (FixNeighHistory *) modify->fix[ifix];
fix_history->pair = this;
}
@ -507,7 +506,7 @@ void PairGranHookeHistory::init_style()
// set fix which stores history info
if (history) {
int ifix = modify->find_fix("NEIGH_HISTORY_HH");
int ifix = modify->find_fix("NEIGH_HISTORY_HH"+std::to_string(instance_me));
if (ifix < 0) error->all(FLERR,"Could not find pair fix neigh history ID");
fix_history = (FixNeighHistory *) modify->fix[ifix];
}