Homec4science
Diffusion COSOlver (new_coulomb) Browse

cosolver/matlab_scripts/check_GKCoulomb.m
b3f94e22e87cnew_coulomb

Last Change
Raw File

check_GKCoulomb.m
View Options

%% Check GK Coulomb
%% Clean up current workspace
clear all; close all;
%% Load file
dir = '/home/bjfrei/Documents/MoliSolver/COSOlver';
%filename_Coulomb1 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_6_JE_12_PMMAXX_16_NFLR_1_kperp_0.0500.h5';
%filename_Coulomb1 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_3_Jmaxi_3_JE_4_PMMAXX_13_NFLR_2_kperp_0.8000.h5';
file = 'iCa/self_Coll_GKE_0_GKI_0_ESELF_0_ISELF_7_Pmaxe_-1_Jmaxe_-1_Pmaxi_80_Jmaxi_20_JE_20.h5'
% Load collision matrix
Cpj = h5read(fullfile(dir,file),'/Caapj/Ciipj');
% Load collision matrices
%C_Coulomb1 = h5read(fullfile(dir,filename_Coulomb1),'/Caapj/Ciipj');
%C_CoulombJs = sum(sum(C_Coulomb1(:,:,:,:,:),5),4);
%idx = [1 2 3 4 5 6 7 8 9 10 11 12 13];
%C_Coulomb1 = sum(C_CoulombJs(:,:,idx),3);
%C_Coulomb1(C_Coulomb1 == 0) = NaN;
figure;
pcolor(Cpj(1:100,1:100)); axis ij;
colorbar
%% Eigenvalues of collisional matrices
% Coulomb eig
[V,D] = eig(Cpj);
eigCpj = diag(Cpj);
% Coulomb eig
%[V,D] = eig(C_Abel);
%eig_Abel = diag(D);
figure; hold all;
plot(real(eigCpj),imag(eigCpj),'b*')
%% Load *.h5 file
dir = '/home/bjfrei/Documents/MoliSolver/COSOlver/iCa/';
f0 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_0_kperp_0.5000.h5';
f1 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_1_kperp_0.5000.h5';
f2 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_2_kperp_0.5000.h5';
f3 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_3_kperp_0.5000.h5';
f4 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_4_kperp_0.5000.h5';
f5 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_5_kperp_0.5000.h5';
f6 = 'self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_6_kperp_0.5000.h5';
% Load collision matrices
C0 = h5read(fullfile(dir,f0),'/Caapj/Ciipj');
C1 = h5read(fullfile(dir,f1),'/Caapj/Ciipj');
C2 = h5read(fullfile(dir,f2),'/Caapj/Ciipj');
C3 = h5read(fullfile(dir,f3),'/Caapj/Ciipj');
C4 = h5read(fullfile(dir,f4),'/Caapj/Ciipj');
C5 = h5read(fullfile(dir,f5),'/Caapj/Ciipj');
C6 = h5read(fullfile(dir,f6),'/Caapj/Ciipj');
%% Contract energy and FLR dimension of GKColl
C0 = sum( sum(sum(C0(:,:,:,:,:),5),4),3);
C1 = sum( sum(sum(C1(:,:,:,:,:),5),4),3);
C2 = sum( sum(sum(C2(:,:,:,:,:),5),4),3);
C3 = sum( sum(sum(C3(:,:,:,:,:),5),4),3);
C4 = sum( sum(sum(C4(:,:,:,:,:),5),4),3);
C5 = sum( sum(sum(C5(:,:,:,:,:),5),4),3);
C6 = sum( sum(sum(C6(:,:,:,:,:),5),4),3);
%% Inquire convergence
Pmax = 15;
Jmax = 4;
bari = @(pp,jj) (Jmax +1)*pp + 1 + jj ;
NFLR = 0:1:6;
Cpj= zeros(length(NFLR),1);
p =1;
j =0;
for iNFLR=1:length(NFLR)
fi = ['self_Coll_GKE_0_GKI_1_ESELF_0_ISELF_1_Pmaxe_-1_Jmaxe_-1_Pmaxi_15_Jmaxi_4_pamaxx_10_NFLR_',num2str(NFLR(iNFLR)),'_kperp_0.5000.h5'];
CiNFLR = h5read(fullfile(dir,fi),'/Caapj/Ciipj');
CiNFLR = sum( sum(sum(CiNFLR(:,:,:,:,:),5),4),3);
Cpj(iNFLR) = CiNFLR(bari(p,j),bari(p,j));
end
figure; hold all;
h1 = plot(NFLR,Cpj);
c4science · Help