示例#1
0
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)
示例#2
0
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)
示例#3
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)))