search

ElmerCSC / elmerfem

Official git repository of Elmer FEM software
http://www.elmerfem.org
Other
1.14k stars 310 forks source link

VtuOutputSolver.F90 change DGIndexes #442

Closed fgillet closed 1 week ago

fgillet commented 5 months ago

setting "Discontinuous Galerkin = Logical True" in the vtu output solver change the solution; see e.g. the elmerice test elmerice/Tests/DGsolver/density.sif comment/uncomment the vtu output file.

Tracking what's happening, VtuOutputSolver calls CheckAndCreateDGIndexes in MainUtils; and as I understand this routine overwrite the previously existing DGIndexes (created by CreateDGPerm in MainUtils during initialisation).

DGIndexes in CheckAndCreateDGIndexes should be updated only if not previously existing? I think this was the intetion with the boolean Failed; but it is never used in the second part of the routine.

raback commented 1 week ago

Thanx Fabien for the report. This was fixed already in May, I close the ticket now.