Hi Daniel,

The ‘combine’ function doesn’t add the created Molecule to chimera.openModels.  You would have to do that yourself with chimera.openModels.add().  Without having added ‘combination' to openModels, openModels.remove(combination) is a no-op, and therefore the C++ memory being used by the Molecule is not reclaimed.  You either have to use openModels.add() or you have to explicitly reclaim the C++ memory with combination.destroy().

—Eric

Eric Pettersen

UCSF Computer Graphics Lab

On Oct 4, 2018, at 6:51 AM, Daniel VH <dviladrich95@gmail.com> wrote:

Hi, In my python script I'm using the combine() function to merge a protein and a molecule into one model and since I'm doing this multiple times and only use the combined result to save it with pdbWrite I want to free the memory of the object combination right after that, but after the instructions chimera.openModels.close(combination), del combination, chimera.update.checkForChanges(), gc.collect(), chimera.closeSession() the memory is not freed. This is the code that reproduces the problem, it uses the molecules in the zipfile I have attatched.

import chimera

from Combine import combine

import zipfile

import time 

import gc

import os 

#from memory_profiler import profile

for itemnum in range(3):

#@profile

def pathmaker(itemnum):

dirpath = os.path.dirname(os.path.realpath(__file__))

zippath=dirpath+'/'+'test1_000.zip'

torsionpath=zipfile.ZipFile(zippath).extract('Pathway_000_Pathway.txt')

proteinpath=zipfile.ZipFile(zippath).extract('Protein_000_Protein.mol2')

ligandpath=zipfile.ZipFile(zippath).extract('Ligand_000_Ligand.mol2')

opened = chimera.openModels.open(proteinpath)

protein = opened[0]

opened1 = chimera.openModels.open(ligandpath)

mol=opened1[0]

combination=combine([mol,protein],protein)

totalnum=str(itemnum).zfill(3)+'_'+str(i).zfill(3)

chimera.pdbWrite([combination], chimera.Xform(), 'testpdb%s.pdb' %totalnum)

print('writing testpdb%s.pdb' %totalnum)

chimera.openModels.close(combination)

del combination

chimera.update.checkForChanges() 

gc.collect()

chimera.closeSession()

for i in range(20):

pathmaker(itemnum)

#this stops the code and shows memory increase for each iteration

print('pause 6 seconds...')

time.sleep(6)

<testing.py><test1_000.zip>_______________________________________________
Chimera-dev mailing list
Chimera-dev@cgl.ucsf.edu
http://plato.cgl.ucsf.edu/mailman/listinfo/chimera-dev