1function [w1]=NCAF(K,w,order,n,mcAdapt)
2
3 if ~isempty(K)
4
5 Omega = sum(w.^2);
6 OmegaMat = ones(order,1) * Omega;
7 inds = OmegaMat > K;
8 if ~isempty(mcAdapt)
9 amcIndMat = ones(order,1) * mcAdapt(n,:);
10 inds = OmegaMat > K & amcIndMat;
11 end
12 w(inds) = w(inds) .* sqrt(K./OmegaMat(inds));
13 end
14w1=w;
15end