[deal.II] Re: Exception when using FETools::extrapolate with ghost hanging nodes

2025-04-20 Thread Guilhem Poy
I have designed a simple patch that fixes the bug for simple situations that were previously failing with deal.II 9.6.1's FETools::extrapolate. I will continue the testing and propose a PR together with a bug report on github once I have time. Given the lack of answers here, I guess this is not

[deal.II] Re: Exception when using FETools::extrapolate with ghost hanging nodes

2025-04-16 Thread Guilhem Poy
Dear all, I think I have pinpointed the problem, which comes from l. 724-733 in deal.II/fe/fe_tools_extrapolate.templates.h: when a cell which has a more refined ghost neighbour owns a DoF that lives at the interface between this cell and this refined neighbour, the associated DoF value will

[deal.II] Re: Exception when using FETools::extrapolate with ghost hanging nodes

2025-04-15 Thread Guilhem Poy
Hi everyone, I have continued trying to fix the issue below on my own. First, I discovered that the extrapolate function is thoroughly tested inside the testsuite, so I went there to see how it is used to fix my MWE, but I still get the exception. This led me to conclude that some triangulati

[deal.II] Re: Exception when using FETools::extrapolate with ghost hanging nodes

2025-04-10 Thread Guilhem Poy
As a complement (I realized I should have put this in the first email), here is (one of) the exception that is thrown with the attached MWE from my first mail: /An error occurred in line <738> of file in function     void dealii::Utilities::MPI::Partitioner::import_from_ghosted_array_finish