function [Hs,gs,Fs,bs,Js,Ls] = mySoftPadding(H,F,bb,J,L,S,rho,m)
% S = weight for quadratic cost of constraint violations
% rho = scalar weight for 1-norm of constraint violations
% m = number of inputs

% your code here

Nm = size(H,1);
N = Nm/m;
ss = size(S,1);
Nss = N*ss;

Zmn = zeros(Nm,Nss);

Hs = [H,Zmn;Zmn',kron(2*eye(N),S)];

gs = rho*(ones(Nss, 1));

Iss = eye(ss);
Zmmss = zeros(2*m,ss);
It = [-Iss;-Iss;Zmmss];

Id = kron(eye(N),It);   %2N(ss+m)xNss

Fs = [F,Id;zeros(Nss,Nm),-eye(Nss)]; %2N(ss+m)+Nss x Nm+Nss

bs = [bb;zeros(Nss,size(bb,2))]; %2N(m+ss)+ Nss x1

Zj = zeros(Nss,size(J,2));
Js = [J;Zj];  %2N(m+ss)+ Nss xn
Ls = [L;Zj];  %2N(m+ss)+ Nss xn

end