Re: [deal.II] mapping of global dof indices after solution transfer

2024-06-04 Thread Wolfgang Bangerth
On 6/4/24 13:11, Simon wrote: I am storing a std::vector of size N holding global DOF indices living at some nodes of my mesh where we have measured experimental values. After doing a SolutionTransfer motivated by AMR, I want to have the new global DOF indices at the same nodes, so the new vec

[deal.II] mapping of global dof indices after solution transfer

2024-06-04 Thread Simon
Dear all, I am storing a std::vector of size N holding global DOF indices living at some nodes of my mesh where we have measured experimental values. After doing a SolutionTransfer motivated by AMR, I want to have the new global DOF indices at the same nodes, so the new vector of indices has th

Re: [deal.II] Re: gmsh imported into dealii

2024-06-04 Thread Bruno Turcksin
Wang, Do you have a geo file that describes the geometry? Can you remove the "Physical Surface" lines in the geo file and regenerate the mesh. Best, Bruno Le mar. 4 juin 2024 à 10:13, 王艺康 a écrit : > Hello, > Thank you very much for your reply. That's what I thought at first, but I > was wron

[deal.II] Re: gmsh imported into dealii

2024-06-04 Thread 王艺康
Hello, Thank you very much for your reply. That's what I thought at first, but I was wrong, when I removed the corresponding physical group it was still the same error. I think gmsh automatically sets the coupling interface to solid when drawing the model body mesh, but I don't know how to chang

[deal.II] Re: gmsh imported into dealii

2024-06-04 Thread Bruno Turcksin
Hello, Deal.II associates the physical group of a face with a boundary id. You need to remove the physical group that are associated with interior faces. Best, Bruno On Tuesday, June 4, 2024 at 3:28:05 AM UTC-4 yikangw...@gmail.com wrote: > I drew the model in gmsh but importing to dealii is

[deal.II] Block preconditioner for multigrid

2024-06-04 Thread Chayapol Chaoveeraprasit
Dear all, I am trying to use a multigrid preconditioner to solve a PDE with two domains, in which the resulting matrix will be a block matrix. I am trying to replicate what is in step-16 tutorial program, but turns out it doesn't work since there are separate dof handlers for two domains, and I