From 988df2b5f3e449bf49178cf42983121a360f32ce Mon Sep 17 00:00:00 2001 From: Karthik Chandrashekara Date: Fri, 22 Nov 2024 16:25:42 +0100 Subject: [PATCH] Added a line to save eigen-values, vectors to a mat file. --- .../+BdGSolver2D/solveBogoliubovdeGennesIn2D.m | 2 +- Dipolar-Gas-Simulator/+Scripts/run_locally.m | 3 +++ Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m | 3 +++ 3 files changed, 7 insertions(+), 1 deletion(-) diff --git a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m index 404d771..e68a59d 100644 --- a/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m +++ b/Dipolar-Gas-Simulator/+BdGSolver2D/solveBogoliubovdeGennesIn2D.m @@ -1,6 +1,6 @@ function [evals, modes] = solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem) -Size = length(psi(:)); +Size = length(psi(:)); % 2-D matrix will be unwrapped and the corresponding matrix of (N^2)^2 elements solved for. Neigs = length(psi(:)); opts.tol = 1e-16; opts.disp = 1; diff --git a/Dipolar-Gas-Simulator/+Scripts/run_locally.m b/Dipolar-Gas-Simulator/+Scripts/run_locally.m index 4f08fb3..683d5d0 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_locally.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_locally.m @@ -128,6 +128,9 @@ muchem = solver.Calculator.calculateChemicalPotential(psi,Params,VParams,Transf, [evals, modes] = BdGSolver2D.solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem); +% Save the eigenvalues and eigenvectors to a .mat file +save(sprintf(strcat(solver.SaveDirectory, '/Run_%03i/bdg_eigen_data.mat'),solver.JobNumber), 'evals', 'modes'); + %% - Create Variational2D and Calculator object with specified options OptionsStruct = struct; diff --git a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m index ff0cf5e..afc632c 100644 --- a/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m +++ b/Dipolar-Gas-Simulator/+Scripts/run_on_cluster.m @@ -204,3 +204,6 @@ V = 0.0*(Params.gx.*X.^2+Params.gy.*Y.^2); muchem = solver.Calculator.calculateChemicalPotential(psi,Params,VParams,Transf,VDk,V); [evals, modes] = BdGSolver2D.solveBogoliubovdeGennesIn2D(psi, Params, VDk, VParams, Transf, muchem); + +% Save the eigenvalues and eigenvectors to a .mat file +save(sprintf(strcat(solver.SaveDirectory, '/Run_%03i/bdg_eigen_data.mat'),Params.njob), 'evals', 'modes');