gmmbvl_sqdist - computes pairwise squared Euclidean distances between points
0001 function d = gmmbvl_sqdist(a,b) 0002 % gmmbvl_sqdist - computes pairwise squared Euclidean distances between points 0003 0004 % original version by Roland Bunschoten, 1999 0005 % 0006 % $Name: $ 0007 0008 if size(a,1)==1 0009 d = repmat(a',1,length(b)) - repmat(b,length(a),1); 0010 d = d.^2; 0011 else 0012 aa = sum(a.*a); bb = sum(b.*b); ab = a'*b; 0013 d = abs(repmat(aa',[1 size(bb,2)]) + repmat(bb,[size(aa,2) 1]) - 2*ab); 0014 end