From 8f96079429f465cf51b21b7cca3d8898da38df34 Mon Sep 17 00:00:00 2001 From: Ericka Florio Date: Fri, 29 Nov 2024 11:23:47 +0000 Subject: [PATCH] Config for convergence off 256 with kstar~64 --- Examples/ScalarField/params.txt | 10 +++++----- .../ScalarFields/RandomField.impl.hpp | 2 +- 2 files changed, 6 insertions(+), 6 deletions(-) diff --git a/Examples/ScalarField/params.txt b/Examples/ScalarField/params.txt index 3352adf23..7fcdab05d 100644 --- a/Examples/ScalarField/params.txt +++ b/Examples/ScalarField/params.txt @@ -7,7 +7,7 @@ verbosity = 0 # location / naming of output files -output_path = "/home/eaf49/rds/hpc-work/dump/" # Main path for all files. Must exist! +output_path = "/home/eaf49/rds/hpc-work/paper1-data-sets/convergence-tests/kstar-32/N64" # Main path for all files. Must exist! chk_prefix = ScalarField_ plot_prefix = ScalarFieldp_ # restart_file = ScalarField_001000.3d.hdf5 @@ -16,9 +16,9 @@ plot_prefix = ScalarFieldp_ checkpoint_interval = 1000 # set to 0 to turn off plot files (except at t=0 and t=stop_time) # set to -1 to never ever print plotfiles -plot_interval = 0 -num_plot_vars = 6 -plot_vars = h11 h12 h13 h22 h23 h33 +plot_interval = 75 +num_plot_vars = 7 +plot_vars = h11 h12 h13 h22 h23 h33 Ham # subpaths - specific directories for hdf5, pout, extraction data # (these are created at runtime) @@ -44,7 +44,7 @@ E = 1.0 # Energy scale (in units of Mp!) scalar_amplitude = 4.0 scalar_velocity = -0.00001627328880468423 scalar_mass = 1.e-4 -# N_fine = 128 +N_fine = 256 tensor_amplitude = 1. # Choose a number between 0 and 9 diff --git a/Source/InitialConditions/ScalarFields/RandomField.impl.hpp b/Source/InitialConditions/ScalarFields/RandomField.impl.hpp index 7a5f4be3e..ac1f6e7b2 100644 --- a/Source/InitialConditions/ScalarFields/RandomField.impl.hpp +++ b/Source/InitialConditions/ScalarFields/RandomField.impl.hpp @@ -17,7 +17,7 @@ const double Mp = 1./m_bkgd_params.E; // Setting physical and window fn parameters - kstar = M_PI*((double) m_params.Nf)/m_params.L * 1./2.; + kstar = M_PI*((double) 64.)/m_params.L * 3./4.; epsilon = m_params.L/30.; H0 = sqrt((8.0 * M_PI/3.0/pow(Mp, 2.)) * (0.5*m_bkgd_params.velocity*m_bkgd_params.velocity