function [C, Ceq] = strong_stat_constr(P)
% forces a process to be w-strongly stationary

global m n w

S = 0;
for k=1:n-w+1
  Pk = getmargX(P,m,n,[k:k+w-1]);
  %[b,del] = is_indep(Pk,m,n,w);
  [hh,Ht,H] = gethhn(Pk,m,w);
  S = S+H-1;
end
Ceq = 0;
C = S;

return