Homec4science
Diffusion GYACOMO (OBSOLETE: LATEST VERSION ON gitlab.epfl.ch) (master) Browse

gyacomo/matlab/plot/plot_cosol_mat.m
0597cb3c91dcmaster

Last Change
Raw File

plot_cosol_mat.m
View Options

% MAT = results.iCa;
% figure
% suptitle('FCGK,P=6,J=3');
% subplot(221)
% imagesc(log(abs(MAT)));
% title('log abs')
% subplot(222)
% imagesc(imag(MAT)); colormap(bluewhitered)
% title('imag'); colorbar
% subplot(223)
% imagesc(imag(MAT)>0);
% title('imag$>$0');
% subplot(224)
% imagesc(imag(MAT)<0);
% title('imag$<$0');
%
%
%% SGGK
P_ = 20; J_ = 10;
mat_file_name = '/home/ahoffman/HeLaZ/iCa/gk_sugama_P_20_J_10_N_150_kpm_8.0.h5';
SGGK_self = h5read(mat_file_name,'/00000/Caapj/Ciipj');
SGGK_ei = h5read(mat_file_name,'/00000/Ceipj/CeipjF')+h5read(mat_file_name,'/00000/Ceipj/CeipjT');
SGGK_ie = h5read(mat_file_name,'/00000/Ciepj/CiepjF')+h5read(mat_file_name,'/00000/Ciepj/CiepjT');
figure
MAT = 1i*SGGK_self; suptitle('SGGK ii,P=20,J=10, k=0');
% MAT = 1i*SGGK_ei; suptitle('SGGK ei,P=20,J=10');
% MAT = 1i*SGGK_ie; suptitle('SGGK ie,P=20,J=10');
subplot(221)
imagesc(abs(MAT));
title('log abs')
subplot(222)
imagesc(imag(MAT)); colormap(bluewhitered)
title('imag'); colorbar
subplot(223)
imagesc(imag(MAT)>0);
title('imag$>$0');
subplot(224)
imagesc(imag(MAT)<0);
title('imag$<$0');
%% PAGK
P_ = 20; J_ = 10;
mat_file_name = '/home/ahoffman/HeLaZ/iCa/gk_pitchangle_8_P_20_J_10_N_150_kpm_8.0.h5';
PAGK_self = h5read(mat_file_name,'/00000/Caapj/Ceepj');
% PAGK_ei = h5read(mat_file_name,'/00000/Ceipj/CeipjF')+h5read(mat_file_name,'/00000/Ceipj/CeipjT');
% PAGK_ie = h5read(mat_file_name,'/00000/Ciepj/CiepjF')+h5read(mat_file_name,'/00000/Ciepj/CiepjT');
figure
MAT = 1i*PAGK_self; suptitle('PAGK ii,P=20,J=10, k=0');
subplot(221)
imagesc(abs(MAT));
title('log abs')
subplot(222)
imagesc(imag(MAT)); colormap(bluewhitered)
title('imag'); colorbar
subplot(223)
imagesc(imag(MAT)>0);
title('imag$>$0');
subplot(224)
imagesc(imag(MAT)<0);
title('imag$<$0');
%% FCGK
P_ = 4; J_ = 2;
mat_file_name = '/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_12_P_4_J_2_N_50_kpm_4.0.h5';
% mat_file_name = '/home/ahoffman/gyacomo/iCa/LDGK_P6_J3_dk_5e-2_km_2.5_NFLR_20.h5';
% mat_file_name = '/home/ahoffman/gyacomo/iCa/LDGK_P10_J5_dk_5e-2_km_5_NFLR_12.h5';
kp = 0.0;
kp_a = h5read(mat_file_name,'/coordkperp');
[~,matidx] = min(abs(kp_a-kp));
matidx = sprintf('%5.5i',matidx);
FCGK_self = h5read(mat_file_name,['/',matidx,'/Caapj/Ciipj']);
FCGK_ei = h5read(mat_file_name,['/',matidx,'/Ceipj/CeipjF'])+h5read(mat_file_name,'/00000/Ceipj/CeipjT');
FCGK_ie = h5read(mat_file_name,['/',matidx,'/Ciepj/CiepjF'])+h5read(mat_file_name,'/00000/Ciepj/CiepjT');
figure
MAT = 1i*FCGK_self; suptitle(['FCGK ii,P=6,J=3, k=',num2str(kp),' (',matidx,')']);
% MAT = 1i*FCGK_ei; suptitle('FCGK ei,P=20,J=10');
% MAT = 1i*FCGK_ie; suptitle('FCGK ie,P=20,J=10');
subplot(221)
imagesc(abs(MAT));
title('log abs')
subplot(222)
imagesc(imag(MAT)); colormap(bluewhitered)
title('imag'); colorbar
subplot(223)
imagesc(imag(MAT)>0);
title('imag$>$0');
subplot(224)
imagesc(imag(MAT)<0);
title('imag$<$0');
%% Eigenvalue spectrum analysis
if 0
%%
mfns = {...
'/home/ahoffman/gyacomo/iCa/gk_sugama_P_20_J_10_N_150_kpm_8.0.h5',...
% '/home/ahoffman/gyacomo/iCa/gk_pitchangle_8_P_20_J_10_N_150_kpm_8.0.h5',...
% '/home/ahoffman/gyacomo/iCa/LDGKii_P10_J5_dk_5e-2_km_5_NFLR_4.h5',...
% '/home/ahoffman/gyacomo/iCa/LDGKii_P10_J5_dk_5e-2_km_5_NFLR_12.h5',...
% '/home/ahoffman/gyacomo/iCa/LDGKii_P10_J5_dk_5e-2_km_5_NFLR_12_k2trunc.h5',...
% '/home/ahoffman/gyacomo/iCa/LDGKii_P10_J5_dk_5e-2_km_5_NFLR_30.h5',...
'/home/ahoffman/gyacomo/iCa/LDGK_P6_J3_dk_5e-2_km_2.5_NFLR_20.h5',...
% '/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_6_P_4_J_2_N_50_kpm_4.0.h5',...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_12_P_4_J_2_N_50_kpm_4.0.h5',...
% '/home/ahoffman/gyacomo/iCa/gk.hacked_sugama_P_10_J_5_N_150_kpm_8.0.h5',...
% '/home/ahoffman/gyacomo/iCa/gk.hacked_sugama_P_4_J_2_N_75_kpm_5.0.h5',...
};
CONAME_A = {...
'SG 20 10',...
% 'PA 20 10',...
% 'FC 10 5 NFLR 4',...
% 'FC 10 5 NFLR 12',...
% 'FC 10 5 NFLR 12 k<2', ...
'LD 6 3 NFLR 20', ...
% 'FC 4 2 NFLR 6',...
'FC 4 2 NFLR 12', ...
% 'Hacked SG A',...
% 'Hacked SG B',...
};
figure
for j_ = 1:numel(mfns)
mat_file_name = mfns{j_}; disp(mat_file_name);
kp_a = h5read(mat_file_name,'/coordkperp');
% kp_a = kp_a(kp_a<=3);
gmax = zeros(size(kp_a));
wmax = zeros(size(kp_a));
for idx_ = 0:numel(kp_a)-1
matidx = sprintf('%5.5i',idx_);
MAT = h5read(mat_file_name,['/',matidx,'/Caapj/Ciipj']);
gmax(idx_+1) = max(abs(real(eig(MAT))));
wmax(idx_+1) = max(abs(imag(eig(MAT))));
end
subplot(121)
plot(kp_a,gmax,'DisplayName',CONAME_A{j_}); hold on;
subplot(122)
plot(kp_a,wmax,'DisplayName',CONAME_A{j_}); hold on;
end
subplot(121)
legend('show'); grid on;
ylim([0,100]);
xlabel('$k_\perp$'); ylabel('$\gamma_{max}$ from Eig(iCa)')
subplot(122)
legend('show'); grid on;
xlabel('$k_\perp$'); ylabel('$\omega_{max}$ from Eig(iCa)')
end
%% Van Kampen plot
if 0
%%
kperp= 1.5;
mfns = {'/home/ahoffman/gyacomo/iCa/gk_sugama_P_20_J_10_N_150_kpm_8.0.h5',...
'/home/ahoffman/gyacomo/iCa/gk_pitchangle_8_P_20_J_10_N_150_kpm_8.0.h5',...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_6_P_4_J_2_N_50_kpm_4.0.h5',...
'/home/ahoffman/gyacomo/iCa/gk_coulomb_NFLR_12_P_4_J_2_N_50_kpm_4.0.h5',...
'/home/ahoffman/gyacomo/iCa/LDGK_P6_J3_dk_5e-2_km_2.5_NFLR_20.h5',...
};
CONAME_A = {'SG 20 10','PA 20 10', 'FC 4 2 NFLR 6', 'FC 4 2 NFLR 12', 'LD 6 3 NFLR 12'};
grow = {};
puls = {};
for j_ = 1:numel(mfns)
mat_file_name = mfns{j_}; disp(mat_file_name);
kp_a = h5read(mat_file_name,'/coordkperp');
[~,idx_] = min(abs(kp_a-kperp));
matidx = sprintf('%5.5i',idx_);
MAT = h5read(mat_file_name,['/',matidx,'/Caapj/Ciipj']);
grow{j_} = real(eig(MAT));
puls{j_} = imag(eig(MAT));
end
figure
for j_ = 1:numel(mfns)
% plot(puls{j_}, grow{j_},'o','DisplayName',CONAME_A{j_}); hold on;
plot(grow{j_},'o','DisplayName',CONAME_A{j_}); hold on;
end
legend('show'); grid on; title(['$k_\perp=$',num2str(kperp)]);
xlabel('$\omega$ from Eig(iCa)'); ylabel('$\gamma$ from Eig(iCa)')
end
c4science · Help