Hello everyone.
I am trying to apply a particular boundary condition at one single point in
the domain for a problem very similar to step-40.
I am using the following code to do it:
for (const auto &cell : dof_handler.active_cell_iterators())
if (cell->is_locally_owned())
{
for (const auto &face : cell->face_iterators())
for (const auto vertex_number : cell->vertex_indices())
{
const auto vert = cell->vertex(vertex_number);
const Point<2>& node = vert;
if (std::fabs(node(0)<1e-12) )
{
const unsigned int dof = cell->vertex_dof_index(vertex_number, 0);
constraints.add_line(dof);
constraints.set_inhomogeneity(dof,1);
} // This set a scalar at node to value 1
}
}
However, I am getting a segmentation fault. Kindly provide some help in
this regard.
Thanks
Wasim
--
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].
To view this discussion on the web visit
https://groups.google.com/d/msgid/dealii/eea6e885-5e7a-4330-b769-4272a70af626n%40googlegroups.com.
