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

Minor corrections

Browse Source
This commit is contained in:
Karthik 2025-05-15 18:01:39 +02:00
parent 40915f3913
commit 883a5bc453
2 changed files with 13 additions and 14 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Dipolar-Gas-Simulator/+Scripts/editPhaseDiagram.m
View File

@ -22,8 +22,8 @@ function editPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST)
yticks(1:ylen);
xticklabels(string(NUM_ATOMS_LIST));
yticklabels(string(SCATTERING_LENGTH_RANGE));
xlabel('Number of Atoms');
ylabel('Scattering Length (a₀)');
xlabel('Number of Atoms', 'Interpreter', 'tex', 'FontSize', 16);
ylabel('Scattering Length (\times a_o)', 'Interpreter', 'tex', 'FontSize', 16);
grid on;
title('Click a cell to modify value. Press "q" to quit.');

23
Dipolar-Gas-Simulator/+Scripts/run_locally.m
View File

@ -613,8 +613,7 @@ for i = 1:length(SCATTERING_LENGTH_RANGE)
end
%% Visualize phase diagram
% load('./Results/phase_diagram_matrix.mat')
load('./Results/edited_phase_diagram.mat')
load('./Results/PhaseDiagramUntilted.mat')
PhaseDiagramMatrix = M;
ScatteringLengths = SCATTERING_LENGTH_RANGE;
NumberOfAtoms = round(NUM_ATOMS_LIST * 1E-5,2);
@ -623,32 +622,32 @@ ylen = length(ScatteringLengths);
figure(1);
set(gcf,'Position',[100 100 950 750])
fig.WindowState = 'maximized';
clf
set(gcf,'Position', [100, 100, 1600, 900])
set(gca,'FontSize',16,'Box','On','Linewidth',2);
hImg = imagesc(M);
set(gca, 'YDir', 'normal');
colormap(parula);
ax = gca; % Get current axes
ax.FontSize = 16; % Set tick label font size (adjust as needed)
colorbar;
axis equal tight;
xticks(1:xlen);
yticks(1:ylen);
xticklabels(string(NumberOfAtoms));
yticklabels(string(ScatteringLengths));
xlabel('Number of Atoms (x 1E5)', FontSize=16);
ylabel('Scattering Length (a0)', FontSize=16);
% title('Zero-temperature Phase Diagram for \theta = 0', 'Interpreter', 'tex', FontSize=16);
xticklabels(string(NUM_ATOMS_LIST));
yticklabels(string(SCATTERING_LENGTH_RANGE));
xlabel('Number of Atoms', 'Interpreter', 'tex', 'FontSize', 16);
ylabel('Scattering Length (\times a_o)', 'Interpreter', 'tex', 'FontSize', 16);
grid on;
%% Edit phase diagram
% load('./Results/phase_diagram_matrix.mat')
load('./Results/edited_phase_diagram.mat')
load('./Results/PhaseDiagramUntilted.mat')
PhaseDiagramMatrix = M;
ScatteringLengths = SCATTERING_LENGTH_RANGE;
NumberOfAtoms = NUM_ATOMS_LIST;
Scripts.editPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST)
%% Smoothen phase diagram
load('./Results/edited_phase_diagram.mat'); % Load M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST
load('./Results/PhaseDiagramUntilted.mat'); % Load M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST
titleString = "[ \omega_x, \omega_y, \omega_z ] = 2 \pi \times [ 50, 20, 150 ] Hz; \theta = 0";
Scripts.plotSmoothedPhaseDiagram(M, SCATTERING_LENGTH_RANGE, NUM_ATOMS_LIST, titleString);