Exemplo n.º 1
0
for pw in [p[223]]:    #[0.58888888888888891, 0.72222222222222221]
    iw = [0.10, 0.4555555555555555555]
    runs.append(pw+iw) 
    




runner = []
prefix = "cont"
for px,py,ix,iy in runs:    
    rname = "%s_px%lfpy%lfix%lfiy%lfe" % (prefix,px,py,ix,iy)    
    runner.append("2dres %s.cfg > /dev/null  \n" % (rname))
    
    gencfg(rname, rname+".1.amdba")
    genpoly(rname+".poly", ix,iy,px,py,0.01,pointsB, segsB, 5)
    system("triangle -pqa0.0002 %s.poly" % (rname))
    system("triangle2amdba.pl %s.1" % (rname))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.poly" % (prefix,px,py,ix,iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.poly" % (prefix,px,py,ix,iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.ele" % (prefix,px,py,ix,iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.node" % (prefix,px,py,ix,iy)) 

    # call triangle
    # call triangle2amdba

#plot [0:1][0:1]"wloc 
#
nruns = 50
# for a,b,c,d in runs:
#     print c,d
Exemplo n.º 2
0
#print p[223]


for pw in [p[223]]:    #[0.58888888888888891, 0.72222222222222221]
    iw = [0.10, 0.4555555555555555555]
    runs.append(pw+iw) 
    




runner = []
prefix = "cont2"
for px,py,ix,iy in runs:    
    rname = "%s_px%lfpy%lfix%lfiy%lfe" % (prefix,px,py,ix,iy)    
    runner.append("2dres %s.cfg > /dev/null  \n" % (rname))
    
    gencfg(rname, rname+".1.amdba")
    genpoly(rname+".poly", ix,iy,px,py,0.01)
    system("triangle -pqa0.0002 %s.poly" % (rname))
    system("triangle2amdba.pl %s.1" % (rname))
    #remove("%s_px%lfpy%lfix%lfiy%lfe.1.poly" % (prefix,px,py,ix,iy))
    #remove("%s_px%lfpy%lfix%lfiy%lfe.poly" % (prefix,px,py,ix,iy))
    #remove("%s_px%lfpy%lfix%lfiy%lfe.1.ele" % (prefix,px,py,ix,iy))
    #remove("%s_px%lfpy%lfix%lfiy%lfe.1.node" % (prefix,px,py,ix,iy)) 

    # call triangle
    # call triangle2amdba

Exemplo n.º 3
0
        #print x,y

#print p[223]

for pw in [p[223]]:  #[0.58888888888888891, 0.72222222222222221]
    iw = [0.10, 0.4555555555555555555]
    runs.append(pw + iw)

runner = []
prefix = "cont"
for px, py, ix, iy in runs:
    rname = "%s_px%lfpy%lfix%lfiy%lfe" % (prefix, px, py, ix, iy)
    runner.append("2dres %s.cfg > /dev/null  \n" % (rname))

    gencfg(rname, rname + ".1.amdba")
    genpoly(rname + ".poly", ix, iy, px, py, 0.01, pointsB, segsB, 5)
    system("triangle -pqa0.0002 %s.poly" % (rname))
    system("triangle2amdba.pl %s.1" % (rname))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.poly" % (prefix, px, py, ix, iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.poly" % (prefix, px, py, ix, iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.ele" % (prefix, px, py, ix, iy))
    remove("%s_px%lfpy%lfix%lfiy%lfe.1.node" % (prefix, px, py, ix, iy))

    # call triangle
    # call triangle2amdba

#plot [0:1][0:1]"wloc
#
nruns = 50
# for a,b,c,d in runs:
#     print c,d