#also take the last bin
restorefilename= '2dboost'+str(logMh)+'_'+str(nas)+'_'+str(17)+'.sav'
restorefile=open(restorefilename,'r')
data=pickle.load(restorefile)
S= data['S']
#Load the boosts from S
ii= 0
jj= nas-1
while ii < nms:
    boostS[ii,jj]= S[ii,jj]
    ii= ii+1
restorefile.close()


#Then calculate the substructure boost for this halo mass, without the Sommerfeld enhancement
subboost= boost(Mh,1,1,1)

#Then calculate the reduced Sommerfeld enhancement for each grid point
#filename
savefilename= '2dredboost'+str(logMh)+'_'+str(nas)+'.sav'

if os.path.exists(savefilename):
    savefile=open(savefilename,'r+')
    data=pickle.load(savefile)
    ii= data['ii']
    jj= data['jj']
    if data['nms'] != nms:
        print "Warning: suppplied nms does not equal savefile nms"
        sys.exit(-1)
    if data['nas'] != nas:
        print "Warning: suppplied nas does not equal savefile nas"
Example #2
0
     sys.stdout.write('\r'+str(ii+1)+'/'+str(VLdata.shape[0]))
     sys.stdout.flush()
     if VLdata[ii,1] > main[6] or VLdata[ii,1] == 0.0:#Use 800 kpc?
         continue
     nlums+= 1
     lumdata.append(VLdata[ii,7])
     lumdata.append(VLdata[ii,8])
     lumdata.append(VLdata[ii,9])
     rs= VLdata[ii,4]/2.16258
     if rs > VLdata[ii,6]:
         print "Warning: rs outside of tidal radius"
     rhos= 8.56*10**4/pow(rs,2.)*pow(VLdata[ii,3],2.)
     lumdata.append(4.63*10**-24*pow(rhos,2.)*pow(rs,3.)*7.*pi/6.)
     #calculate total mass within rs
     mass= 4.*pi*rhos*pow(rs,3.)*massfactor
     boo= boost(mass,mv,a,0,1,0,mo,bo)
     mo= mass
     bo= boo
     lumdata.append(boo)
     lumdata.append(boost(mass,mv,a,1))
     vdisp= 2.5*10**-2*pow(mass,1./3.)/(3*10**5)
     lumdata.append(avg_enhance(mv,1.,a,vdisp))
     lumdata.append(rs)
     lumdata.append(rhos)
     lumdata.append(mass)
 sys.stdout.write('\n')
 #Now put them in a nice matrix form
 print "Number of bound subhalos:"
 print nlums
 sublum= zeros((nlums,natt))
 for ii in range(nlums):
Example #3
0
# -*- coding: utf-8 -*-
"""
Created on Thu Aug 30 16:39:43 2018

@author: Administrator
"""

import boost
import tf_idf
import selectData

if __name__ == '__main__':
    data = selectData.clusData(varname="title", tablename="qa_original_data")
    text = "我想知道该如何释放UIM卡异常状态"
    test = boost.boost(data)
    result_boost = test.getoriginData(text)
    for each in result_boost["match_sentence1"]:
        print(each)
    print("==================================================================")
    for each in result_boost["match_snetence2"]:
        print(each)

#    data = selectData.clusData(varname="title" , tablename="qa_original_data")
#    test_test = "CRM宽带资源的地址在哪里"
    result_idf = tf_idf.ifidf(data=data, text=text)
    print("==================================================================")
    for each in result_idf:
        print(each)

    for each in result_idf:
        if each in result_boost["match_snetence2"]:
Example #4
0
    if data['nas'] != nas:
        print "Warning: suppplied nas does not equal savefile nas"
        sys.exit(-1)
    mv= data['mv']
    a= data['a']
    S= data['S']
    print 'restarting at ii= '+str(ii)+', jj= '+str(jj)+', at iteration '+str(ii*nas+jj+1)+'/'+str(nms*nas)
else:
    S= zeros((nms,nas))
    ii= 0
    jj= 0
    savefile=open(savefilename,'w')

while ii < nms:
    while jj < nas:
        boo= boost(Mh,mv[ii],a[jj])
        so= avg_enhance(mv[ii],m,a[jj],vdisp)
#        print mv[ii]
#        print a[jj]
#        print boo
#        print so
        S[ii,jj]= log10(so+boo)
#        print S[ii,jj]
        sys.stdout.write('\r'+str(ii*nas+jj+1)+'/'+str(nms*nas))
        sys.stdout.flush()
        jj= jj+1
        data={'mv':mv,'a':a,'S':S,'nms':nms,'nas':nas,'ii':ii,'jj':jj}
        savefile.seek(0,0)
        pickle.dump(data,savefile)
        savefile.flush()
    ii= ii+1
Example #5
0
    a= data['a']
    S= data['S']
    print 'restarting at ii= '+str(ii)+', jj= '+str(jj)+', at iteration '+str(ii+(jj-(bin-1)*(nas-1)/nbins)*nms+1)+'/'+str(nms*nas/nbins)
else:
    S= zeros((nms,nas))
    ii= 0
    if bin != 17:
        jj= (bin-1)*(nas-1)/nbins
    else:
        jj= nas-1
    savefile=open(savefilename,'w')

if bin != 17:
    while jj < bin*(nas-1)/nbins:
        while ii < nms:
            boo= boost(M,mv[ii],a[jj],0,gamma,einasto)
            so= avg_enhance(mv[ii],m,a[jj],vdisp)
            S[ii,jj]= log10((so+boo)*7.*pi/6.*pow(rhos,2.)*pow(rs,3.)/pow(D,2)*4.63*10**-24)
            sys.stdout.write('\r'+str(ii+(jj-(bin-1)*(nas-1)/nbins)*nms+1)+'/'+str(nms*(nas-1)/nbins))
            sys.stdout.flush()
            ii= ii+1
            data={'mv':mv,'a':a,'S':S,'nms':nms,'nas':nas,'ii':ii,'jj':jj}
            savefile.seek(0,0)
            pickle.dump(data,savefile)
            savefile.flush()
        jj= jj+1
        ii= 0
    sys.stdout.write('\n')
else:
    while ii < nms:
        boo= boost(M,mv[ii],a[jj],0,gamma,einasto)