nothing

% finame = 'J1023_B_2017_Bary.fits';

finame = 'EPN_0744840201_bary.fits';

% t_raw = fitsread([pathfi,finame],"binarytable");

t_raw = fitsread('C:\Users\Filippo\Desktop\EPN_0744840201_bary.fits','binarytable');

t_raw = fitsread([pathfi,finame],"binarytable");

% t_raw = fitsread('C:\Users\Filippo\Desktop\EPN_0744840201_bary.fits','binarytable');

t_raw = t_raw{1};

t_raw=t_raw./86400+50814;

MJDREF=t_raw(1);

toc

tic

Tseg=128; %segments' length in seconds

Tseg=256; %segments' length in seconds

M=fix((t(end)-t(1))/Tseg); %number of segments

%con fix prendo la parte intera, scarto l'ultimo segmento che tanto non

%sarà mai di lunghezza Tseg (molto improbabile)

% % M2015 eq. 15

% singal = zeros(4,1);

% singah = zeros(4,1);

% nismin = zeros(4,1);

% nismax = zeros(4,1);

nismin = zeros(4,1);

nismax = zeros(4,1);

% for s = 1:4

%     singah(s) = max(sin(gam + s*pi/2),[],'all');

%     singal(s) = min(sin(gam + s*pi/2),[],'all');

% Cercheremo un modo per combinare in tutti i possibili modi vettori di 

% vari ni1,ni2,...,nis_s

%%%[kung,fu,fight] = ndgrid(nis{1},nino{2},nino{3});

tm=gpuArray(tm);

nibank = gpuArray(combinations(nis{:}).Variables);

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

% tm=gpuArray(tm);

nibank = combinations(nis{:}).Variables;

Lambda = zeros(length(nibank),length(f_gr),M);

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

toc 

%Fourier transform on original time-series --------------------------------

%per mantenere l'informazione di fase, non faccio il valore assoluto al quadrato della fft

            Lambda(i,n,m)=sum(abs(Y1).^2)/sum(x(m,:)); %CREDO (oppure prendono la potenza massima?)

            % toc

        end

        disp('1ni')

        % disp('1ni')

    end

    toc