function [Lex, llrpOut,llrsOut]=RSC2decoder(FrameL,wnum,wl,Lsys,Lpari,LexIn)

%input:Lsys,Lpari1,Lpari2,Lu
%output:Lsys,Lpari1,Lpari2,Le

% FrameL=4096;
% wnum=64;
% wl=64;
% 
% Lsys=zeros(1,FrameL);
% Lpari=zeros(1,FrameL);
% % Lpari2=zeros(1,FrameL);
% LexIn=zeros(1,FrameL);
% 
% fid = fopen('E:/turbo/mfile/float_window/decoder_source_Lsysns.txt','r');
% Lsys=fscanf(fid,'%d\n '); 
% fclose(fid);
% 
% fid = fopen('E:/turbo/mfile/float_window/decoder_Lpari1ns.txt','r');
% Lpari=fscanf(fid,'%d\n '); 
% fclose(fid);
% 
% % fid = fopen('E:/turbo/mfile/float_window/decoder_Lpari2ns.txt','r');
% % Lpari2=fscanf(fid,'%d\n '); 
% % fclose(fid);

Gs=zeros(1,FrameL);
Gp=zeros(1,FrameL);
G0=zeros(1,FrameL);
G1=zeros(1,FrameL);
Gadd=zeros(1,FrameL);
Gsub=zeros(1,FrameL);

Gadd_A=zeros(1,FrameL);
Gsub_A=zeros(1,FrameL);
Gadd_B1=zeros(1,FrameL);
Gsub_B1=zeros(1,FrameL);
Gadd_B2=zeros(1,FrameL);
Gsub_B2=zeros(1,FrameL);


for i=1:FrameL
    Gs(i)=Lsys(i)+LexIn(i);
    Gp(i)=Lpari(i);
    G0(i)=Gs(i)/2;
    G1(i)=Gp(i)/2;
    Gadd(i)=G0(i)+G1(i);
    Gsub(i)=G0(i)-G1(i);
end

addrrA=gen_Aaddrr(FrameL);
addrrB1=gen_B1addrr(FrameL);
addrrB2=gen_B2addrr(FrameL);

for i=1:FrameL
Gadd_A(i)=Gadd(addrrA(i));
Gadd_B1(i)=Gadd(addrrB1(i));
Gadd_B2(i)=Gadd(addrrB2(i));
Gsub_A(i)=Gsub(addrrA(i));
Gsub_B1(i)=Gsub(addrrB1(i));
Gsub_B2(i)=Gsub(addrrB2(i));
end

A_metric=gen_A_metric(wnum,wl,Gadd_A,Gsub_A);

%   fid = fopen('E:/turbo/mfile/float_window/A_metric.txt','w');
%   fprintf(fid,'%f\n ',A_metric); 
%   fclose(fid);

B1_metric=gen_B1_metric(wnum,wl,Gadd_B1,Gsub_B1);

% fid = fopen('E:/turbo/mfile/float_window/B1_metric.txt','w');
%  fprintf(fid,'%f\n ',B1_metric); 
%  fclose(fid);

B2_metric=gen_B2_metric(wnum,wl,Gadd_B2,Gsub_B2);

% fid = fopen('E:/turbo/mfile/float_window/B2_metric.txt','w');
%  fprintf(fid,'%f\n ',B2_metric); 
%  fclose(fid);

%combine B1 and B2
B_metric=zeros(4,FrameL);

for i=1:wnum
    for j=1:wl
        if(mod(i,2))
          B_metric(1,(i-1)*wl+j)=B1_metric(1,(i-1)*wl+129-j);
          B_metric(2,(i-1)*wl+j)=B1_metric(2,(i-1)*wl+129-j);
          B_metric(3,(i-1)*wl+j)=B1_metric(3,(i-1)*wl+129-j);
          B_metric(4,(i-1)*wl+j)=B1_metric(4,(i-1)*wl+129-j);
        else
          B_metric(1,(i-1)*wl+j)=B2_metric(1,(i-2)*wl+129-j);
          B_metric(2,(i-1)*wl+j)=B2_metric(2,(i-2)*wl+129-j);
          B_metric(3,(i-1)*wl+j)=B2_metric(3,(i-2)*wl+129-j);
          B_metric(4,(i-1)*wl+j)=B2_metric(4,(i-2)*wl+129-j);
        end
    end
end

for j=1:64
B_metric(1,(64-1)*wl+j)=B2_metric(1,(64-2)*wl+65-j);
B_metric(2,(64-1)*wl+j)=B2_metric(2,(64-2)*wl+65-j);
B_metric(3,(64-1)*wl+j)=B2_metric(3,(64-2)*wl+65-j);
B_metric(4,(64-1)*wl+j)=B2_metric(4,(64-2)*wl+65-j);
end

% fid = fopen('E:/turbo/mfile/float_window/B_metric.txt','w');
%  fprintf(fid,'%f\n ',B_metric); 
%  fclose(fid);

 
%comput Le
Lex=zeros(1,FrameL);

Lex=LeOut(A_metric,B_metric,Gp,FrameL);

%  fid = fopen('E:/turbo/mfile/float_window/le2_laststage1.txt','w');
%  fprintf(fid,'%f\n ',Lex); 
%  fclose(fid);

LepOut=zeros(1,FrameL);

LepOut=Lep(A_metric,B_metric,Gs,FrameL);

llrpOut=zeros(1,FrameL);

for i=1:FrameL
llrpOut(i)=LepOut(i)+Gp(i);
end

% fid = fopen('E:/turbo/mfile/float_window/llrp_laststage1.txt','w');
%  fprintf(fid,'%f\n ',llrpOut); 
%  fclose(fid);

llrsOut=zeros(1,FrameL);

for i=1:FrameL
llrsOut(i)=Lex(i)+Gs(i);
end

%  fid = fopen('E:/turbo/mfile/float_window/llrs_laststage1.txt','w');
%  fprintf(fid,'%f\n ',llrsOut); 
%  fclose(fid);