Hello everyone,
I was using for a project the bond package (LIGGGHTS_FLEXIBLE_FIBERS) and I was trying to write out the number of bonds in a specific region of the domain. Does anybody know how one could do that?
With "thermo_style numbonds" one can print the number of bounds in the whole domain to the screen. However, when i try to assign numbonds keyword to some variable it does not work. It gives a completely different float number. Eventhough it works fine with other thermo_style keywords like "step" for example.
And when i try with:
compute "compute ID" all reduce/region "Region Name" numbonds
It comes a compute error. I'd appreciate any hint / help.