// velocity and pressure block via <code>block_component</code>.
const std::vector<types::global_dof_index> dofs_per_block =
DoFTools::count_dofs_per_fe_block(dof_handler, block_component);
- const unsigned int n_u = dofs_per_block[0];
- const unsigned int n_p = dofs_per_block[1];
+ const types::global_dof_index n_u = dofs_per_block[0];
+ const types::global_dof_index n_p = dofs_per_block[1];
std::cout << " Number of active cells: " << triangulation.n_active_cells()
<< std::endl
const std::vector<types::global_dof_index> stokes_dofs_per_block =
DoFTools::count_dofs_per_fe_block(stokes_dof_handler, stokes_sub_blocks);
- const unsigned int n_u = stokes_dofs_per_block[0],
- n_p = stokes_dofs_per_block[1],
- n_T = temperature_dof_handler.n_dofs();
- const std::vector<unsigned int> stokes_block_sizes = {n_u, n_p};
+ const types::global_dof_index n_u = stokes_dofs_per_block[0],
+ n_p = stokes_dofs_per_block[1],
+ n_T = temperature_dof_handler.n_dofs();
std::cout << "Number of active cells: " << triangulation.n_active_cells()
<< " (on " << triangulation.n_levels() << " levels)" << std::endl
{
stokes_matrix.clear();
- BlockDynamicSparsityPattern dsp(stokes_block_sizes, stokes_block_sizes);
+ BlockDynamicSparsityPattern dsp(stokes_dofs_per_block,
+ stokes_dofs_per_block);
Table<2, DoFTools::Coupling> coupling(dim + 1, dim + 1);
Mp_preconditioner.reset();
stokes_preconditioner_matrix.clear();
- BlockDynamicSparsityPattern dsp(stokes_block_sizes, stokes_block_sizes);
+ BlockDynamicSparsityPattern dsp(stokes_dofs_per_block,
+ stokes_dofs_per_block);
Table<2, DoFTools::Coupling> coupling(dim + 1, dim + 1);
for (unsigned int c = 0; c < dim + 1; ++c)
const std::vector<types::global_dof_index> stokes_dofs_per_block =
DoFTools::count_dofs_per_fe_block(stokes_dof_handler, stokes_sub_blocks);
- const unsigned int n_u = stokes_dofs_per_block[0],
- n_p = stokes_dofs_per_block[1],
- n_T = temperature_dof_handler.n_dofs();
+ const types::global_dof_index n_u = stokes_dofs_per_block[0],
+ n_p = stokes_dofs_per_block[1],
+ n_T = temperature_dof_handler.n_dofs();
std::locale s = pcout.get_stream().getloc();
pcout.get_stream().imbue(std::locale(""));
const std::vector<types::global_dof_index> darcy_dofs_per_block =
DoFTools::count_dofs_per_fe_block(darcy_dof_handler,
darcy_block_component);
- const unsigned int n_u = darcy_dofs_per_block[0],
- n_p = darcy_dofs_per_block[1],
- n_s = saturation_dof_handler.n_dofs();
- const std::vector<unsigned int> darcy_block_sizes = {n_u, n_p};
+ const types::global_dof_index n_u = darcy_dofs_per_block[0],
+ n_p = darcy_dofs_per_block[1],
+ n_s = saturation_dof_handler.n_dofs();
std::cout << "Number of active cells: " << triangulation.n_active_cells()
<< " (on " << triangulation.n_levels() << " levels)" << std::endl
{
darcy_matrix.clear();
- BlockDynamicSparsityPattern dsp(darcy_block_sizes, darcy_block_sizes);
+ BlockDynamicSparsityPattern dsp(darcy_dofs_per_block,
+ darcy_dofs_per_block);
Table<2, DoFTools::Coupling> coupling(dim + 1, dim + 1);
for (unsigned int c = 0; c < dim + 1; ++c)
Mp_preconditioner.reset();
darcy_preconditioner_matrix.clear();
- BlockDynamicSparsityPattern dsp(darcy_block_sizes, darcy_block_sizes);
+ BlockDynamicSparsityPattern dsp(darcy_dofs_per_block,
+ darcy_dofs_per_block);
Table<2, DoFTools::Coupling> coupling(dim + 1, dim + 1);
for (unsigned int c = 0; c < dim + 1; ++c)