The commit that I am using is: 489c9d198d4ff0654f7d1f69f09404a2613b7954
I will work on a minimal example, but in the mean time here are the
functions I use to create the sparsity pattern and assemble the system. The
code works perfectly if either I turn off mesh refinement but keep periodic
Andrew,
Still hard to say what the problem is. Which commit are you using?
We have quite a number of tests that make sure that adaptive grid
refinement also works with periodic boundary conditions.
Can you provide a full minimal example that shows the problem you need your
workaround for?
How
I just say a slight problem and changed this line the declaration of
periodicityVector to:
std::vector::cell_iterator>
> periodicityVector;
But I'm getting the same behavior.
On Friday, February 28, 2020 at 8:08:19 AM UTC-5, Andrew Davis wrote:
>
> Yes, here is the code I used to tell the
Yes, here is the code I used to tell the triangulation about the periodic
boundaries:
// create a mesh on [0,1] with colorized boundaries
GridGenerator::hyper_cube(triangulation, 0.0, 1.0, true);
// collect the periodic boundary pairs
std::vector::cell_iterator> >
periodicityVector;
Andrew,
Did you tell the triangulation that it should take periodic boundaries into
account? My best bet (with the information given) would be that you didn't
call parallel::distributed::Triangulation::add_periodicity()
HI all,
I'm trying to implement a time-dependent convection equation using DG
elements with periodic boundary conditions on and adaptive mesh.
However, when I try to adapt the mesh using periodic using this code:
// estimate the error in each cell
Vector