Commit 8153b08b authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

If unspecified, false() produces a square matrix, fixed it in nimask

parent 0239a660

scsearch.m

+2 −2
Original line number Diff line number Diff line
@@ -345,7 +345,7 @@ parbank = combinations(porb_gr,a_gr,tasc_gr).Variables; 
toc





nimask = false(length(nibank));

nimask = false(length(nibank),1);

totlam = zeros(length(parbank),length(f_gr));

tic

% for n=1:length(f_gr)
@@ -448,7 +448,7 @@ sigmastar = 2*gammaincinv((1-(1-pfa)^(1/(sum(nimask)*length(f_gr)))),M,'upper'); 
[bru,tto] = max(totlam);

[maxtotlam,maxinde] = max(bru);

bestpar = [f_gr(maxinde), parbank(tto(maxinde),:), maxtotlam]; 

bparch(:) = num2str(bestpar(:));

bparch = string(bestpar);

disp('Our most significant peak has parameters:');

disp(['f_spin = '+bparch(1)+' Hz, p_orb = '+bparch(2)+' s, a*sin(i)/c = '+bparch(3)+' s, t_asc = '+bparch(4)+' s.']);

disp(['Its detection statistics is Σ_max = ',num2str(maxtotlam),',']);