Dear Roberto,
I think the problem is that you are using GridTools::transform (via
GridTools::rotate), in a distributed triangulation. This should not be
done, as noted
here:
https://www.dealii.org/8.5.0/doxygen/deal.II/namespaceGridTools.html#a212e99cf0d923cebfa04f1d23fa60b04
If this is indee
Hi Matthias,
Thanks! I'll pass along the advice to fetch and build the development
version and otherwise wait for the 18.04 release.
Best,
Steve
On Thursday, April 12, 2018 at 4:36:53 PM UTC-4, Matthias Maier wrote:
>
> Status update:
>
> - We missed the opportunity to do a rebuild for 17.10 [
I've found what was wrong in this piece of code. It turned out that
index.first < dim does not guarantee this dof is a velocity dof. Changing
it to* if (i < fe.base_element(0).dofs_per_cell && index.first < dim)*
solves the problem.
But there are two things that I can not understand:
1. using i
Dear Lucas,
thanks very much for replying.
I tried to solve the same problem without rotating the domain. I rotated
the problem.
Unfortunately the issue is still there.
I attach the new images.
I also attach the source code again because I realized that the previous
one had a misplaced bracke
Wolfgang,
Thank you for your reply. If I understand it correctly, in step 4 of the
ParallelDistributed::distribute_dofs() function in the dof_handler.cc file
you compute and flag the cells that have ghost neighbours by looking at
only the vertices by calling the function
compute_vertices_with_