wolffd@0: function [dis_s12,dis_s13] = simple_dist(index_list,vector_list,index_map)
wolffd@0: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
wolffd@0: % compute the distance vector from 2 vector                                     %
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: 
wolffd@0: if all(all(checker(:,1:3)))
wolffd@0:   dis_s12 = abs(vector_list(index_list(:,1),:)-vector_list(index_list(:,2),:));
wolffd@0:   dis_s13 = abs(vector_list(index_list(:,1),:)-vector_list(index_list(:,3),:));
wolffd@0: else
wolffd@0:   fprintf('Some indices are not mapped\n');
wolffd@0: end
wolffd@0: 
wolffd@0: end
wolffd@0: % simple_dist(sontest{1},sonfeatures,sonindex)