Re: [deal.II] Particles and field interpolation error

2020-07-20 Thread luca.heltai
Dear Franco, what type of triangulation are you using? If you look at the code that is run with `update_cached_numbers` https://github.com/dealii/dealii/blob/9e50f3ac04d088a9d2ffa74a81f08348e7f98a65/source/particles/particle_handler.cc#L173 you will see that the numbers are updated (just as y

Re: [deal.II] ABOUT READING GMSH

2020-07-20 Thread Daniel Arndt
Xiang, deal.II interprets physical tags as boundary ids when reading .msh files using GridIn::read_msh ( https://www.dealii.org/current/doxygen/deal.II/classGridIn.html#a83872db02e04f52ac52d578912f6da5e ). If your mesh consists of multiple parts in separate msh files you can create triangulation f

[deal.II] ABOUT READING GMSH

2020-07-20 Thread xs...@berkeley.edu
Hi, can Dealii read the inside boundary from .gmsh file, which can distinguish different parts of model geometry? If it can not, how can dealii import the model with multiple parts? Thanks a lot. Best, Xiang -- The deal.II project is located at http://www.dealii.org/ For mailing list/forum op

Re: [deal.II] Particles and field interpolation error

2020-07-20 Thread franco.m...@gmail.com
Dear Luca, thanks for your answer, I've moved along but I am still encountering problems in the interpolation. This seems to be a problem with the way I use PETSc. So, as far as I understand, I need to initialize a PETSc vector of length n_global_particles * components (in my case, 2, the spac