Hello Sambit,

Can you  try doing this?
Also move the vertices of the ghost cell and avoid calling   
dftPtr->triangulation.communicate_locally_moved_vertices(locally_owned_vertices);

When I tried to use 
triangulation.communicate_locally_moved_vertices(locally_owned_vertices) 
last time, something wierd (Dont exactly remember what, Code was not doing 
what I was expecting) happened with me. But Since then, I have avoided 
using it.



Thanks.

On Monday, December 4, 2017 at 1:22:37 PM UTC-5, Sambit Das wrote:
>
> Hello Dr. Arndt,
>
> Thank you for your reply. My apologies for not being clear on the " breaks 
> the periodic face pairs match". Following is the error message I get when I 
> run on parallel.
>  
> *An error occurred in line <3699> of file 
> </home/vikramg/DFT-FE-softwares/softwareCentos/dealii/dealii-8.5.1/source/grid/grid_tools.cc>
>  
> in function*
> *    void 
> dealii::GridTools::match_periodic_face_pairs(std::set<std::pair<CellIterator, 
> unsigned int>, std::less<std::pair<CellIterator, unsigned int>>, 
> std::allocator<std::pair<CellIterator, unsigned int>>> &, 
> std::set<std::pair<dealii::identity<CellIterator>::type, unsigned int>, 
> std::less<std::pair<dealii::identity<CellIterator>::type, unsigned int>>, 
> std::allocator<std::pair<dealii::identity<CellIterator>::type, unsigned 
> int>>> &, int, 
> std::vector<dealii::GridTools::PeriodicFacePair<CellIterator>, 
> std::allocator<dealii::GridTools::PeriodicFacePair<CellIterator>>> &, const 
> dealii::Tensor<1, CellIterator::AccessorType::space_dimension, double> &, 
> const dealii::FullMatrix<double> &) [with CellIterator = 
> dealii::TriaIterator<dealii::CellAccessor<3, 3>>]*
> *The violated condition was: *
> *    n_matches == pairs1.size() && pairs2.size() == 0*
> *Additional information: *
> *    Unmatched faces on periodic boundaries*
>
>
> I suspect this has something to do with the ghost nodes across the 
> periodic boundary not being handled correctly. I am right now creating a 
> minimal working example of my bug. I will post that soon.
>
> Best,
> Sambit
>
>
> On Monday, December 4, 2017 at 8:53:51 AM UTC-5, Daniel Arndt wrote:
>>
>> Sambit,
>>
>> I am trying to move all parallel triangulation nodes by a constant 
>>> displacement, but that breaks the periodic face pairs match when I 
>>> call GridTools::collect_periodic_faces(...). I use the following code for 
>>> the mesh movement. The dftPtr->triangulation has periodicity constraints 
>>> using add_periodicity(...). 
>>>
>> What exactly do you mean by "that breaks the periodic face pairs match"? 
>> What is the error you are observing?
>> Can you provide us with a minimal example that shows the problem so we 
>> can check?
>>
>> Best,
>> Daniel
>>
>

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

Reply via email to