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");
