我已经在molecules_group1中给出了一组分子,我想通过我的函数将这些分子移动到其他新的空molecules_group2上。它只在空的molecule_group2中添加了一个分子。每当我上传包含一个或多个分子的任何分子组时,它总是打印“排序的前五个分子存在于新组”,而不添加分子。此外,如果由于耐受性标准而没有添加分子,那么我如何使函数再次运行,直到成功添加分子。
def initial_configurations(newemptygroup,backgroundmols,tolerance,boxspace):
p = PointRef(backgroundmols.molecule().evaluate().centerOfGeometry())
c = CloseMols(p,newemptygroup,5)
g_random = c.closeMolecules()
if len(g_random) == 0:
newemptygroup.add(backgroundmols)
print("first molecules to be added")
else:
g = sorted(g_random.items(), key=operator.itemgetter(1)) # sort the dict by value which give tuple
print("sorted first five molecules present in newgroup =%s" %g)
t = [x[0] for x in g]
selected = newemptygroup[t[0]].molecule().evaluate().centerOfGeometry()
for i in range(1,len(t)):
closemols = newemptygroup[t[i]].molecule().evaluate().centerOfGeometry()
if (selected-closemols).length() > tolerance:
newemptygroup.add(backgroundmols)
else:
print("not added")
return newemptygroup发布于 2015-01-16 10:37:57
你可以在newemptygroup.add(backgroundmols)后面插入一条break语句
...
if (selected-closemols).length() > tolerance:
newemptygroup.add(backgroundmols)
break
else:
print("not added")
...这假设您会在for循环中的某个时候找到正确的背景工具。
https://stackoverflow.com/questions/27976353
复制相似问题