def onShrink(self, event): nstart,nstop=mys.get_range() # have xs,xs_star set in the end for nc in range(nstop-nstart+1): nsnap = nc + nstart if(not mys.snap_exists(nsnap)): continue print(mys.get_nhalo(nsnap)) for i in range(3):#range(mys.get_nhalo(nsnap)): if(not mys.exists_snap(nsnap)): print("snapshot "+str(nsnap)+" missing") continue cmd = "shrink_sphere.py "+str(nsnap)+" "+str(i+1) my.run(cmd) if(not mys.is_dmonly()): cmd = "shrink_sphere_stars.py "+str(nsnap)+" "+str(i+1) my.run(cmd) # TODO: if empty xs (shrinking sphere algorithm did not succeed: use xc) my.done()
def onShrink(self, event): nstart, nstop = mys.get_range() # have xs,xs_star set in the end for nc in range(nstop - nstart + 1): nsnap = nc + nstart if (not mys.snap_exists(nsnap)): continue print(mys.get_nhalo(nsnap)) for i in range(3): #range(mys.get_nhalo(nsnap)): if (not mys.exists_snap(nsnap)): print("snapshot " + str(nsnap) + " missing") continue cmd = "shrink_sphere.py " + str(nsnap) + " " + str(i + 1) my.run(cmd) if (not mys.is_dmonly()): cmd = "shrink_sphere_stars.py " + str(nsnap) + " " + str( i + 1) my.run(cmd) # TODO: if empty xs (shrinking sphere algorithm did not succeed: use xc) my.done()
from numpy import array, inner, sqrt import lib.mysql as mys import lib.initialize as my i = len(sys.argv) if (i != 3): print("shrink_sphere.py snap hid") exit(0) snap = int(sys.argv[1]) hid = int(sys.argv[2]) eps = 1e-5 frad = 0.90 # shrinked sphere has radius frad*maxr, i.e. frad = 0.90 means 10% smaller if (not mys.exists_snap(snap)): print("snapshot " + str(snap) + " missing") exit(0) xc, yc, zc, mvir, rvir = mys.getxyzmr(snap, 1) halo = my.open_file(mys.d(snap) + "dm/dm_" + str(hid) + ".dat", "r") x = [] y = [] z = [] m = [] for line in halo: val = line.split() m.append(float(val[0])) x.append(float(val[1])) y.append(float(val[2])) z.append(float(val[3]))
import lib.mysql as mys import lib.initialize as my i=len(sys.argv) if(i!=3): print("shrink_sphere.py snap hid") exit(0) snap=int(sys.argv[1]) hid=int(sys.argv[2]) eps=1e-5 frad = 0.9 # shrinked sphere has radius frad*maxr, i.e. frad = 0.90 means 10% smaller #print("missing file") if(not mys.exists_snap(snap)): print("snapshot "+str(snap)+" missing") exit(0) xc,yc,zc,mvir,rvir=mys.getxyzmr(snap,1) halo = my.open_file(mys.d(snap)+"stars/stars_"+str(hid)+".dat","r") x=[];y=[];z=[];m=[] for line in halo: val=line.split() m.append(float(val[0])) x.append(float(val[1])) y.append(float(val[2])) z.append(float(val[3])) x = array(x); y = array(y); z = array(z); m = array(m) #print(x)