// TODO: fix this (do not run compress here, but fail)
if (last_action == Insert)
{
+# ifdef DEBUG
int ierr;
- ierr = matrix->GlobalAssemble(*column_space_map,
- matrix->RowMap(),
- false);
+ ierr =
+# endif
+ matrix->GlobalAssemble(*column_space_map,
+ matrix->RowMap(),
+ false);
Assert(ierr == 0, ExcTrilinosError(ierr));
- (void)ierr; // removes -Wunused-but-set-variable in optimized mode
}
last_action = Add;
irregular_cells.back() = task_info.n_ghost_cells % n_lanes;
}
+#ifdef DEBUG
{
unsigned int n_cells = 0;
for (unsigned int i = 0; i < task_info.cell_partition_data.back();
n_cells += irregular_cells[i] > 0 ? irregular_cells[i] : n_lanes;
AssertDimension(n_cells, task_info.n_ghost_cells);
}
+#endif
task_info.cell_partition_data.push_back(
task_info.cell_partition_data.back() + n_ghost_slots);
}
unsigned int j_indices[dim];
internal::FE_Q_Base::zero_indices<dim>(j_indices);
+#ifdef DEBUG
double sum_check = 0;
+#endif
for (unsigned int j = 0; j < q_dofs_per_cell; j += dofs1d)
{
double val_extra_dim = 1.;
my_restriction(mother_dof, child_dof) = 1.;
else if (std::fabs(val) > eps)
my_restriction(mother_dof, child_dof) = val;
+#ifdef DEBUG
sum_check += val;
+#endif
}
internal::FE_Q_Base::increment_indices<dim>(j_indices,
dofs1d);
std::reverse(step_sizes[i].begin(), step_sizes[i].end());
}
+# ifdef DEBUG
double x = 0;
for (unsigned int j = 0; j < step_sizes.at(i).size(); ++j)
x += step_sizes[i][j];
Utilities::int_to_string(i) +
" must be equal to the distance of the two given "
"points in this coordinate direction."));
+# endif
}