Fixed nismax nismin

    % parameters over their respective ranges (given in Table 2). This

    % is done with the exception of ν which is held fixed at its

    % maximum value within sub-bands over the frequency search space.

    if singal(s)>0

        nismin(s) = f_min*a_min*(Omega_min^s)*singal(s);

        nismax(s) = f_max*a_max*(Omega_max^s)*singah(s);

    elseif singah(s)>0

        nismin(s) = f_max*a_max*(Omega_max^s)*singal(s);

        nismax(s) = f_max*a_max*(Omega_max^s)*singah(s);

    else

        nismin(s) = f_max*a_max*(Omega_max^s)*singal(s);

        nismax(s) = f_min*a_min*(Omega_min^s)*singah(s);

    end

end

toc

% Try s* and check \nu_s range

g_jj=((pi*Tseg)^2)/3.*[1; (Tseg^2)/60; (Tseg^3)/1344; (Tseg^4)/172800]; %eq. 22 M2015 + calcoli da eq. 21 M2015

mu_s=0.5; %massimo mismatch sulla griglia coerente da scegliere

mu_s=0.1; %massimo mismatch sulla griglia coerente da scegliere

%% Troppo alto μ ma non abbiamo RAM adesso

%% Magari spezzare a blocchi i nibank e poi cercare la minima distanza 

%% in ciascuno?

    end

end

%%  Create (...) a coherent template bank in the ν_s space

%%  Create a coherent template bank in the ν_s space

%%  Save this template bank, we'll need it later

% Nni+1 is the number of ni_s points in the grid (counting borders, 

% Nni is the number of steps) for each s couple.

        tau = sum((nibank(i,1:s_s)./(nizero*factorial(1:s_s))).*((tm(m,1:N)-tmid(m)).^((1:s_s).').'),2);

        X1 = interp1(tm(m,:),x(m,:),tau,'linear',0);