Commit daa0299a authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

tweaks parameters for kraichnan test

parent 024712e0
......@@ -452,11 +452,6 @@ if (BUILD_TESTING)
### compare DNS output to stored results
add_test(
NAME test_NSVEparticles
COMMAND turtle.test_NSVEparticles
WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
### compare DNS output to stored results for two threads per process
add_test(
NAME test_NSVEparticles_ntpp2
COMMAND turtle.test_NSVEparticles --ntpp 2
WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
### simple runs of post-processing tools
......@@ -495,7 +490,7 @@ if (BUILD_TESTING)
WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
add_test(
NAME test_kraichnan_field
COMMAND turtle DNS kraichnan_field --simname dns_kraichnan_field --nparticles 10000
COMMAND turtle DNS kraichnan_field --simname dns_kraichnan_field --dtfactor 0.05 --nparticles 10000
WORKING_DIRECTORY ${TEST_OUTPUT_DIRECTORY})
add_test(
NAME test_NSVE_Stokes_particles
......
......@@ -702,7 +702,7 @@ class DNS(_code):
pars['field_random_seed'] = int(1)
pars['spectrum_slope'] = float(-5./3)
pars['spectrum_k_cutoff'] = float(16)
pars['spectrum_coefficient'] = float(1)
pars['spectrum_coefficient'] = float(0.1)
if dns_type == 'NSVE_Stokes_particles':
pars['initial_field_amplitude'] = float(0.0)
pars['initial_particle_vel'] = float(0.05)
......
......@@ -139,12 +139,13 @@ kspace<be, dt>::kspace(
// single thread because we're filling up unstructured map
DEBUG_MSG("before computing dealias_exp_factor\n");
DEBUG_MSG("dkx = %g, dky = %g, dkz = %g\n", this->dkx, this->dky, this->dkz);
if (dt == SMOOTH)
{
double max_kk2 = pow(std::numeric_limits<double>::max(), 0.95/18);
for (hsize_t yindex = 0; yindex < this->layout->subsizes[0]; yindex++)
for (hsize_t zindex = 0; zindex < this->layout->subsizes[1]; zindex++)
for (hsize_t xindex = 0; xindex < this->layout->subsizes[1]; xindex++)
for (hsize_t xindex = 0; xindex < this->layout->subsizes[2]; xindex++)
{
double kk2 = (pow(this->kx[xindex]/this->kMx, 2) +
pow(this->ky[yindex]/this->kMy, 2) +
......@@ -160,7 +161,7 @@ kspace<be, dt>::kspace(
this->dealias_exp_factor[int64_t(round(k2))] = value;
}
}
DEBUG_MSG_WAIT(this->layout->comm, "before computing dealias_exp_factor\n");
DEBUG_MSG_WAIT(this->layout->comm, "after computing dealias_exp_factor\n");
MPI_Allreduce(
nshell_local_thread.getMasterData(),
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment