cd ../string_proc/combinatorics/
addpcurdir
cd ../../fsa/


global ALF
M=bool2fsa_lt('((ab)+(ba))-(((ab)^(ba)))','ab')
ALF = M.ALF;

%lambda = 10;
%ntrain = 100;
%ntrain = 10;
%ntest = 20;

[X,Y] = reglang_posneg_dense(0,8,M);
%[X,Y] = reglang_posneg_sparse(M,lambda,ntrain);

H = gram_matrix(X,@ngram_kernel,Y);
[nsv, alpha, b0] = svc_gram(H,Y);

ii = find(alpha>1e-5);
Yii = Y(ii);
aa = alpha(ii).*Yii;
Xii = X(ii);


[Xt,Yt] = reglang_posneg_dense(9,9,M);
y_hat = compute_yhat_c(Xii,Xt,aa,b0,@ngram_kernel);

nerr = length(y_hat~=Yt)