Parallel config added.

Author: Alex

config.m

@@ -4,7 +4,7 @@
 %% Configure parameters
 %Polar-Code values
 capacity = 0.5; %I(W), Channel's W Capacity
-n_values = [4,5,6,7,8,9,10];     %value of N
+n_values = [4,5,6,7,8];     %value of N
 code_rate = 1/4;
 %EbNo
 EbNo_dB = 0:5; %AWGN -4:2 %Fading 0:2:10
@@ -28,6 +28,30 @@
 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
+parcore_nums = 4;
+if FLAG_Enable_parpool
+    % determine the number of physical cores
+    corenum = feature('numcores');
+    parcorenum = parcore_nums;
+    % parcorenum = 15; %  # of workers
+    p = gcp('nocreate');   % Not create new pool if it does not exist
+    if isempty(p)
+        p = parcluster('local');
+        p.NumWorkers = parcorenum;
+        parpool(p, p.NumWorkers);
+    end
+    if p.NumWorkers ~= parcorenum
+        delete(p);
+        p = parcluster('local');
+        p.NumWorkers = parcorenum;
+        parpool(p, p.NumWorkers);
+    end
+    parallel_frames = 100*parcorenum;
+else
+    parallel_frames = 10;
+end
 %% Variable Initializations
 bit_error_rate = zeros(length(n_values),length(snrdb_values));
 fer_error_rate = zeros(length(n_values),length(snrdb_values));

main.m

@@ -15,7 +15,10 @@
         codewords_tmp = 0;fer_errors = 0;bit_errors = 0;
         snr = snrdb_values(i_index);
         while (fer_errors<min_fer_errors || codewords_tmp<min_codewords)
-            codewords_tmp = codewords_tmp + 1;
+            bit_errors_parfor = zeros(1,parallel_frames);
+            fer_errors_parfor = zeros(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);
@@ -38,8 +41,12 @@
             %Calculate temporary bit/frame errors
             final_outputs = transform_outputs(outputs,frozen_bits,N);
             temp_bit_errors = size(find(final_outputs ~= inputs),2);
-            bit_errors = bit_errors + temp_bit_errors;
-            fer_errors = fer_errors + (temp_bit_errors>0);
+            bit_errors_parfor(frame) = temp_bit_errors;
+            fer_errors_parfor(frame) = (temp_bit_errors>0);
+            end
+            codewords_tmp = codewords_tmp + parallel_frames;
+            bit_errors = bit_errors + sum(bit_errors_parfor);
+            fer_errors = fer_errors + sum(fer_errors_parfor);
         end
         %Calculate snr bit/frame errors
         bit_error_rate(index_n,i_index) = bit_errors/(codewords_tmp*K);