On 11/09/2016 08:24 AM, Ashkan Dorostkar wrote:
void
dealii::TrilinosWrappers::VectorBase::add(dealii::TrilinosWrappers::VectorBase::size_type,
const size_type*, const double*)
The violated condition was:
!has_ghost_elements()
The name and call sequence of the exception was:
ExcGhostsPresent()
Additional Information:
You are trying an operation on a vector that is only allowed if the
vector has no ghost elements, but the vector you are operating on does
have ghost elements. Specifically, vectors with ghost elements are
read-only and cannot appear in operations that write into these vectors.
See the glossary entry on 'Ghosted vectors' for more information.
Stacktrace:
-----------
#0 1 libdeal_II.g.8.4.2.dylib 0x0000000110560188
_ZN6dealii16TrilinosWrappers10VectorBase3addEjPKjPKd + 236: 1
libdeal_II.g.8.4.2.dylib 0x0000000110560188
_ZN6dealii16TrilinosWrappers10VectorBase3addEjPKjPKd
#1 2 libdeal_II.g.8.4.2.dylib 0x00000001121830f0
_ZNK6dealii16ConstraintMatrix26distribute_local_to_globalINS_16TrilinosWrappers17BlockSparseMatrixENS2_3MPI11BlockVectorEEEvRKNS_10FullMatrixINT_10value_typeEEERKNS_6VectorINT0_10value_typeEEERKSt6vectorIjSaIjEERS7_RSD_bNS_8internal9bool2typeILb1EEE
+ 3124: 2 libdeal_II.g.8.4.2.dylib 0x00000001121830f0
_ZNK6*dealii16ConstraintMatrix26distribute_local_to_global*INS_16TrilinosWrappers17BlockSparseMatrixENS2_3MPI11BlockVectorEEEvRKNS_10FullMatrixINT_10value_typeEEERKNS_6VectorINT0_10value_typeEEERKSt6vectorIjSaIjEERS7_RSD_bNS_8internal9bool2typeILb1EEE
#2 3 gia 0x000000010cb1a548
_ZN10GeoPhysics3gia8assembleEv + 13366: 3 gia
0x000000010cb1a548 _ZN10GeoPhysics3gia8assembleEv
#3 4 gia 0x000000010cb1bcfd
_ZN10GeoPhysics3gia3runEv + 421: 4 gia
0x000000010cb1bcfd _ZN10GeoPhysics3gia3runEv
system_rhs is not touched anywhere between.
So if you output immediately before the call to
distribute_local_to_global what
system_rhs.has_ghost_elements()
says, what do you find?
Best
W.
--
------------------------------------------------------------------------
Wolfgang Bangerth email: [email protected]
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 [email protected].
For more options, visit https://groups.google.com/d/optout.
