#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"
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):
# -*- 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"]:
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
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)