Issues related to usage of parallel::fullydistributed::triangulation

[Aditya Pratap Singh]

Hi,

I am working with dea...@9.2.0 for one of my research projects and one of the applications includes reading from a large mesh file. So I was using the parallel::distributed::triangulation in my code which was working perfectly but now as the size of my mesh files has increased more, I have to switch to parallel::fullydistributed::triangulation but getting errors in compilation of the code itself! I have written a basic example using p:f:t and tried to compile it (attached to the mail) but I am getting an undefined reference error (as shown below) although the code compiles well for the p:d:t case!

error: undefined reference to 'void dealii::GridTools::collect_periodic_faces<dealii::parallel::fullydistributed::Triangulation<2, 2> >(dealii::parallel::fullydistributed::Triangulation<2, 2> const&, unsigned int, unsigned int, int, std::vector<dealii::GridTools::PeriodicFacePair<dealii::parallel::fullydistributed::Triangulation<2, 2>::cell_iterator>, std::allocator<dealii::GridTools::PeriodicFacePair<dealii::parallel::fullydistributed::Triangulation<2, 2>::cell_iterator> > >&, dealii::Tensor<1, dealii::parallel::fullydistributed::Triangulation<2, 2>::space_dimension, double> const&, dealii::FullMatrix<double> const&)'


Can you please help in what is probably going wrong from my side in implementation of p:f:t?

-- 
Regards,
Aditya Pratap Singh
Master's Student (FAU Erlangen-Nurnberg)

[peterrum]

Just replace `dealii::parallel::fullydistributed::Triangulation<dim>` by `dealii::Triangulation<dim>` within `Check_meshPeriodicity()` . I think that should work.

https://www.dealii.org/developer/doxygen/deal.II/namespaceTriangulationDescription_1_1Utilities.html#af575881c2cf233fe6f85d1a3a65a73f6 and https://www.dealii.org/developer/doxygen/deal.II/namespaceTriangulationDescription_1_1Utilities.html#aefc3e841bcfd37714a07d04e42c8ffca show you how to construct p:f:T based on an external file.

Hope that helps!

PM

[Aditya Pratap Singh]

Hi,

This may help in solving the problem for this particular example but for my main code, as I said I want to work with p:f:t triangulation for efficiently working with large meshes. Is that we can't use p:f:t in member functions like collect_periodic_faces and have to use dealii::Triangulation only?  Also, as I mentioned earlier, if I use p:d:t triangulation otherwise, it doesn't throw any error, but p:f:t does which is quite a concern for me.

As I am a new user of dealii and still a learner, kindly correct me if I am interpreting the things wrongly.

Regards.

[Aditya Pratap Singh]

I have tried to make changes according to the link provided in past replies but still getting the undefined reference error ! I have attached the modified code herein for reference. If anyone can suggest something, it will be great help. 

[Praveen C]

In your check periodicity function you must replace p:f:t with Triangulation. 

Attached code has those changes and compiles.

best

praveen

[Aditya Pratap Singh]

Thanks for the reply.

Sadly my issue still remain unsolved.

This may help in compiling this particular code but lead to some other errors when I tried to run a test example for my original code. May be there are some other changes with respect to p:f:T that has to be made.

Can anyone share possible relevant source where I can find the information about how p:f:t behaves different to p:d:t / dealii::triangulation when given as an input to the member functions (get_periodic_face as an example)?