本文整理汇总了Python中Sire.System.monitors方法的典型用法代码示例。如果您正苦于以下问题:Python System.monitors方法的具体用法?Python System.monitors怎么用?Python System.monitors使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类Sire.System
的用法示例。
在下文中一共展示了System.monitors方法的2个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: print
# 需要导入模块: from Sire import System [as 别名]
# 或者: from Sire.System import monitors [as 别名]
new_molpro.setMolproExe(qmff.molproExe())
new_molpro.setEnergyOrigin(qmff.energyOrigin())
system_qmff = system.forceFields().forceField(qmff.ID())
new_molpro.addToQM(system_qmff.molecules(qmff.groups().qm()))
new_molpro.addToMM(system_qmff.molecules(qmff.groups().mm()), {new_molpro.parameters().coulomb(): "charges"})
print("QM energy = %f, new molproff = %f" % (system.forceFields().forceField(qmff.ID()).energy(), new_molpro.energy()))
qmff.change(system.forceFields().molecules())
print("qmff == %f" % qmff.energy())
# get the RDFs
rdfmonitor = system.monitors().monitor(Symbol("RDF"))
rdf = rdfmonitor.getRDF(oxygen, oxygen)
print("OXYGEN-OXYGEN")
for point in rdf.normalise():
print("%f %f" % point)
print("\nOXYGEN-HYDROGEN")
rdf = rdfmonitor.getRDF(oxygen, hydrogen)
for point in rdf.normalise():
print("%f %f" % point)
示例2: System
# 需要导入模块: from Sire import System [as 别名]
# 或者: from Sire.System import monitors [as 别名]
system = System(groups, ffields, monitors)
mc = RigidBodyMC( UniformSampler(solute) )
mc.setTemperature( 25 * celsius )
mc.setMaximumTranslation( 0.3 * angstrom )
#for i in range(0,50):
# moves = system.run(mc, 10000)
#
# mc = moves.moves()[0].clone()
# print "%d accepted, %d rejected, ratio = %f %%" % \
# (mc.nAccepted(), mc.nRejected(), 100 * mc.acceptanceRatio())
#
# print system.monitors().monitor(ffields.total()).average()
mtsmc = MTSMC(mc, e_fast.function(), 500)
mtsmc.setEnergyComponent(e_slow.function())
for i in range(0,5000):
moves = system.run(mtsmc, 20)
mtsmc = moves.moves()[0].clone()
print("%d accepted, %d rejected, ratio = %f %%" % \
(mtsmc.nAccepted(), mtsmc.nRejected(), 100 * mtsmc.acceptanceRatio()))
print("AVGENERGY: %f" % system.monitors().monitor(ffields.total()).average())
print("ENERGY: %f" % system.forceFields().energy())
sys.stdout.flush()