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

Latest script.

Browse Source
This commit is contained in:
Karthik 2025-05-20 17:40:09 +02:00
parent 8ae97eb267
commit 443e533e31
1 changed files with 32 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -557,7 +557,7 @@ SaveDirectory = './Results/Data_Full3D/PhaseDiagram/ImagTimePropagation/aS_8.312
JobNumber = 0;
Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
%%
SaveDirectory = 'D:/Results - Numerics/Data_Full3D/PhaseDiagram/ImagTimePropagation/Theta0/HighN/aS_8.729000e+01_theta_000_phi_000_N_100000';
SaveDirectory = 'D:/Results - Numerics/Data_Full3D/PhaseDiagram/ImagTimePropagation/Theta20/aS_8.833000e+01_theta_020_phi_000_N_2550000';
JobNumber = 0;
Plotter.visualizeGSWavefunction(SaveDirectory, JobNumber)
%% Identify and count droplets
@ -686,6 +686,36 @@ for i = 1:length(SCATTERING_LENGTH_RANGE)
pause(1.5)
end
end
%% Missing files
data = {
'aS_8.208000e+01_theta_000_phi_000_N_90909'
'aS_8.729000e+01_theta_000_phi_000_N_77273'
'aS_8.729000e+01_theta_000_phi_000_N_95455'
'aS_8.729000e+01_theta_000_phi_000_N_3775000'
'aS_9.250000e+01_theta_000_phi_000_N_5000000'
};
% Preallocate arrays
numEntries = numel(data);
aS_values = zeros(numEntries, 1);
N_values = zeros(numEntries, 1);
% Regular expression pattern
pattern = 'aS_([\d.eE+-]+)_.*?_N_(\d+)';
% Extract values using regular expressions
for i = 1:numEntries
tokens = regexp(data{i}, pattern, 'tokens');
if ~isempty(tokens)
aS_values(i) = str2double(tokens{1}{1});
N_values(i) = str2double(tokens{1}{2});
end
end
% Create and display table
T = table(aS_values, N_values, 'VariableNames', {'aS', 'N'});
disp(T);
%% Visualize phase diagram
load('./Results/Data_Full3D/PhaseDiagramUntilted.mat')
@ -761,3 +791,4 @@ else
end