test_hottensor_jd

 script test_hottensor_jd

0001 % script test_hottensor_jd
0002 %
0003 clear all
0004 close all
0005 %
0006 var = 'beta';
0007 %
0008 base.a = 0.1;
0009 base.y = 0.55;
0010 base.beta = 0.1;
0011 base.Npar = 0.3;
0012 %
0013 nx = 11;
0014 %
0015 list.a = linspace(0,1,nx);
0016 list.y = linspace(0,1,nx);
0017 list.beta = [5e-3,linspace(0.1/(nx-1),0.1,(nx-1))];
0018 list.Npar = linspace(0,1,nx);
0019 %
0020 zZi = [1,1,1];
0021 zmi = [1,2,3];
0022 zfi = [0,1,0];
0023 %
0024 [qe,me,mp] = pc_dke_yp;
0025 %
0026 x = list.(var);
0027 %
0028 Nperpp = zeros(1,nx);
0029 Nperpm = zeros(1,nx);
0030 X_cold = zeros(3,3,nx);
0031 X_hotp = zeros(3,3,nx);
0032 X_hotm = zeros(3,3,nx);
0033 dXdNpar_hotp = zeros(3,3,nx);
0034 dXdNpar_hotm = zeros(3,3,nx);
0035 dXdNperp_hotp = zeros(3,3,nx);
0036 dXdNperp_hotm = zeros(3,3,nx);
0037 %
0038 for ix = 1:nx,
0039     %
0040     calc = base;
0041     calc.(var) = x(ix);
0042     %
0043     calc.sa = calc.a*[1;(zZi.^2.*zfi./zmi/(mp/me)).'];
0044     calc.sy = calc.y*[-1;(zZi./zmi/(mp/me)).'];
0045     calc.sbeta = calc.beta*[1;([1,1,1]).'];
0046     %
0047     [Nperpp(ix),Nperpm(ix),Kxyz] = coldtensor_jd(calc.sa,calc.sy,calc.Npar);
0048     %
0049     X_cold(:,:,ix) = Kxyz - eye(3);
0050     %
0051     [X_hotp(:,:,ix),dXdNpar_hotp(:,:,ix),dXdNperp_hotp(:,:,ix)] = hottensor_jd(calc.sa,calc.sy,calc.sbeta,calc.Npar,Nperpp(ix),1);
0052     [X_hotm(:,:,ix),dXdNpar_hotm(:,:,ix),dXdNperp_hotm(:,:,ix)] = hottensor_jd(calc.sa,calc.sy,calc.sbeta,calc.Npar,Nperpm(ix),1);
0053     %
0054 end
0055     
0056     
0057     
0058     
0059     
0060     
0061     
0062     
0063     
0064