Hi All,

I'm having problem when applying periodicity on meshes generated with Gmsh 
or Abaqus. The problem arises when I run in parallel (in serial it works) 
getting the following error:

 void 
dealii::GridTools::match_periodic_face_pairs(std::set<std::pair<CellIterator, 
unsigned int> >&, std::set<std::pair<typename 
dealii::identity<RangeType>::type, unsigned int> >&, int, 
std::vector<dealii::GridTools::PeriodicFacePair<CellIterator> >&, const 
dealii::Tensor<1, typename FaceIterator::AccessorType:: space_dimension>&, 
const dealii::FullMatrix<double>&) [with CellIterator = 
dealii::TriaIterator<dealii::CellAccessor<3, 3> >; typename 
dealii::identity<RangeType>::type = 
dealii::TriaIterator<dealii::CellAccessor<3, 3> >; typename 
FaceIterator::AccessorType = dealii::CellAccessor<3, 3>]
The violated condition was:
    n_matches == pairs1.size() && pairs2.size() == 0
The name and call sequence of the exception was:
    ExcMessage ("Unmatched faces on periodic boundaries")
Additional Information:
Unmatched faces on periodic boundaries 

I suspect that match_periodic_face_pairs function does not have access to 
the whole faces in parallel. 

I am wondering if you have encountered this issue before, or any clue what 
is the reason for this issue.

Thanks in advance.

Hamed

-- 
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.

Reply via email to