From 8f801c67a588ae9fa2e24a86cbbd909438d3d110 Mon Sep 17 00:00:00 2001 From: Bert Vandenbroucke <bert.vandenbroucke@gmail.com> Date: Tue, 21 Nov 2017 09:40:05 +0000 Subject: [PATCH] Changed theta parameter for 3D Evrard collapse test. --- examples/EvrardCollapse_3D/evrard.yml | 6 +----- examples/EvrardCollapse_3D/makeIC.py | 2 +- 2 files changed, 2 insertions(+), 6 deletions(-) diff --git a/examples/EvrardCollapse_3D/evrard.yml b/examples/EvrardCollapse_3D/evrard.yml index 72543ae40c..db69623044 100644 --- a/examples/EvrardCollapse_3D/evrard.yml +++ b/examples/EvrardCollapse_3D/evrard.yml @@ -32,14 +32,10 @@ SPH: Gravity: eta: 0.025 # Constant dimensionless multiplier for time integration. epsilon: 0.001 # Softening length (in internal units). - theta: 0.1 + theta: 0.9 a_smooth: 1.25 # (Optional) Smoothing scale in top-level cell sizes to smooth the long-range forces over (this is the default value). r_cut: 4.5 # (Optional) Cut-off in number of top-level cells beyond which no FMM forces are computed (this is the default value). # Parameters related to the initial conditions InitialConditions: file_name: ./evrard.hdf5 # The file to read - -Scheduler: - tasks_per_cell: 400 - cell_split_size: 20 # or so diff --git a/examples/EvrardCollapse_3D/makeIC.py b/examples/EvrardCollapse_3D/makeIC.py index 15886f88ea..beb133cca8 100644 --- a/examples/EvrardCollapse_3D/makeIC.py +++ b/examples/EvrardCollapse_3D/makeIC.py @@ -45,7 +45,7 @@ pos[:,2] = r * cos_theta pos += array([R, R, R]) -h = ones(numPart) * R / numPart**(1. / 3.) +h = ones(numPart) * 2. * R / numPart**(1. / 3.) # Generate extra arrays v = zeros((numPart, 3)) -- GitLab