From 443e533e31309dd0750771a3915a8180fe67ae53 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Tue, 20 May 2025 17:40:09 +0200 Subject: [PATCH] Latest script. --- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 33 +++++++++++++++++++- 1 file changed, 32 insertions(+), 1 deletion(-) diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index d8660c4..428fcd0 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -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 +