From b7ea0139f0ead25cb942b544de9e1b47773fd10d Mon Sep 17 00:00:00 2001
From: Matthieu Schaller <matthieu.schaller@durham.ac.uk>
Date: Sat, 30 Jun 2018 22:49:50 +0100
Subject: [PATCH] Better parameters for the Zel'dovich pancake problem.
---
examples/ZeldovichPancake_3D/plotSolution.py | 6 ++----
examples/ZeldovichPancake_3D/run.sh | 7 +++++--
examples/ZeldovichPancake_3D/zeldovichPancake.yml | 9 +++++----
3 files changed, 12 insertions(+), 10 deletions(-)
diff --git a/examples/ZeldovichPancake_3D/plotSolution.py b/examples/ZeldovichPancake_3D/plotSolution.py
index dd1c6243e4..163e61eb9d 100644
--- a/examples/ZeldovichPancake_3D/plotSolution.py
+++ b/examples/ZeldovichPancake_3D/plotSolution.py
@@ -86,10 +86,8 @@ phi = sim["/PartType0/Potential"][:]
x -= 0.5 * boxSize
-if snap > 40:
- filename_g = "snapshot_%03d.hdf5"%(snap+1)
-else:
- filename_g = "snapshot_%03d.hdf5"%(snap)
+# Check for Gadget solution
+filename_g = "snapshot_%03d.hdf5"%(snap)
if os.path.exists(filename_g):
sim_g = h5py.File(filename_g, "r")
x_g = sim_g["/PartType0/Coordinates"][:,0]
diff --git a/examples/ZeldovichPancake_3D/run.sh b/examples/ZeldovichPancake_3D/run.sh
index 196cd1f41a..9b6b8166ac 100755
--- a/examples/ZeldovichPancake_3D/run.sh
+++ b/examples/ZeldovichPancake_3D/run.sh
@@ -8,7 +8,10 @@ then
fi
# Run SWIFT
-../swift -s -c -G -t 4 zeldovichPancake.yml 2>&1 | tee output.log
+../swift -a -s -c -G -t 8 zeldovichPancake.yml 2>&1 | tee output.log
# Plot the result
-python plotSolution.py 2
+for i in {0..119}
+do
+ python plotSolution.py $i
+done
diff --git a/examples/ZeldovichPancake_3D/zeldovichPancake.yml b/examples/ZeldovichPancake_3D/zeldovichPancake.yml
index d0bbe55222..03445aa666 100644
--- a/examples/ZeldovichPancake_3D/zeldovichPancake.yml
+++ b/examples/ZeldovichPancake_3D/zeldovichPancake.yml
@@ -33,7 +33,8 @@ InitialConditions:
Scheduler:
max_top_level_cells: 8
- cell_split_size: 400
+ cell_split_size: 50
+ tasks_per_cell: 125
Cosmology:
Omega_m: 1.
@@ -46,7 +47,7 @@ Cosmology:
Gravity:
mesh_side_length: 16
eta: 0.025
- theta: 0.85
+ theta: 0.3
r_cut_max: 5.
- comoving_softening: 0.0001
- max_physical_softening: 0.0001
+ comoving_softening: 0.001
+ max_physical_softening: 0.001
--
GitLab