lass="comment">%}</span>
0126 
0127 <span class="comment">% initial clustering, DSS</span>
0128 <span class="keyword">if</span> isempty(init)
0129     [C,A,score]=<a href="nt_cluster1D.html" class="code" title="function [C,A,score]=nt_cluster1D_b(x);">nt_cluster1D</a>(xx); <span class="comment">% cluster all columns of cross products</span>
0130     [~,idx]=min(score); <span class="comment">% select column with best score (tightest clusters)</span>
0131     A=A(:,idx); 
0132         
0133     <span class="c