Skip to content
Snippets Groups Projects
Commit b7ea0139 authored by Matthieu Schaller's avatar Matthieu Schaller
Browse files

Better parameters for the Zel'dovich pancake problem.

parent d0a7508f
No related branches found
No related tags found
1 merge request!566Periodic gravity calculation
...@@ -86,10 +86,8 @@ phi = sim["/PartType0/Potential"][:] ...@@ -86,10 +86,8 @@ phi = sim["/PartType0/Potential"][:]
x -= 0.5 * boxSize x -= 0.5 * boxSize
if snap > 40: # Check for Gadget solution
filename_g = "snapshot_%03d.hdf5"%(snap+1) filename_g = "snapshot_%03d.hdf5"%(snap)
else:
filename_g = "snapshot_%03d.hdf5"%(snap)
if os.path.exists(filename_g): if os.path.exists(filename_g):
sim_g = h5py.File(filename_g, "r") sim_g = h5py.File(filename_g, "r")
x_g = sim_g["/PartType0/Coordinates"][:,0] x_g = sim_g["/PartType0/Coordinates"][:,0]
......
...@@ -8,7 +8,10 @@ then ...@@ -8,7 +8,10 @@ then
fi fi
# Run SWIFT # 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 # Plot the result
python plotSolution.py 2 for i in {0..119}
do
python plotSolution.py $i
done
...@@ -33,7 +33,8 @@ InitialConditions: ...@@ -33,7 +33,8 @@ InitialConditions:
Scheduler: Scheduler:
max_top_level_cells: 8 max_top_level_cells: 8
cell_split_size: 400 cell_split_size: 50
tasks_per_cell: 125
Cosmology: Cosmology:
Omega_m: 1. Omega_m: 1.
...@@ -46,7 +47,7 @@ Cosmology: ...@@ -46,7 +47,7 @@ Cosmology:
Gravity: Gravity:
mesh_side_length: 16 mesh_side_length: 16
eta: 0.025 eta: 0.025
theta: 0.85 theta: 0.3
r_cut_max: 5. r_cut_max: 5.
comoving_softening: 0.0001 comoving_softening: 0.001
max_physical_softening: 0.0001 max_physical_softening: 0.001
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment