f1.write("\n") json.dump("-----------------------------", f1) f1.write("\n") json.dump("N=" + str(N), f1) f1.write("\n") json.dump("p_decode=" + str(design_p), f1) f1.write("\n") json.dump("capacity*N:" + str(C), f1) f1.write("\n") json.dump("sim ran :" + str(runsim), f1) f1.write("\n") try: I_order = pcon.getGCHsim('ZK', N, design_p, N) except: (I_order, E) = pcon.getGChZCK(design_p, N, N) used_rate = [] FER = [] Ratelist = [i * float(C) / N for i in deratelist] for R in Ratelist: if R <= (C / N): K = int(R * N) I = I_order[:K] #print R block_error = pch.polarchannelsim(N, design_p, design_p, I, runsim, False) used_rate.append(R) FER.append(block_error)
#---------------------------------------------------------------design: C = pl.CapacityBSC(N, design_p) deratepercentage = 10 tolerable_error = -2 if designbyrate: #K=int((100-deratepercentage)*C/100) K = int(C) try: I = pcon.getGCHsim('ZK', N, design_p, K) RI = pcon.getGChsim('ZK', design_p, N, N) except: (I, E) = pcon.getGChZCK(design_p, N, K) (RI, E) = pcon.getGChZCK(design_p, N, N) else: try: I = pcon.getGCHsim('ZL', N, design_p, tolerable_error) except: (I, E) = pcon.getGChZCL(design_p, N, tolerable_error) #==========================================================Construction print "LLR Lambda Channel REPORT" print "---------------------------" print "N=" + str(N) print "design_p=" + str(design_p) PotGCh = int(ma.floor(pl.CapacityBSC(N, design_p)))
def getreliability_order(N): return pcon.getGChZCK(0.01, N, N)[0]
# # Created: 19/08/2017 #---------------------------------------- import numpy as np import math as ma import problib as pl import polarencdec as ec import polarconstruct as pcon from datetime import datetime import json N=8 K=4 design_p=0.1 I=pcon.getGChZCK(design_p,N,K)[0] #print I G=len(I) #number of good channels FD=np.zeros(N-G,dtype=int).tolist()#frozen data YN=np.random.randint(2,size=N) YN=np.array([0,1,1,1,1,1,1,1]) #N=8,K=4 print "YN:" print YN.tolist()
channel_p=design_p runsim=100 stamp=datetime.now().strftime("%y-%m-%d_%H-%M-%S") f1=open("./simresults/timing_"+str(design_p)+"_"+stamp+".txt",'w') f2=open("./simresults/timplot_"+str(design_p)+"_"+stamp+".txt",'w') OneMBlist=[] #-------------------------------------------------------construction for n in nlist: N=2**n #bits start = timer() for i in range(runsim): I=pcon.getGChZCK(design_p,N,N/2)[0] end = timer() G=len(I) #number of good channels D=np.zeros(N-G,dtype=int).tolist()#frozen data print "Timing REPORT for each string (in Secs)" print "---------------------------------------" print "N="+str(N) print "p_channel="+str(channel_p) print "sim ran :"+str(runsim) print "Rate:"+str(float(G)/N) json.dump( "Timing REPORT for each string (in Secs)",f1) ;f1.write("\n") json.dump( "---------------------------",f1) ;f1.write("\n") json.dump( "N="+str(N),f1) ;f1.write("\n");