Hi Konrad, I have the following suggestion. You can use (1) https://www.dealii.org/current/doxygen/deal.II/p4est__wrappers_8cc.html check *vertices_with_ghost_neighbors*. (2)*p4est_search_local *check link https://github.com/cburstedde/p4est/blob/master/src/p4est_search.h
Dear all, I finally had time to go back to my code. Quick problem summary: I am currently working on a semi-Lagrangian method for an advection-diffusion equation with a multiscale method that I designed for advection dominated processes. During runtime I must evaluate the (known) multiscale
Thank you, Bruno. :-) On Tuesday, May 12, 2020 at 2:50:35 PM UTC+2, Bruno Turcksin wrote: > > Konrad, > > There is nothing out of the box. However, deal.II uses p4est which can > use more that one layer of ghost cells. So you should take a look > there to see how hard it is to change the code.
Konrad, There is nothing out of the box. However, deal.II uses p4est which can use more that one layer of ghost cells. So you should take a look there to see how hard it is to change the code. If that's the way you want to go, I am sure that someone will be able to provide you with some guidance.
Thank you, Bruno. What is sometimes done for semi-Lagrangian methods is that one defines a halo region around the locally owned cells (like ghost cells) that contains information from previous time steps and then limits the time step so that one never leaves the halo region when asking for
Konrad, Unfortunately, you will need to do the communication yourself. You can only evaluate the solution on cells that locally owned or on ghost cells. Best, Bruno On Monday, May 11, 2020 at 6:26:47 AM UTC-4, Konrad Simon wrote: > > Dear all, > > I am currently working on a semi-Lagrangian