I realized my vector is not an MPI Vector and I am trying to use the
compress function. This won't work I assume. So I changed my global
variable to
LA::MPI::Vector configurational_forces;
This means I have to initialize it by the following command:
configurational_forces.reinit(locally_owned_dofs, mpi_com);
Now I have a problem. If I use the above approach and synchronize.
Everything works, but when I try to output or store, I think something goes
wrong. In my opinion, I just need to initialize the configurational_forces
vector with a larger IndexSet and not just to a locally owned one. Is that
possible?
When I run the problem I receive the following error:
----------------------------------------------------
Exception on processing:
--------------------------------------------------------
An error occurred in line <131> of file
</home/seyedali/programming/c++/projects/dealii/source/lac/petsc_vector_base.cc>
in function
dealii::PETScWrappers::internal::VectorReference::operator double()
const
The violated condition was:
(index >= static_cast<size_type>(begin)) && (index <
static_cast<size_type>(end))
Additional information:
You tried to access element 18 of a distributed vector, but only
elements 0 through 17 are stored locally and can be accessed.
--------------------------------------------------------
Aborting!
----------------------------------------------------
Any ideas?
Kind regards,
S. A. Mohseni
--
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.
