Hamed,
An error occurred in line <1480> of file
</shared/hpc/deal.ii-candi/deal.II-v8.4.1/include/deal.II/base/index_set.h> in
function
dealii::IndexSet::size_type
dealii::IndexSet::index_within_set(dealii::IndexSet::size_type) const
The violated condition was:
is_element(n) == true
The name and call sequence of the exception was:
ExcIndexNotPresent (n)
Additional Information:
The global index 0 is not an element of this set.
Stacktrace:
-----------
#0 ./parallel_large-strain_PT: dealii::IndexSet::index_within_set(unsigned
int) const
#1 /shared/hpc/deal.ii-candi/deal.II-v8.4.1/lib/libdeal_II.g.so.8.4.1:
dealii::PETScWrappers::internal::VectorReference::operator double() const
#2 ./parallel_large-strain_PT: PhaseField::Solid<3>::output_results() const
#3 ./parallel_large-strain_PT: PhaseField::Solid<3>::run()
#4 ./parallel_large-strain_PT: main
--------------------------------------------------------
Following is my out_put object. problem is related to the MappingQEulerian
since when I uncomment it (the red lines) I don't get error.
template <int dim>
void Solid<dim>::output_results() const
{
> [...]
LA::MPI::Vector soln ;
soln.reinit(locally_owned_dofs,mpi_communicator);
for (unsigned int i = 0; i < soln.size(); ++i)
soln(i) = solution(i);
You are trying to read *all* elements of the vector 'solution' here. But in
parallel, you do not have access to all elements but only those that are
stored on the processor on which this code is run. You need to read up on the
documentation of programs such as step-40 to understand which elements of a
parallel vector are stored on each processor.
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: bange...@colostate.edu
www: http://www.math.colostate.edu/~bangerth/
--
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 dealii+unsubscr...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.
