Merge branch 'master' of github.com:grasingerm/polymer-stats

grasingerm · Oct 29, 2021 · 431d5a9 · 431d5a9
2 parents 6337812 + b0c5f26
commit 431d5a9
Show file tree

Hide file tree

Showing 3 changed files with 149 additions and 0 deletions.
diff --git a/2D/run/Ising_2021-10-28.jl b/2D/run/Ising_2021-10-28.jl
@@ -0,0 +1,50 @@
+println("Hello world!");
+
+@everywhere using Printf;
+
+@show workdir = if length(ARGS) > 0
+  ARGS[1];
+else
+  ".";
+end
+
+@everywhere fmt(x) = @sprintf("%07d", round(Int, 1e3*x));
+@everywhere fmt_int(x) = @sprintf("%03d", x);
+
+@everywhere function prefix(case)
+  "E0-$(fmt(case[:E0]))_K1-$(fmt(case[:K1]))_K2-$(fmt(case[:K2]))_kT-$(fmt(case[:kT]))_Fz-$(fmt(case[:Fz]))_Fx-$(fmt(case[:Fx]))_n-$(fmt(case[:n]))_b-$(fmt(case[:b]))";
+end
+
+cases = Any[];
+E0s = [1e-1; 0.5; 1.0];
+K1s = [0.0]; 
+K2s = [1.0]; 
+kTs = vcat(0.25:0.25:1.0, 2.0:1.0:5.0, 10.0:10.0:100.0);
+Fzs = [1.0; 5.0; 10.0; 15.0; 20.0; 25.0];
+Fxs = [0.0];
+ns = Int[100];
+bs = [0.5; 1.0; 2.0];
+#bs = [1.0];
+for b in bs, n in ns, Fx in Fxs, Fz in Fzs, kT in kTs, E0 in E0s, K1 in K1s, K2 in K2s
+  push!(cases, Dict(:E0 => E0, :K1 => K1, :K2 => K2,
+                    :kT => kT, :Fz => Fz,
+                    :Fx => Fx, :n => n, :b => b, :run => run));
+end
+
+@info "total number of cases to run: $(length(cases))";
+
+mkpath(workdir);
+
+pmap(case -> begin;
+
+  outfile = joinpath(workdir, "$(prefix(case)).out");
+  if !isfile(outfile)
+    println("Running case: $case.");
+    command = `julia -O 3 mcmc_clustering_eap_chain.jl --chain-type dielectric --energy-type Ising -b $(case[:b]) --E0 $(case[:E0]) --K1 $(case[:K1]) --K2 $(case[:K2]) --kT $(case[:kT]) --Fz $(case[:Fz]) --Fx $(case[:Fx]) -n $(case[:n]) --num-steps 1000000 -v 2 --prefix $(joinpath(workdir, prefix(case)))`;
+    output = read(command, String);
+    write(outfile, output); 
+  else
+    println("Case: $case has already been run.");
+  end
+
+end, cases);
diff --git a/run/interacting-dielectric-study_2021-10-27.jl b/run/interacting-dielectric-study_2021-10-27.jl
@@ -0,0 +1,50 @@
+println("Hello world!");
+
+@everywhere using Printf;
+
+@show workdir = if length(ARGS) > 0
+  ARGS[1];
+else
+  ".";
+end
+
+@everywhere fmt(x) = @sprintf("%07d", round(Int, 1e3*x));
+@everywhere fmt_int(x) = @sprintf("%03d", x);
+
+@everywhere function prefix(case)
+  "E0-$(fmt(case[:E0]))_K1-$(fmt(case[:K1]))_K2-$(fmt(case[:K2]))_kT-$(fmt(case[:kT]))_Fz-$(fmt(case[:Fz]))_Fx-$(fmt(case[:Fx]))_n-$(fmt(case[:n]))_b-$(fmt(case[:b]))";
+end
+
+cases = Any[];
+E0s = [0.0; 1.0; 2.0; 3.0; 5.0];
+K1s = [1.0]; 
+K2s = [0.0]; 
+kTs = [1.0];
+Fzs = vcat(0.0:0.05:1.0, 1.5:0.5:5.0);
+Fxs = [0.0];
+ns = Int[100];
+#bs = [0.5; 1.0; 2.0];
+bs = [1.0];
+for b in bs, n in ns, Fx in Fxs, Fz in Fzs, kT in kTs, E0 in E0s, K1 in K1s, K2 in K2s
+  push!(cases, Dict(:E0 => E0, :K1 => K1, :K2 => K2,
+                    :kT => kT, :Fz => Fz,
+                    :Fx => Fx, :n => n, :b => b, :run => run));
+end
+
+@info "total number of cases to run: $(length(cases))";
+
+mkpath(workdir);
+
+pmap(case -> begin;
+
+  outfile = joinpath(workdir, "$(prefix(case)).out");
+  if !isfile(outfile)
+    println("Running case: $case.");
+    command = `julia -O 3 mcmc_clustering_eap_chain.jl --chain-type dielectric --energy-type Ising -b $(case[:b]) --E0 $(case[:E0]) --K1 $(case[:K1]) --K2 $(case[:K2]) --kT $(case[:kT]) --Fz $(case[:Fz]) --Fx $(case[:Fx]) -n $(case[:n]) --num-steps 1000000 -v 2 --prefix $(joinpath(workdir, prefix(case)))`;
+    output = read(command, String);
+    write(outfile, output); 
+  else
+    println("Case: $case has already been run.");
+  end
+
+end, cases);
diff --git a/run/interacting-dielectric-study_2021-10-29.jl b/run/interacting-dielectric-study_2021-10-29.jl
@@ -0,0 +1,49 @@
+println("Hello world!");
+
+@everywhere using Printf;
+
+@show workdir = if length(ARGS) > 0
+  ARGS[1];
+else
+  ".";
+end
+
+@everywhere fmt(x) = @sprintf("%07d", round(Int, 1e3*x));
+@everywhere fmt_int(x) = @sprintf("%03d", x);
+
+@everywhere function prefix(case)
+  "E0-$(fmt(case[:E0]))_K1-$(fmt(case[:K1]))_K2-$(fmt(case[:K2]))_kT-$(fmt(case[:kT]))_Fz-$(fmt(case[:Fz]))_Fx-$(fmt(case[:Fx]))_n-$(fmt(case[:n]))_b-$(fmt(case[:b]))";
+end
+
+cases = Any[];
+E0s = [0.0; 1.0; 3.0; 5.0];
+K1s = [1.0]; 
+K2s = [0.0]; 
+kTs = vcat(0.05:0.05:1.0, 1.5:0.5:5.0, 10.0, 100.0);
+Fs = [(0.0, 0.0), (0.0, 1.0), (1.0, 0.0), (1.0, 1.0)];
+ns = Int[100];
+#bs = [0.5; 1.0; 2.0];
+bs = [1.0];
+for b in bs, n in ns, (Fx, Fz) in Fs, kT in kTs, E0 in E0s, K1 in K1s, K2 in K2s
+  push!(cases, Dict(:E0 => E0, :K1 => K1, :K2 => K2,
+                    :kT => kT, :Fz => Fz,
+                    :Fx => Fx, :n => n, :b => b, :run => run));
+end
+
+@info "total number of cases to run: $(length(cases))";
+
+mkpath(workdir);
+
+pmap(case -> begin;
+
+  outfile = joinpath(workdir, "$(prefix(case)).out");
+  if !isfile(outfile)
+    println("Running case: $case.");
+    command = `julia -O 3 mcmc_clustering_eap_chain.jl --chain-type dielectric --energy-type Ising -b $(case[:b]) --E0 $(case[:E0]) --K1 $(case[:K1]) --K2 $(case[:K2]) --kT $(case[:kT]) --Fz $(case[:Fz]) --Fx $(case[:Fx]) -n $(case[:n]) --num-steps 1000000 -v 2 --prefix $(joinpath(workdir, prefix(case)))`;
+    output = read(command, String);
+    write(outfile, output); 
+  else
+    println("Case: $case has already been run.");
+  end
+
+end, cases);