Commit 5dc01cda authored by riccardo's avatar riccardo
Browse files

Optimizations?

parent 869ef806

scsearch.m

+6 −6
Original line number Diff line number Diff line
@@ -389,7 +389,7 @@ tm=gpuArray(tm); 
nibank = gpuArray(combinations(nis{:}).Variables);

% Lambda = zeros(length(f_gr),length(nibank),'single','gpuArray'); % MATLAB

% stores array elements in a column-major layout

Lambda = zeros(length(f_gr),length(nibank),'single','gpuArray');

Lambda = ones(length(f_gr),length(nibank),'single','gpuArray');

F1=gpuArray((0:N-1)./(N*dt_psd)).';

f_ind = zeros(1,length(f_gr),'uint32');

for n = 1:length(f_gr)
@@ -448,7 +448,7 @@ for m=1:M 
    Y0 = fft(xtemp);

    % norman = 2./mean((abs(Y0(F1>1000 & F1<Nyq)).^2)./sumx);

    norman = 2./(Nphot.*var(xtemp)./mean(xtemp));

    tau = ones("like",ttemp);

    tau = ones(size(ttemp),"like",ttemp);





    % toc
@@ -479,8 +479,8 @@ for m=1:M 
        % pause

        % Careful! ni_0 still needs to be defined (it's the spin freq. 

        % currently being considered)

        Yenne(:)=Y1(f_ind(:)); 

        Lambda(:,i) = (abs(Yenne(:)));

        % Yenne(:)=Y1(f_ind(:)); 

        Lambda(:,i) = Y1(f_ind(:)); 

        % Lambda(:,i) = norman.*(abs(Yenne(:)).^2)./sumo;

        % 

        % for n = 1:length(f_gr)
@@ -491,7 +491,7 @@ for m=1:M 
        % end

        % disp('1ni')

    end

    Lambda = (Lambda.^2).*norman;

    Lambda = (abs(Lambda).^2).*norman;

    % disp(max(Lambdacp(:,:)-Lambda(:,:)))

    % disp(min(Lambdacp(:,:)-Lambda(:,:)))

    % toc
@@ -500,7 +500,7 @@ for m=1:M 


    lamfiname = sprintf('Lambda_fmax_%d_seg_%d.mat', f_max, m);

    save([pathfi,lamfiname],"Lambda","-v7.3","-nocompression");

    Lambda = zeros(length(f_gr),length(nibank),'single');

    Lambda = ones(length(f_gr),length(nibank),'single','gpuArray');

    % Lambda = zeros(length(f_gr),length(nibank),'single','gpuArray');

    segend = toc(segstart);

    avetime = avetime + segend;