From 9577bafff44d10cd9c3eea28e321797de7b59d71 Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Wed, 20 Nov 2024 23:24:45 +0100 Subject: [PATCH] Latest working version - seemingly reproduces Blair Blakie results. --- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 1 + .../+Scripts/run_on_cluster.m | 20 +++++++++---------- .../bwhpc_matlab_gpe_sim_gpu.slurm | 2 +- 3 files changed, 11 insertions(+), 12 deletions(-) diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index c979301..a4bafa6 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -122,6 +122,7 @@ SaveDirectory = './Data_StripePhase'; JobNumber = 2; % Plotter.visualizeGSWavefunction2D(SaveDirectory, JobNumber) Scripts.analyzeRun2D(SaveDirectory, JobNumber) + %% - Cluster Runs - Analysis SaveDirectory = './Data_HoneycombPhase'; JobNumber = 3; diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m index bc5d723..3d409ab 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m @@ -9,7 +9,7 @@ % as = ((as/add)*Params.add)/Params.a0 % Critical point: 102.5133; Triangular phase: 98.0676; Stripe phase: 100.0289; Honeycomb phase: 101.9903 -%{ + %% - Create Variational2D and Calculator object with specified options OptionsStruct = struct; @@ -24,7 +24,7 @@ OptionsStruct.TrapPotentialType = 'None'; OptionsStruct.NumberOfGridPoints = [256, 256]; OptionsStruct.Dimensions = [100, 100]; -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 100E-6; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; @@ -34,8 +34,8 @@ OptionsStruct.NoiseScaleFactor = 4; OptionsStruct.MaxIterations = 20; OptionsStruct.VariationalWidth = 5.7; OptionsStruct.WidthLowerBound = 0.2; -OptionsStruct.WidthUpperBound = 12; -OptionsStruct.WidthCutoff = 1e-2; +OptionsStruct.WidthUpperBound = 20; +OptionsStruct.WidthCutoff = 1e-3; OptionsStruct.PlotLive = false; OptionsStruct.JobNumber = 1; @@ -51,7 +51,6 @@ solver.Potential = pot.trap(); %-% Run Solver %-% [Params, Transf, psi, V, VDk] = solver.run(); -%} %% - Create Variational2D and Calculator object with specified options @@ -94,7 +93,7 @@ solver.Potential = pot.trap(); %-% Run Solver %-% [Params, Transf, psi, V, VDk] = solver.run(); -%{ + %% - Create Variational2D and Calculator object with specified options OptionsStruct = struct; @@ -109,7 +108,7 @@ OptionsStruct.TrapPotentialType = 'None'; OptionsStruct.NumberOfGridPoints = [256, 256]; OptionsStruct.Dimensions = [100, 100]; -OptionsStruct.TimeStepSize = 500E-6; % in s +OptionsStruct.TimeStepSize = 100E-6; % in s OptionsStruct.MinimumTimeStepSize = 1E-5; % in s OptionsStruct.TimeCutOff = 2E6; % in s OptionsStruct.EnergyTolerance = 5E-10; @@ -119,8 +118,8 @@ OptionsStruct.NoiseScaleFactor = 4; OptionsStruct.MaxIterations = 20; OptionsStruct.VariationalWidth = 5.7; OptionsStruct.WidthLowerBound = 0.2; -OptionsStruct.WidthUpperBound = 12; -OptionsStruct.WidthCutoff = 1e-2; +OptionsStruct.WidthUpperBound = 20; +OptionsStruct.WidthCutoff = 1e-3; OptionsStruct.PlotLive = false; OptionsStruct.JobNumber = 3; @@ -135,5 +134,4 @@ pot = VariationalSolver2D.Potentials(options{: solver.Potential = pot.trap(); %-% Run Solver %-% -[Params, Transf, psi, V, VDk] = solver.run(); -%} \ No newline at end of file +[Params, Transf, psi, V, VDk] = solver.run(); \ No newline at end of file diff --git a/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm b/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm index cc582cc..3bf3d60 100644 --- a/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm +++ b/Dipolar-Gas-Simulator/bwhpc_matlab_gpe_sim_gpu.slurm @@ -8,7 +8,7 @@ #SBATCH --gres=gpu:A100:1 #SBATCH --mem=8G # Estimated wallclock time for job -#SBATCH --time=5:00:00 +#SBATCH --time=16:00:00 #SBATCH --job-name=simulation #SBATCH --error=simulation.err #SBATCH --output=simulation.out