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

Fixed issue with creation of save directories.

Browse Source
This commit is contained in:
Karthik 2025-04-07 21:44:03 +02:00
parent 178350bdc2
commit 5053e3ffba
2 changed files with 9 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -82,6 +82,7 @@ OptionsStruct.PlotLive = true;
OptionsStruct.JobNumber = 0; OptionsStruct.JobNumber = 0;
OptionsStruct.RunOnGPU = false; OptionsStruct.RunOnGPU = false;
OptionsStruct.SaveData = false; OptionsStruct.SaveData = false;
OptionsStruct.SaveDirectory = './Results/Data_3D/RealTimeDynamics';
options = Helper.convertstruct2cell(OptionsStruct); options = Helper.convertstruct2cell(OptionsStruct);
sim = Simulator.DipolarGas(options{:}); sim = Simulator.DipolarGas(options{:});
@ -91,7 +92,8 @@ sim.Potential = pot.trap();
%-% Run Simulation %-% %-% Run Simulation %-%
[Params, Transf, psi, V, VDk] = sim.run(); [Params, Transf, psi, V, VDk] = sim.run();
SaveDirectory = './Results/Data_3D/RealTimeDynamics'; % Save only final state as initial wavefunction for real-time propagation
mkdir(sprintf(OptionsStruct.SaveDirectory))
save(sprintf(strcat(SaveDirectory, '/psi_init.mat'),Params.njob),'psi','Transf','Params','VDk','V'); save(sprintf(strcat(SaveDirectory, '/psi_init.mat'),Params.njob),'psi','Transf','Params','VDk','V');
close all close all
clear clear

13
Dipolar-Gas-Simulator/+Simulator/@DipolarGas/run.m
View File

@ -3,17 +3,12 @@ function [Params, Transf, psi,V,VDk] = run(this)
[Params] = this.setupParameters(); [Params] = this.setupParameters();
this.SimulationParameters = Params; this.SimulationParameters = Params;
if strcmp(this.SimulationMode, 'RealTimeEvolution')
this.QuenchSettings = this.setupQuenchSettings(Params);
end
% --- Set up spatial grids and transforms --- % --- Set up spatial grids and transforms ---
[Transf] = this.setupSpace(Params); [Transf] = this.setupSpace(Params);
[TransfRad] = this.setupSpaceRadial(Params); [TransfRad] = this.setupSpaceRadial(Params);
% --- Initialize --- % --- Initialize ---
mkdir(sprintf(this.SaveDirectory)) mkdir(sprintf(this.SaveDirectory))
[psi,V,VDk] = this.initialize(Params,Transf,TransfRad); [psi,V,VDk] = this.initialize(Params,Transf,TransfRad);
Observ.EVec = []; Observ.NormVec = []; Observ.PCVec = []; Observ.tVecPlot = []; Observ.mucVec = []; Observ.EVec = []; Observ.NormVec = []; Observ.PCVec = []; Observ.tVecPlot = []; Observ.mucVec = [];
@ -21,11 +16,15 @@ function [Params, Transf, psi,V,VDk] = run(this)
Observ.res_idx = 1; Observ.res_idx = 1;
% --- Run Simulation --- % --- Run Simulation ---
mkdir(sprintf(this.SaveDirectory)) if this.DoSave
mkdir(sprintf(strcat(this.SaveDirectory, '/Run_%03i'),Params.njob)) mkdir(sprintf(strcat(this.SaveDirectory, '/Run_%03i'),Params.njob))
end
if strcmp(this.SimulationMode, 'EnergyMinimization') if strcmp(this.SimulationMode, 'EnergyMinimization')
[psi] = this.runGradientDescent(psi,Params,Transf,VDk,V,Observ); [psi] = this.runGradientDescent(psi,Params,Transf,VDk,V,Observ);
else else
if strcmp(this.SimulationMode, 'RealTimeEvolution')
this.QuenchSettings = this.setupQuenchSettings(Params);
end
[psi] = this.propagateWavefunction(psi,Params,Transf,VDk,V,t_idx,Observ); [psi] = this.propagateWavefunction(psi,Params,Transf,VDk,V,t_idx,Observ);
end end
end end