logger = hoomd.logging.Logger()
St = hoomd.md.force.Force()
thermodynamic_properties = hoomd.md.compute.ThermodynamicQuantities(filter=all_)
sim.operations.computes.append(thermodynamic_properties)
logger.add(thermodynamic_properties,quantities=['kinetic_temperature','pressure_tensor','virial_ind_tensor','potential_energy'])
logger.add(sim,quantities=['tps'])
logger.add(morse,quantities=['virials'])
gsd_writer2 = hoomd.write.GSD(filename='p0.003_gel_poiseuille_12kT-6kT_LX=300_smallTrigger.gsd',trigger=1000,filter=all_,dynamic=['property','momentum','attribute'],mode='wb')
sim.operations.writers.append(gsd_writer2)
gsd_writer2.log = logger
'''
Hello,
If you want to average over every single frame (which is what it seems like you are suggesting) I would suggest creating a custom writer (documentation) that would run every frame and then average after being called n times. I would caution that averaging over every frame is likely not be useful given the high degree of correlations between consecutive frames. As for the binning you could use local snapshots to bin the particles (documentation).
Best,
Brandon
--
You received this message because you are subscribed to the Google Groups "hoomd-users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to hoomd-users...@googlegroups.com.
To view this discussion on the web visit https://groups.google.com/d/msgid/hoomd-users/ae5919b8-a187-43c4-9ace-4206590caf12n%40googlegroups.com.