Commit e8f86c31 authored by Riccardo La Placa's avatar Riccardo La Placa
Browse files

Checking username settings

parent 67f5a6a3

scsearch.m

+13 −8
Original line number Diff line number Diff line
@@ -6,15 +6,18 @@ pathfi = '/data/Sorgenti/3FGLJ1544.6-1125/'; 
reset(gpuDevice(2));

gpuDevice(2);

%%%%%%%%%%%%%%%%

daterella = ['Date and time now ',char(datetime('now')),' (UTC), maybe'];

disp(daterella); clear daterella



tic



finame = 'SiFAP2_20220429_3FGLJ1544_N1_Bary_N3_NOTCORRECTED.fits';

finame = 'SiFAP2_20220429_3FGLJ1544_N1_Bary.fits';

% finame = 'J1023_B_2017_Bary.fits';

% finame = 'EPN_0744840201_bary.fits';

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

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

t_raw = t_raw{1};

MJDREF=0.943055555602768;

MJDREF=59698.943055555602768;

% t_raw=t_raw./86400+50814;

% MJDREF=fix(t_raw(1));

% t_raw=(t_raw-MJDREF).*86400;
@@ -56,9 +59,9 @@ Tseg=490; %segments' length in seconds 
f_gr = (768:1/Tseg:780).';

%f_gr=(772:1/Tseg:774).';

% porb_gr = (18000:10:21600).';

porb_gr=((20868.72-3*0.31):0.1:(20868.72+3*0.31));

porb_gr=(20868.72).';

a_gr = (0.01:8.0e-4:0.26).';

tasc_gr = (((0.724:8e-4:0.922)-MJDREF).*86400).';

tasc_gr = ((MJDREF - (59698.7231:8e-4:59698.9232)).*86400).';

% for j = 1:5

%     %f_gr(j) = f_tru*((j^2)/9);

%     porb_gr(j) = porb_tru*((j^2)/9);
@@ -191,7 +194,7 @@ tic 


% Try s* and check \nu_s range

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

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

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

%s_s=uint8(4);

s_s = 4;

while(1)
@@ -255,6 +258,8 @@ toc 
% dtau = zeros(1,N);

% taul = zeros(N+1,1);

taul = zeros(2,1);

disp(length(nibank));

disp(s_s);

disp(M);

for m=1:M
@@ -315,7 +320,7 @@ for m=1:M 
    lamfiname = sprintf('Lambda_seg_%d.mat', m);

    % save(lamfiname,"Lambda");

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

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

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

    disp(m);



end
@@ -433,7 +438,7 @@ toc 


% disp(bestpar);

% save('C:\Users\Filippo\Desktop\XMM_Jxxx\risultelli.mat');

save([pathfi,'risultelli_07_03_',num2str(Tseg),'s.mat'],"-v7.3");

save([pathfi,'risultelli_',char(datetime('now','Format','dd_MM')),'_',num2str(Tseg),'s.mat'],"-v7.3");

%%