From: Martin Kronbichler Date: Thu, 11 May 2017 19:30:53 +0000 (+0200) Subject: Fix bug in matrix free thread graph creation. X-Git-Tag: v9.0.0-rc1~1605^2 X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=cb5c93d8c694f580533c41d83503fa99081235a1;p=dealii.git Fix bug in matrix free thread graph creation. --- diff --git a/include/deal.II/matrix_free/dof_info.templates.h b/include/deal.II/matrix_free/dof_info.templates.h index 0edceea0d5..d9568242db 100644 --- a/include/deal.II/matrix_free/dof_info.templates.h +++ b/include/deal.II/matrix_free/dof_info.templates.h @@ -839,13 +839,11 @@ no_constraint: std::min(((size_info.boundary_cells_end+task_info.block_size-1)/ task_info.block_size)*task_info.block_size, size_info.n_macro_cells); - start_up = start_nonboundary; size_info.boundary_cells_end = start_nonboundary; } else { start_nonboundary = size_info.n_macro_cells; - start_up = size_info.n_macro_cells; size_info.boundary_cells_start = 0; size_info.boundary_cells_end = size_info.n_macro_cells; } @@ -948,10 +946,11 @@ no_constraint: { // put all cells up to begin_inner_cells into first partition. if // the numbers do not add up exactly, assign an additional block - if (start_nonboundary>0 && start_up == start_nonboundary) + if (start_nonboundary>0) { unsigned int n_blocks = ((start_nonboundary+task_info.block_size-1) /task_info.block_size); + start_nonboundary = 0; for (unsigned int cell=0; cell