wolffd@0: function out = conv_subspace_delta(a, b, w, norm)
wolffd@0: % out = conv_subspace_delta(a, b, w, norm)
wolffd@0: % compute the distance vector from 2 vectors a,b using convolution euclidean         %
wolffd@0: % w = convolution window 
wolffd@0: %
wolffd@0: % DW adapted from sontran2013
wolffd@0: 
wolffd@0: % window for convolution
wolffd@0: 
wolffd@0: s_len = numel(a);
wolffd@0: if w == 0
wolffd@0:     w = s_len;
wolffd@0: end
wolffd@0: cv_vec = ones(w,1);
wolffd@0: 
wolffd@0: out = sqrt(conv(cv_vec,(a - b).^2));
wolffd@0: % out = out(:,w:end-w);
wolffd@0: if norm 
wolffd@0:   out = out./(conv(cv_vec,ones(s_len,1)));
wolffd@0: end
wolffd@0: 
wolffd@0: % how to get to 1 dimension?
wolffd@0: % out = sum(out);
wolffd@0: 
wolffd@0: