filep.close() xc = numpy.mean(xlist) yc = numpy.mean(ylist) zc = numpy.mean(zlist) for i in range(len(xlist)): xlist[i] = xlist[i] - xc ylist[i] = ylist[i] - yc zlist[i] = zlist[i] - zc nanoparticles_init = [] botx, topx, boty, topy, botz, topz = \ nanoparticle.file_to_nanoparticle_list(nanopfile, nanoparticles_init) print >> sys.stderr, "Read done" nanoparticles = [] for selectedid in range(len(nanoparticles_init)): nanop = nanoparticles_init[selectedid] pcx, pcy, pcz = nanop.get_center() # seleziona solo particelle interne distx = min(math.fabs(pcx-botx), math.fabs(pcx-topx)) disty = min(math.fabs(pcy-boty), math.fabs(pcy-topy)) distz = min(math.fabs(pcz-botz), math.fabs(pcz-topz)) #print min(distx, disty, distz), 2.0 * common.R
import vtk # no mi interessano le intersezioni nanoparticle.POINTINSIDEDIM = 0 filename = "nanoparticle_final_config.txt" if (len(sys.argv)) == 2: filename = sys.argv[1] nanoparticles = [] print ('Start reading '+filename+' ...') botx, topx, boty, topy, botz, topz = \ nanoparticle.file_to_nanoparticle_list(filename, nanoparticles) print ('Done') if (botx >= topx) or (boty >= topy) or \ (boty >= topy): print ('Error Invalid BOX\n') exit() # nanoparticle minbox_x = 100000.0 maxbox_x =-100000.0 minbox_y = 100000.0 maxbox_y =-100000.0 minbox_z = 100000.0 maxbox_z =-100000.0
plt.savefig("curve_" + id + ".png", bbox_inches='tight') # too big in case of many traps #fullsetoftraps = [] if args.filmfile != "": # non mi inter)essano le intersezioni nanoparticle.POINTINSIDEDIM = 0 nanoparticles = [] print >> sys.stderr, "Reading nanoparticles file" botx, topx, boty, topy, botz, topz = \ nanoparticle.file_to_nanoparticle_list(args.filmfile, \ nanoparticles) if (botx >= topx) or (boty >= topy) or (boty >= topy): print "Error Invalid BOX" exit() sources = [] counter = 0 for nanop in nanoparticles: p1, p2, tetha = nanop.get_rotation_info() cx, cy, cz = nanop.get_center() counter = counter + 1