[deal.II] Exception when using MGTransferMatrixFree::interpolate_to_mg with periodic BCs

2021-09-02 Thread Guilhem Poy
Dear deal-ii users and developers, I am currently implementing a code that minimizes an energy functional depending on a traceless symmetric tensor (5D solution) for liquid crystal physics. The physics details matters not, I will just mention that I am using a combination of trust-region Newto

[deal.II] Re: Exception when using MGTransferMatrixFree::interpolate_to_mg with periodic BCs

2021-09-03 Thread Guilhem Poy
Hi Peter, I just used your PR as a patch in my candi-deal-ii setup, and I am happy to report that indeed the bug was solved, thanks a lot for your reactivity! If you want to further simplify the new test in your PR, I think the AffineConstraints object is not needed. Now I just need to investi

[deal.II] Re: Exception when using MGTransferMatrixFree::interpolate_to_mg with periodic BCs

2021-09-03 Thread Guilhem Poy
ons). Or maybe the loss of performance is due to something else? Best regards, Guilhem On Friday, 3 September 2021 at 12:33:23 UTC+2 Guilhem Poy wrote: > Hi Peter, > > I just used your PR as a patch in my candi-deal-ii setup, and I am happy > to report that indeed the bug was solved, t

[deal.II] Re: Exception when using MGTransferMatrixFree::interpolate_to_mg with periodic BCs

2021-09-06 Thread Guilhem Poy
_constraints with the >> custom-built level_constraints in the level setup loop, but this does not >> work (and makes things worse in terms of iterations). Or maybe the loss of >> performance is due to something else? >> >> Best regards, >> >> Guilhem >>

[deal.II] Re: Exception when using MGTransferMatrixFree::interpolate_to_mg with periodic BCs

2021-09-06 Thread Guilhem Poy
rly visible when running the >>>> program, >>>> O(N) scaling is lost as soon as PBCs are activated, with terrible number >>>> of >>>> iterations (>100) for the solver. >>>> >>>> So probably the problem is that I do not

[deal.II] Non-square block with MatrixFreeOperators

2023-03-23 Thread Guilhem Poy
Hi everyone, I am trying to set-up a non-square MatrixFree operator for which the destination and source vectors do not have the same dimension (in order to use it in a Schur complement). The way I am doing this is by passing two dof_handlers to the MatrixFree object, and select row_block=0 an

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

2025-04-10 Thread Guilhem Poy
Dear deal-ii users and developers, I am trying to use the FETools::extrapolate method on a distributed triangulation that is adaptively refined. However, I get an exception inside this function every time my triangulation contains hanging nodes at ghost interfaces. I have attached a minimal wo

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

2025-04-20 Thread Guilhem POY
: https://www.math.colostate.edu/~bangerth/videos.676.32.8.html Thanks, and I’m sorry you got the impression no-one cared. It’s just that no-one knew, this time, the answer to your question… L. On 20 Apr 2025, at 09:54, Guilhem Poy wrote: I have designed a simple patch that fixes the bug for

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

2025-04-20 Thread Guilhem Poy
is is not the right channel to discuss all this? Best Guilhem On Wednesday, 16 April 2025 at 18:27:15 UTC+2 Guilhem Poy wrote: > 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

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

2025-04-16 Thread Guilhem Poy
thinking behind this dofs_on_refined_neighbors map? Best regards, Guilhem Le 15/04/2025 à 09:04, Guilhem Poy a écrit : 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

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

2025-04-25 Thread Guilhem Poy
Hi everyone, I had a bit of time today and as promised I have opened a bug report (https://github.com/dealii/dealii/issues/18389) and a PR (https://github.com/dealii/dealii/pull/18390). I wish you all a nice week-end! Best Guilhem On Monday, 21 April 2025 at 08:24:07 UTC+2 Guilhem Poy wrote

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

2025-04-10 Thread Guilhem Poy
mentationOutVector>::extrapolate_parallel. Best Guilhem On 10/04/2025 10:13, Guilhem Poy wrote: Dear deal-ii users and developers, I am trying to use the FETools::extrapolate method on a distributed triangulation that is adaptively refined. However, I get an exception inside this function

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

2025-04-15 Thread Guilhem Poy
t triggers the exception, while the one on the right don't: I will continue to investigate by diving into the internal functions of extrapolate, but I would welcome some help on this since I am not really familiar with the p4est code. Best regards, Guilhem On 10/04/2025 14:31, Guilh