h1.setMinMaxProteinSize()
print('Surf', rSurf1.getMinMaxProteinSize())
print('Cyto', rCyto1.getMinMaxProteinSize())
print('o1', o1.getMinMaxProteinSize())
print('smallest', h1.smallestProteinSize)
print('largest', h1.largestProteinSize)
print('Bounding box', h1.boundingBox)
pad = 100.
afviewer.SetHistoVol(h1,pad,display=False)
afviewer.displayPreFill()
(x,y,z), maxi = h1.boundingBox
bb = [[0, y, z], maxi]
h1.buildGrid(boundingBox=bb)#, gridFileOut='1vesicle_new.grid' )
afviewer.displayFillBox(bb)
print('gridSpacing', h1.gridSpacing)
h1.printFillInfo()
#raw_input('press enter to start')
t1 = time()
h1.fill3(seedNum=0)
#
print('time to fill', time()-t1)
h1.printFillInfo()
#
afviewer.displayFill()
#
vertices = geomS.getVertices() vnormals = geomS.getVNormals() o1 = Organelle(vertices, faces, vnormals) h1.addOrganelle(o1) """ # set recipes h1.setExteriorRecipe(rCyto) #o1.setSurfaceRecipe(rSurf1) #o1.setInnerRecipe(rMatrix1) #again how interactively put ingrediants.... #maybe pred_displayFill can load everyavailable ingrediant....or at least create #an empty for everyone of them """print 'Bounding box', h1.boundingBox # add padding bb = h1.boundingBox pad = 100. x,y,z = bb[0] bb[0] = [x-pad, y-pad, z-pad] x,y,z = bb[1] bb[1] = [x+pad, y+pad, z+pad] print 'Bounding box with padding', h1.boundingBox (x,y,z), maxi = h1.boundingBox bb = [[0, y, z], maxi] #bb = [[0, 0, 0], maxi] h1.buildGrid(boundingBox=bb)#, gridFileOut='1vesicle_new.grid' ) #h1.buildGrid(gridFileOut='1vesicle_new.grid' )
vnormals = geomS.getVNormals() o1 = Organelle(vertices, faces, vnormals) h1.addOrganelle(o1) """ # set recipes h1.setExteriorRecipe(rCyto) #o1.setSurfaceRecipe(rSurf1) #o1.setInnerRecipe(rMatrix1) #again how interactively put ingrediants.... #maybe pred_displayFill can load everyavailable ingrediant....or at least create #an empty for everyone of them """print 'Bounding box', h1.boundingBox # add padding bb = h1.boundingBox pad = 100. x,y,z = bb[0] bb[0] = [x-pad, y-pad, z-pad] x,y,z = bb[1] bb[1] = [x+pad, y+pad, z+pad] print 'Bounding box with padding', h1.boundingBox (x,y,z), maxi = h1.boundingBox bb = [[0, y, z], maxi] #bb = [[0, 0, 0], maxi] h1.buildGrid(boundingBox=bb)#, gridFileOut='1vesicle_new.grid' ) #h1.buildGrid(gridFileOut='1vesicle_new.grid' )
print 'smallest', h1.smallestProteinSize
print 'largest', h1.largestProteinSize
print 'Bounding box', h1.boundingBox
# add padding
bb = h1.boundingBox
pad = 200.
x, y, z = bb[0]
bb[0] = [x - pad, y - pad, z - pad]
x, y, z = bb[1]
bb[1] = [x + pad, y + pad, z + pad]
print 'Bounding box with padding', h1.boundingBox
(x, y, z), maxi = h1.boundingBox
bb = [[0, y, z], maxi]
#bb = [[0, 0, 0], maxi]
h1.buildGrid(boundingBox=bb) #, gridFileOut='1vesicle_new.grid' )
#h1.buildGrid(gridFileOut='1vesicle_new.grid' )
execfile(plgDir + '/extension/testAF/c_displayPreFill.py')
print 'gridSpacing', h1.gridSpacing
#h1.printFillInfo()
#raw_input('press enter to start')
t1 = time()
h1.fill3(seedNum=0)
print 'time to fill', time() - t1
#h1.printFillInfo()
t1 = time()
