function lsp = lse_stat_par(sig, krnls, w)
  nr_krnls = size(krnls,1);
  sig_mat = repmat(sig,nr_krnls,1);
  w_mat   = repmat(w  ,nr_krnls,1);
  kk = sum(  krnls.*conj(krnls).*w_mat,2);
  ks = sum(sig_mat.*conj(krnls).*w_mat,2);
  lsp = log(ks/kk*2);
end