Dear Wolfgang,
In the meanwhile that the VectorTools::interpolate_based_on_material_id
function is taking care of, I just did the same through the following code
to change my solution vector in the specific regions with designated
material_id, which may be used by those who have the same issue with this
function until it is fixed. Hope it helps!
///////////////////////////////////
FEValues<dim> fe_values (fe, qf_cell, update_values );
std::vector<types::global_dof_index> local_dof_indices(dofs_per_cell);
typename DoFHandler<dim>::active_cell_iterator
cell = dof_handler.begin_active(),
endc = dof_handler_c.end();
for (; cell!=endc; ++cell)
if (cell->is_locally_owned())
{
fe_values.reinit(cell);
if (cell->material_id()==1)
{
cell->get_dof_indices (local_dof_indices);
for (unsigned int i=0; i<fe.dofs_per_cell; ++i)
solution(local_dof_indices[i])=1;
}
}
////////////////////////////////////////////
Thanks and regards,
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
