pan> C(1)&lt;C(2); A=1-A; <span class="keyword">end</span> <span class="comment">% ensure that first cluster has low amplitude</span>
0188     
0189     A=double(<a href="nt_smooth.html" class="code" title="function x=nt_smooth(x,T,nIterations,nodelayflag)">nt_smooth</a>(A,smooth, [],1)&gt;=1/smooth); <span class="comment">% extend ownership to include effect of smoothing</span>
0190 
0191     <span class="comment">% upsample the cluster ownership index so we can apply it to x</span>
0192     A=repmat(A',[d