Hello everyone.
I am solving two equations iteratively in an MPI framework. The first
equation is the equilibrium equation (slightly modified) solved in step-8
while the second equation (Damage equation) is very similar to the Laplace
equation solved in step-3.
The solution of equilibrium equation depends on the damage solution vector
and vice versa. To begin with, the damage solution vector is created as
follows:
(locally_relevant_solution_damage.reinit(locally_owned_dofs_damage,
locally_relevant_dofs_damage,
mpi_communicator);)
Next, for the first iteration, we need to set the damage solution as 1 for
certain nodes, as shown:
for (const auto &cell : dof_handler_damage.active_cell_iterators())
{
if (cell->is_locally_owned())
{
for (const auto vertex_number : cell->vertex_indices())
{
const auto vert = cell->vertex(vertex_number);
int a = cell->vertex_dof_index(vertex_number, 0);
if (condition is satisfied by the vertex)
{
locally_relevant_solution_damage[a] = 1;
}
else
{
locally_relevant_solution_damage[a] = 0;
}
}
}
}
However, I am getting a segmentation fault when I try to extract any entry
of locally_relevant_solution_damage.
I first thought that I needed to use the compress operation, but then I
read that compressing does not apply to vectors with ghost elements. Now, I
am unable to figure out what could be causing this segmentation fault.
Error message:
Primary job terminated normally, but 1 process returned
a non-zero exit code. Per user-direction, the job has been aborted.
--------------------------------------------------------------------------
--------------------------------------------------------------------------
mpirun noticed that process rank 0 with PID 0 on node wasim-OptiPlex-5080
exited on signal 11 (Segmentation fault).
--
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/3b8125ee-119c-4727-a64d-236f2cecb75dn%40googlegroups.com.
