Extra complexity removed. Mex filed added.

Author: Alex

config.m

@@ -14,23 +14,30 @@
 Fading_Independent = 0; %if Fading_Channel = 1
 quasi = 0;              % quasi channel
 fading_channel = 'TU120';   %custom matlab channel
-
+%seed
 fix_seed = 1;           %1:fix data, 0:random data
-if fix_seed
-    rng(sum(100*clock));
-else
-    rand('seed',123456);
-end
 %Simulation values
-fast_run = 1;           %1:run of optimal-matlab code, 0:run of hardware-code (suboptimal f/g)
+fast_run = 1;           %1:run of optimal-matlab code + mex-files, 0:run of hardware-code (suboptimal f/g)
 min_fer_errors = 100;   %minimum frame errors to count
 min_codewords = 100;    %minimum codewords to count
 NbitsPerSymbol = 1;     %modulation parameter
 constDims = 1;          %modulation parameter
 snrdb_values =EbNo_dB+10*log10(double(code_rate*NbitsPerSymbol*2/constDims));
 %% Parfor configuration
-FLAG_Enable_parpool=1;    % 0: Disable, 1: Enable
+FLAG_Enable_parpool=0;    % 0: Disable, 1: Enable
 parcore_nums = 4;
+%% Initializations
+%add mex-files
+if(fast_run == 1)
+    addpath('./support/mex_files');
+end
+%initialize seed
+if fix_seed
+    rng(sum(100*clock));
+else
+    rand('seed',123456);
+end
+%initialize parcorenum
 if FLAG_Enable_parpool
     % determine the number of physical cores
     corenum = feature('numcores');
@@ -50,7 +57,7 @@
     end
     parallel_frames = 100*parcorenum;
 else
-    parallel_frames = 10;
+    parallel_frames = 100;
 end
 %% Variable Initializations
 bit_error_rate = zeros(length(n_values),length(snrdb_values));

main.m

@@ -7,7 +7,7 @@
     K = N* code_rate;                               %Code keyword length
     mat_file = [result_path 'polar_N' num2str(N) '_K' num2str(K) '_' timestamp '.mat'];
     %Polar-code initializations
-    [Fn,frozen_bits,partial_sum_adders, sc_functions, sc_2nd_indxs] = polar_initialization(N, K, capacity);
+    [Fn,frozen_bits, frozen_indxs, non_frozen_indxs, partial_sum_adders, sc_functions, sc_2nd_indxs] = polar_initialization(N, K, capacity);
     fprintf("Polar Code %d/%d running:\n",N,K);
     %% SNR LOOP
     for i_index = 1:length(snrdb_values)
@@ -17,11 +17,11 @@
         while (fer_errors<min_fer_errors || codewords_tmp<min_codewords)
             bit_errors_parfor = zeros(1,parallel_frames);
             fer_errors_parfor = zeros(1,parallel_frames);
-%             for frame = 1:parallel_frames
-            parfor frame = 1:parallel_frames
+            for frame = 1:parallel_frames
+%             parfor frame = 1:parallel_frames
             inputs = rand(1,K)>0.5;    %write random inputs
             %transform inputs
-            inputs_to_encode = transform_inputs(inputs,frozen_bits,N);
+            inputs_to_encode = transform_inputs(inputs,non_frozen_indxs,N);
             %encode
             encoded_inputs = encode(inputs_to_encode,Fn);      %Reversed Polar Encoding
             %modulate
@@ -38,9 +38,9 @@
                 %hardware-version (suboptimal-calculations of f/g)
                 outputs = decode(llr,frozen_bits,partial_sum_adders,sc_functions,sc_2nd_indxs); %or decode2(llr,frozen_bits); for the other algorithm
             end
+            final_outputs = outputs(non_frozen_indxs);%transform_outputs(outputs,frozen_bits,N);
             %Calculate temporary bit/frame errors
-            final_outputs = transform_outputs(outputs,frozen_bits,N);
-            temp_bit_errors = size(find(final_outputs ~= inputs),2);
+            temp_bit_errors = sum(final_outputs ~= inputs);
             bit_errors_parfor(frame) = temp_bit_errors;
             fer_errors_parfor(frame) = (temp_bit_errors>0);
             end

support/mex_files/l_f.mexw64

@@ -1,8 +1,10 @@
-function [ Fn,frozen_bits,partial_sum_adders, sc_functions, sc_2nd_indxs] = polar_initialization( N, K, capacity )
+function [ Fn,frozen_bits, frozen_indxs, non_frozen_indxs, partial_sum_adders, sc_functions, sc_2nd_indxs] = polar_initialization( N, K, capacity )
 %POLAR_INITIALIZATION Summary of this function goes here
 %   Detailed explanation goes here
 Fn = fkronecker(N);
 frozen_bits = initialize_frozen_bits(N,K,capacity); %0=frozen, 1=not_frozen
+frozen_indxs = find(frozen_bits == 0);
+non_frozen_indxs = find(frozen_bits == 1);
 partial_sum_adders = partial_sums_initialize(N);   %!!!!! NOT bit_reversed array -- %outputs(z,i,l) -- l stage ,bit Ui is added,z is the number of g adder
 [sc_functions,sc_2nd_indxs] = sc_array_initialize(N);
 end

support/mex_files/s_f.mexw64

@@ -10,7 +10,7 @@
     end
 else
     if (mod(floor((j-1)/power(2,l-2)),2) == 0)
-        output = xor(s_f(l-1,j,frozen_bits,estimated),s_f(l-1,j+power(2,l-2),frozen_bits,estimated));
+        output = double(xor(s_f(l-1,j,frozen_bits,estimated),s_f(l-1,j+power(2,l-2),frozen_bits,estimated)));
     else
         output = s_f(l-1,j,frozen_bits,estimated);
     end

support/polar_initialization.m

@@ -1,6 +1,7 @@
 function [sc_function,sc_2nd_indxs] = sc_array_initialize(bits)
 %decoding initialization
-outputs = zeros(bits,log2(bits)+1,3);%outputs(i-1,l-1,dimension) // dimension= 1 values/dimension=2 if 0 then f else g/dimension=3 2nd input
+%sc_function:if 0 then f else g
+%sc_2nd_indxs: index of 2nd input
 sc_function = zeros(bits,log2(bits)+1);
 sc_2nd_indxs = zeros(bits,log2(bits)+1);
 reverse_i = bitrevorder(1:1:bits);

support/s_f.m

@@ -1,7 +1,7 @@
-function outputs = transform_inputs(inputs,frozen_bits,N)
+function outputs = transform_inputs(inputs,non_frozen_indxs,N)
 %takes Xi inputs,the position of frozen bits, and makes the outputs array
 %inputs at not frozen positions, 0 at frozen positions
 outputs = zeros(1,N);
-outputs(frozen_bits == 1) = inputs;
+outputs(non_frozen_indxs) = inputs;
 % correct