Fix massflow/mesh can not use 'delete_atoms' with fix multisphere

Submitted by yidongxia on Wed, 02/01/2017 - 00:27

It appears to be a huge inconvenience if 'fix massflow/mesh' cannot use 'delete_atoms' with 'fix multisphere', as the error shows below:

#########################################################################################
ERROR: Fix massflow/mesh (id massflow): can not use 'delete_atoms' with fix multisphere/* (../fix_massflow_mesh.cpp:304)
#########################################################################################

You will only see this error message if 'fix multisphere' comes before 'fix massflow/mesh' in input script. Strangely, this warning has been commented out in the code; see src/fix_massflow_mesh.cpp, line 332-334:

#########################################################################################
/*
if(!fix_ms_ && static_cast(modify->find_fix_style("multisphere",0)))
error->fix_error(FLERR,this,"fix multisphere must come before fix massflow/mesh in input script");*/
#########################################################################################

Is there a critical difficulty to allow particle deletion in fix_massflow_mesh together with fix_multisphere?

ckloss's picture

ckloss | Mon, 02/27/2017 - 20:59

Hi yidongxia ,

if you find this to be a huge inconvenience, feel free to implement it flawlessly and contribute the code :-)

happy coding and best wishes,
Christoph