def main(directori, fitxer=None): if fitxer is None: fitxer = nom_fitxer_defectes fitxer = os.path.join(directori, fitxer) g = Gnuplot.Gnuplot() for frame_info, (walls, dislos, cincs, sets) in defectFrames(fitxer): title = "frame: %s t = %s strain: %s" % tuple(frame_info.split()) g.reset() g("set title '%s'" % title) if len(walls) > 0: g.replot(walls) if len(dislos) > 0: g.replot(dislos) g.replot(cincs) g.replot(sets)
def main(directori, fitxer=None): pyx.unit.set(defaultunit="cm", uscale=0.4) if fitxer is None: fitxer = nom_fitxer_defectes fitxer = os.path.join(directori, fitxer) foto = 0 yet = False for frame_info, (walls, dislos, cincs, sets) in defectFrames(fitxer): if not yet: if len(dislos) > 0: yet = True if foto % 5 == 0 and yet: sys.stderr.write("processant dislo%04d.eps\n" % foto) canvas = pyx.canvas.canvas() #for dislo in dislos: #gdislo = Dislo(dislo) #gdislo.color = pyx.color.rgb.blue #gdislo.showIn(canvas) for part in walls: gpart = Dislo(part) gpart.color = pyx.color.rgb.black gpart.showIn(canvas) for part in cincs: gpart = Dislo(part) gpart.color = pyx.color.rgb.red gpart.showIn(canvas) for part in sets: gpart = Dislo(part) gpart.color = pyx.color.rgb.green gpart.showIn(canvas) canvas.writeEPSfile("dislo%04d" % foto) foto = foto + 1 sys.exit(0)
return firstNucl def mainPersistent(directori, fitxer=None): if fitxer is None: fitxer = nom_fitxer_defectes fitxer = os.path.join(directori, fitxer) fitxerPersistent = os.path.join(directori, "./../firstNucleationsPersist") fitxerInfo = os.path.join(directori, "./../firstNucleations.Info") try: firstNucl = recupera(fitxerPersistent) except Exception, e: raise AbortDirectoryApply("fitxer: %s\nexception: %s" % (fitxer, e)) firstTime = True # False for frame_info, (walls, dislos, cincs, sets) in defectFrames(fitxer): if len(firstNucl) == 0: for part in walls: firstNucl.append(part) if len(dislos) > 0: if firstTime: for dislo in dislos: firstNucl.append(dislo) break else: firstTime = True open(fitxerInfo,"a").write("%s %s" % (directori, frame_info)) import pickle try: sys.stdout.write("guardo %s %d\n" % (directori, len(firstNucl)))