od_rkhs.m

% computes optimality deficiency for each xk in the rkhs
% od = Z(x) - z(x) + log|x|
%   log2|x| / 8 to put cardinality into bytes
% maybe bytes is wrong? xk is in rkhs units, normalized and from ncd
function dx = od_rkhs(xk)

cx = repmat(mean(xk),size(xk,1),1);
dx = vecnorm(cx - xk, 2, 2) + log2(size(xk,1))/8;
4;