karthik/Calculations
1
0
Fork 0
Code Issues Pull Requests repo.project_board Releases Wiki Activity

Removed unnecessary function.

Browse Source
This commit is contained in:
Karthik 2024-11-15 19:13:24 +01:00
parent 927a6642da
commit c039356cc2
1 changed files with 0 additions and 26 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

26
Dipolar-Gas-Simulator/+VariationalSolver2D/@Calculator/find_nk_fit.m
View File

@ -1,26 +0,0 @@
function [ret] = find_nk_fit(~,VParams)
Lz = 20;
Nz = 4000;
normfac = Lz/Nz;
z = linspace(-0.5*Lz,0.5*Lz-Lz/Nz,Nz);
dz = abs(z(2)-z(1));
Kmax = pi*Nz/Lz;
kz = linspace(-Kmax,Kmax,Nz+1);
kz = kz(1:end-1);
psiz = exp(-0.5*(abs(z)/VParams.ell).^VParams.nu);% VParams.nu
normz = sum(abs(psiz).^2)*dz;
psiz = psiz / sqrt(normz);
psik = fftn(psiz)*normfac/(sqrt(2*pi)); % fftshifted
nk = abs(psik).^2;
nk_forfit = ifftshift(nk);
f = fit(kz',nk_forfit','gauss1'); % f(x) = a1*exp(-((x-b1)/c1)^2)
ret.ell = VParams.ell;
ret.nu = VParams.nu;
ret.ell_eff = 1/f.c1;
end