stamp = datetime.now().strftime("%y-%m-%d_%H-%M-%S")
f1 = open("./simresults/lambda_theta_LT" + stamp + ".txt", 'w')
#---------------------------------------------------------------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"
json.dump("RATE Vs FER REPORT derate pch", f1)
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)
stamp = datetime.now().strftime("%y-%m-%d_%H-%M-%S")
f1 = open("./simresults/lambda_avg" + stamp + ".txt", 'w')
#==========================================================Construction
print "LLR Lambda REPORT"
print "---------------------------"
print "N=" + str(N)
print "design_p=" + str(p)
PotGCh = int(ma.floor(pl.CapacityBSC(N, p)))
print "Capacity for " + str(N) + "channels:" + str(PotGCh)
print "tolerable error exponent:" + str(tolerable_error)
print "sim ran :" + str(runsim)
#-----------------------------------------ZC
#(I,E)=pcon.getGChZCL(p,N,tolerable_error)
#(I,E)=pcon.getGChZCK(p,N,K)
I = pcon.getGCHsim("MK_ALL", N, p, K)
print "Good Channels:"
print I
print "Number of good channels:"
print len(I)
R = float(len(I)) / N
print "R=" + str(R)
print "Frozen channels:"
B = list(set(range(N)) - set(I))
print len(B)
##########################################file
json.dump("LLR Lambda CHANNEL REPORT", f1)
f1.write("\n")
json.dump("---------------------------", f1)
f1.write("\n")
json.dump("N=" + str(N), f1)
