@@ -35,19 +35,19 @@ plot!(left_margin = 5mm) # hide
3535## analyze the signal via HGLET
3636dmatrixH, dmatrixHrw, dmatrixHsym = HGLET_Analysis_All(G, GP)
3737dvec_hglet, BS_hglet, trans_hglet = HGLET_GHWT_BestBasis(GP, dmatrixH = dmatrixH,
38- dmatrixHrw = dmatrixHrw, dmatrixHsym = dmatrixHsym, costfun = 1)
38+ dmatrixHrw = dmatrixHrw, dmatrixHsym = dmatrixHsym, cfspec = 1)
3939
4040## LP-HGLET
4141dmatrixsH, dmatrixsHsym = LPHGLET_Analysis_All(G, GP; ϵ = 0.3)
4242dvec_lphglet, BS_lphglet, trans_lphglet = HGLET_GHWT_BestBasis(GP, dmatrixH = dmatrixsH,
43- dmatrixHsym = dmatrixsHsym, costfun = 1)
43+ dmatrixHsym = dmatrixsHsym, cfspec = 1)
4444
4545# find the top 10 HGLET basis vectors
4646important_idx = sortperm(dvec_hglet[:].^2; rev = true)
4747hglet_top10 = zeros(N, 10)
4848for i in 1:10
4949 w, _ = HGLET_Synthesis(reshape(spike(important_idx[i], N), (N, 1)),
50- GP, BS_hglet, G, method = :L)
50+ GP, BS_hglet, G, gltype = :L)
5151 hglet_top10[:, i] = w[:]
5252end
5353wiggle(hglet_top10; sc = 0.45)
@@ -58,7 +58,7 @@ important_idx = sortperm(dvec_lphglet[:].^2; rev = true)
5858lphglet_top10 = zeros(N, 10)
5959for i in 1:10
6060 w, _ = LPHGLET_Synthesis(reshape(spike(important_idx[i], N), (N, 1)),
61- GP, BS_lphglet, G; method = :L, ϵ = 0.3)
61+ GP, BS_lphglet, G; gltype = :L, ϵ = 0.3)
6262 lphglet_top10[:, i] = w[:]
6363end
6464wiggle(lphglet_top10; sc = 0.45)
0 commit comments