%?CSSA1.m:??CONDENSING?SORT?AND?SELECT?ALGORITHM?（1-d）
%
%?B:?????input?data
%?L:?????number?of?levels?in?wavelet?transform
%?vol:???stopping?volume?of?mask?（default?is?full?volume）
%????????NOTE:?does?NOT?include?scaling?coeff?volume!
%
%?mask:??final?supernode?configurations?ordered?
%
%?Assumes?that?the?underlying?tree?is?the?1-d?dyadic?wavelet?tree
%
%?RGB?INI?September?1998

function?mask?=?cssa1（BLvol）;

%---------------------------------------------------------------------------%
%?INITIALIZATIONS

%?make?input?data?a?column?vector
B????=?B（:）;
N????=?length（B）;
cvol?=?0;
sno??=?2;
%?scaling?coefficient?with?highest?index?in?wavelet?array
parroot?=?N/2^L;

if?nargin???vol?=?（N-N/2^L）;
end
if?vol?>?（N-N/2^L）
??disp（‘ERROR:?volume?parameter?cannot?be?>?number?of?wavelet?coeffs‘）
??mask?=?[];
??return
end

%?set?up?each?node?as?a?supernode.??use?an?Nx5?matrix:
%????col?1:?utp??（uptree?pointer）
%????col?2:?rp???（supernode?root?pointer）
%????col?3:?snv??（supernode?value）
%????col?4:?num??（number?of?internal?nodes）
%????col?5:?mask?（mask/kernel?value）

sn（1:N?1）?=?（1:N）‘;??????%?point?uptree?top?pointers?to?self
sn（1:N?2）?=?-ones（N1）;??%?each?node?is?a?sn?root?=>?-1?flag
sn（1:N?3）?=?B;???????????%?each?snv?is?the?data?itself
sn（1:N?4）?=?ones（N1）;???%?num?=?1?for?each?sn???
sn（1:N?5）?=?zeros（N1）;??%?output?mask?is?zero?at?start
sn（1:N/2^L?5）?=?ones（N/2^L1）;??%?except?for?parents?of?roots?（“zero“?nodes）

%?set?up?SNV?queue??（include?scaling?coeff?values?just?so?that?Q?never?empties
snvQ（1:N?1）?=?（1:N）‘;????%?each?node?is?its?own?sn?topnode
snvQ（1:N?2）?=?B;?????????%?each?snv?is?the?wavelet?data?itself
snvQ（1:N/2^L2）?=?NaN*snvQ（1:N/2^L2）;


%---------------------------------------------------------------------------%
%?MAIN?LOOP?（follow?pseudo?code?on?p.?138?and?steps?on?p.?142）

%?notation:???MAX?refers?to?the?current?sn?we?take?from?snv?queue
%?????????????CHECK?refers?to?the?parent?sn?it?abutts

while?（cvol?
??%?STEP?1
??%?find?supernode?with?largest?SNV?（sn?MAX）
??%?delete?it?from?the?snvQ?as?well
??topmax?=?Qfind（snvQ）;
??snvQ?=?Qdel（snvQtopmax）;
??
??%?STEP?2
??%?find?the?root?node?of?sn?MAX?and?then?its?parent
??%?this?latter?node?is?a?leaf?of?the?CHECK?sn?above?ours
??if?sn（topmax2）?==?-1
????%?then?topmax?is?the?sn?root
????leafcheck?=?ceil（topmax/2）;
??else
????%?go?to?the?sn?root
????leafcheck?=?ceil（sn（topmax2）/2）;
??end
??
??
??%?check?whether?CHECK?sn?has?mask=10?（first?find?top?node?of?CHECK）
%??topcheck?=?Tfind（snleafcheck）;
??[topchecksn]?=?Tfind（snleafcheck）;???????%?need?to?pass?sn?back
?????????????????????????????????????????????%?if?we?do?path?comp?step
????
??if?sn（topcheck5）?>?0?
????%?we?are?in?the?money?-?CHECK?sn?has?mask?>?0.??so?take?sn?MAX
????%?STEP?3?（usual?case?-?not?the?first?sn?taken）:?
????sn?=?umask（sntopmaxsnoN）;
????cvol?=?cvol?+?sn（topmax4）;?
????sno?=?sno?+?1;
?
??else
????%

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????6148??2009-08-05?19:08??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\.DS_Store

?????文件???????3184??2008-07-15?12:35??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\cssa1.m

?????文件???????1960??2008-07-15?12:35??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\greedy1.m

?????文件???????2090??2009-08-05?10:40??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\treemp.m

?????文件???????1747??2009-08-04?16:48??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\treemp_example.m

?????文件???????2022??2009-08-04?14:58??RICE?UNIVERSITY代碼\modelcs_v1.1\1D?trees\treemp_fun.m

?????文件???????6148??2009-08-05?19:08??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\.DS_Store

?????文件???????1565??2008-08-13?19:10??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\cosamp_2D_nf.m

?????文件???????3938??2008-07-30?15:14??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\cssa2.m

?????文件???????2451??2008-07-15?12:35??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\greedy2.m

?????文件?????538749??2008-08-13?17:10??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\peppers.png

?????文件???????1974??2009-08-05?10:34??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\treemp_2D_example.m

?????文件???????2439??2009-08-05?10:42??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\treemp_greedy_2D.m

?????文件???????2611??2009-08-05?10:26??RICE?UNIVERSITY代碼\modelcs_v1.1\2D?trees\treemp_nf_2D.m

?????文件???????6148??2009-08-05?19:08??RICE?UNIVERSITY代碼\modelcs_v1.1\Block?sparsity\.DS_Store

?????文件???????2346??2009-08-05?11:06??RICE?UNIVERSITY代碼\modelcs_v1.1\Block?sparsity\jsmp.m

?????文件???????1301??2009-08-05?11:04??RICE?UNIVERSITY代碼\modelcs_v1.1\Block?sparsity\jsmp_example.m

?????文件???????2454??2009-08-05?11:06??RICE?UNIVERSITY代碼\modelcs_v1.1\Block?sparsity\jsmp_fun.m

?????文件???????6148??2009-08-05?19:09??RICE?UNIVERSITY代碼\modelcs_v1.1\Clustered?sparsity\.DS_Store

?????文件???????1873??2009-08-05?13:40??RICE?UNIVERSITY代碼\modelcs_v1.1\Clustered?sparsity\kc.m

?????文件????????632??2009-08-05?13:37??RICE?UNIVERSITY代碼\modelcs_v1.1\Clustered?sparsity\kc_example.m

?????文件???????6148??2009-08-05?19:09??RICE?UNIVERSITY代碼\modelcs_v1.1\Delta\.DS_Store

?????文件????????659??2009-08-05?14:28??RICE?UNIVERSITY代碼\modelcs_v1.1\Delta\bestdelta.m

?????文件???????2468??2009-08-05?14:44??RICE?UNIVERSITY代碼\modelcs_v1.1\Delta\deltarec.m

?????文件????????760??2009-08-05?15:00??RICE?UNIVERSITY代碼\modelcs_v1.1\Delta\delta_example.m

?????文件???????2678??2012-02-15?09:19??RICE?UNIVERSITY代碼\modelcs_v1.1\README.txt

?????文件???????6148??2009-08-05?19:09??RICE?UNIVERSITY代碼\modelcs_v1.1\Sparsity\.DS_Store

?????文件??????????0??2010-10-01?16:58??RICE?UNIVERSITY代碼\modelcs_v1.1\Sparsity\;.m

?????文件???????1743??2009-07-31?15:14??RICE?UNIVERSITY代碼\modelcs_v1.1\Sparsity\cosamp.m

?????文件???????1127??2009-08-04?11:20??RICE?UNIVERSITY代碼\modelcs_v1.1\Sparsity\cosamp_example.m

............此處省略202個文件信息