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
#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
#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