]> https://gitweb.dealii.org/ - dealii.git/commit
Fixing the compression bug of FETools::extrapolate.
authorGuilhem Poy <guilhem.poy@umontpellier.fr>
Fri, 25 Apr 2025 13:22:10 +0000 (15:22 +0200)
committerGuilhem Poy <guilhem.poy@umontpellier.fr>
Fri, 25 Apr 2025 13:22:10 +0000 (15:22 +0200)
commita9034913f2f7af7c350fe36735b54d34fefdfda9
treef85180593be364d74282787a8f02def6ee05d8af
parent1de1b85718a3a46f32af84817e6024315ad56d1d
Fixing the compression bug of FETools::extrapolate.

Due to the different compress(insert) behaviour of Petsc/Trilinos and
LinearAlgebra::distributed vectors, the previous version of the
extrapolate method was failing with the latter type of vector on
distributed simulations with locally owned DoF at the interface between
a locally owned coarse cell and a more refined neighbouring ghost cell.

This bug was not detected before by the fe_tools_extrapolate_*** tests
of the testsuite because the tested triangulation included no such DoFs.
A simple mirror inversion of the refinement pattern allows to detect the
bug without having to change the tested outputs (the extrapolated
function is invariant under the inversion symmmetry), and was also
included in this commit in addition to the fix.
include/deal.II/fe/fe_tools_extrapolate.templates.h
tests/mpi/fe_tools_extrapolate_common.h

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.