Commit 8628541d authored by Cristian Lalescu's avatar Cristian Lalescu
Browse files

hide some debug messages

to simplify debugging of stokes drag code.
parent 23cb8afc
......@@ -799,8 +799,8 @@ int field<rnumber, be, fc>::write_filtered(
}
}
}
DEBUG_MSG("count[0] = %ld, offset[0] = %ld\n",
count[0], offset[0]);
//DEBUG_MSG("count[0] = %ld, offset[0] = %ld\n",
// count[0], offset[0]);
if (nz>=2)
{
assert(nz%2==0);
......@@ -812,8 +812,8 @@ int field<rnumber, be, fc>::write_filtered(
offset[1] = cz*nz/2;
memshape [1] = this->clayout->sizes[1];
memoffset[1] = cz*(this->clayout->sizes[1] - nz/2);
DEBUG_MSG("cz = %d, count[1] + offset[1] = %ld\n",
cz, count[1] + offset[1]);
//DEBUG_MSG("cz = %d, count[1] + offset[1] = %ld\n",
// cz, count[1] + offset[1]);
//now write data
mspace = H5Screate_simple(ndim(fc), memshape, NULL);
......@@ -830,7 +830,7 @@ int field<rnumber, be, fc>::write_filtered(
offset[1] = 0;
memshape [1] = this->clayout->sizes[1];
memoffset[1] = 0;
DEBUG_MSG("filtered write nz=1\n");
//DEBUG_MSG("filtered write nz=1\n");
//now write data
mspace = H5Screate_simple(ndim(fc), memshape, NULL);
......@@ -921,9 +921,9 @@ void field<rnumber, be, fc>::compute_rspace_stats(
// and possible component indices that are NOT averaged over
// HDF5 dataset has 2 extra indices for time and order of moment,
// and possible component indices
DEBUG_MSG("ndims = %d, ndim(fc) = %d, dsetname = %s\n",
ndims, ndim(fc),
dset_name.c_str());
//DEBUG_MSG("ndims = %d, ndim(fc) = %d, dsetname = %s\n",
// ndims, ndim(fc),
// dset_name.c_str());
assert(ndims == int(ndim(fc))-1);
assert(dims[1] == nmoments);
switch(ndims)
......@@ -2130,12 +2130,12 @@ int make_gaussian_random_field(
rseed*omp_get_max_threads()*output_field->clayout->nprocs +
output_field->clayout->myrank*omp_get_max_threads() +
thread_id);
DEBUG_MSG("in make_gaussian_random_field, thread_id = %d, current_seed = %d\n", thread_id, current_seed);
//DEBUG_MSG("in make_gaussian_random_field, thread_id = %d, current_seed = %d\n", thread_id, current_seed);
rgen[thread_id].seed(current_seed);
}
output_field->real_space_representation = false;
*output_field = 0.0;
DEBUG_MSG("slope: %g\n", slope);
//DEBUG_MSG("slope: %g\n", slope);
// inside loop use only thread-local random number generator
kk->CLOOP_K2([&](
ptrdiff_t cindex,
......
......@@ -47,7 +47,7 @@ int NSVE_Stokes_particles<rnumber>::initialize(void)
particles_inner_computer_2nd_order_Stokes<double, long long int> current_particles_inner_computer;
current_particles_inner_computer.set_drag_coefficient(0.1);
DEBUG_MSG_WAIT(MPI_COMM_WORLD, "before call to particles_system_builder\n");
//DEBUG_MSG_WAIT(MPI_COMM_WORLD, "before call to particles_system_builder\n");
this->ps = particles_system_builder_with_p2p(
this->fs->cvelocity, // (field object)
this->fs->kk, // (kspace object, contains dkx, dky, dkz)
......@@ -63,7 +63,7 @@ int NSVE_Stokes_particles<rnumber>::initialize(void)
std::move(p2p_ghost_collisions<double, long long int>()),
std::move(particles_inner_computer_2nd_order_Stokes<double, long long int>()),
this->tracers0_cutoff);
DEBUG_MSG_WAIT(MPI_COMM_WORLD, "after call to particles_system_builder\n");
//DEBUG_MSG_WAIT(MPI_COMM_WORLD, "after call to particles_system_builder\n");
this->particles_output_writer_mpi = new particles_output_hdf5<
long long int, double, 6>(
MPI_COMM_WORLD,
......
......@@ -461,7 +461,7 @@ template <class rnumber,
void vorticity_equation<rnumber, be>::omega_nonlin(
int src)
{
DEBUG_MSG("vorticity_equation::omega_nonlin(%d)\n", src);
//DEBUG_MSG("vorticity_equation::omega_nonlin(%d)\n", src);
assert(src >= 0 && src < 3);
this->compute_velocity(this->v[src]);
/* get fields from Fourier space to real space */
......@@ -515,7 +515,7 @@ template <class rnumber,
field_backend be>
void vorticity_equation<rnumber, be>::step(double dt)
{
DEBUG_MSG("vorticity_equation::step\n");
//DEBUG_MSG("vorticity_equation::step\n");
TIMEZONE("vorticity_equation::step");
*this->v[1] = 0.0;
this->omega_nonlin(0);
......
Markdown is supported
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