cd /afs/cs.cmu.edu/project/lemur-1/lkontor/Archive/CURRENT/matlab_work/string_proc/combinatorics
addpcurdir
cd ../../fsa/
M = read_dfsa('corinna2.txt')
M = bool2fsa_pt('(@)-(ab)','ab')
M = bool2fsa_pt('((@)-(aba))^((@)-(bab))','ab')
M = bool2fsa_pt('((@)-(aa))+((@)-(bb))','ab')
M = bool2fsa_pt(' (((@)-(ab)) + ((@)-(aa))) + (((@)-(ba)) + ((@)-(bb))) ','ab')
M = bool2fsa_pt('(((@)-(aa))+((@)-(bb)))-((aba)+(bab))','ab')
M = bool2fsa_pt('(((@)-(aa))+(ab))^(((@)-(bb))+(ba))','ab')

M = bool2fsa_pt(' ( ((@)-(aa)) + ((ba)+((@)-( (bb)+(ab) )))) ^ ( ((@)-(bb)) + ((ba)+((@)-( (aa)+(ba) ))))','ab')

M = bool2fsa_pt('((@)-(a))^((@)-(b))','abcdef')

M = bool2fsa_pt(' ((a)+(b))-((a)^(b)) ','ab')

global ALF
ALF = M.ALF

%Mcut = cutallreps_dfst(M.ALF);

%X = simple_paths(M);
%bb = run_dfsa_cell(M,X');


[X,Y] = reglang_posneg_dense(0,5,M);

if 1
  H = gram_matrix(X,@ptkermx,Y);
  [nsv, alpha, b0] = svc_gram(H,Y);
  
  ii=find(alpha>0.001)
  aa=alpha(ii).*Y(ii)
  Xii=X(ii)
  print_xa(Xii,aa)
end
  
%[Xt,Yt] = reglang_posneg_sparse(Mcor);
[Xt,Yt] = reglang_posneg_dense(6,7,M);
y_hat = compute_yhat_c(Xii,Xt,aa,b0,@ptkermx);
%y_hat = compute_yhat_vpcell(X,Y,Xt,v,@ptker2);

totdif(y_hat,Yt)

if 0
  
  X = lensort(X);
  Y = b2y(run_dfsa_cell(M,X(:)));
  
  good = 0;
  alpha = zeros(size(Y));
  while ~good
    %v = votedperceptron_cell(X,Y,1,@ptker2,v);
    %y_hat = compute_yhat_vpcell(X,Y,X,v,@ptker2);
    alpha = votedperceptron_mycell(X,Y,1,@ptkermx,alpha);  
    ii = find(alpha);
    Xii = X(ii);
    Yii = Y(ii);
    aa = alpha(ii).*Yii;
    y_hat = compute_yhat_c(Xii,X,aa,0,@ptkermx);
    
    nerr = length(find(y_hat~=Y))
    good = ~nerr;
  end

  ii = find(alpha);
  Xii = X(ii);
  Yii = Y(ii);
  aa = alpha(ii).*Yii;
  [Xt,Yt] = reglang_posneg_dense(8,9,M);  y_hat = compute_yhat_c(Xii,Xt,aa,0,@ptker2);
  
  
end