wolffd@0: function [dis_s12 dis_s13] = conv_euclidean_dist(index_list,vector_list,index_map,w,norm)
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wolffd@0: % compute the distance vector from 2 vector using convolution euclidean         %
wolffd@0: % w = convolution window                                                        %
wolffd@0: % sontran2013                                                                   %
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wolffd@0: s_num = size(index_list,1);
wolffd@0: s_len = size(vector_list,2);
wolffd@0: 
wolffd@0: dis_s12 = zeros(s_num,s_len);
wolffd@0: dis_s13 = zeros(s_num,s_len);
wolffd@0: 
wolffd@0: % Convert list of assigned indices to ordered indices
wolffd@0: [checker index_list] = ismember(index_list,index_map);
wolffd@0: cv_vec = ones(1,w);
wolffd@0: if all(all(checker(:,1:3)))
wolffd@0:   dis_s12 = sqrt(conv2([1],cv_vec,(vector_list(index_list(:,1),:) - vector_list(index_list(:,2),:)).^2));
wolffd@0:   dis_s13 = sqrt(conv2([1],cv_vec,(vector_list(index_list(:,1),:) - vector_list(index_list(:,3),:)).^2));
wolffd@0:   % dis_s12 = dis_s12(:,w:end-w);
wolffd@0:   % dis_s13 = dis_s13(:,w:end-w);
wolffd@0:   if norm      
wolffd@0:       dis_s12 = dis_s12./(conv2([1],cv_vec,ones(s_num,s_len)));
wolffd@0:       dis_s13 = dis_s13./(conv2([1],cv_vec,ones(s_num,s_len)));
wolffd@0:   end
wolffd@0: else
wolffd@0:   fprintf('Some indices are not mapped\n');
wolffd@0: end
wolffd@0: end