[Libmesh-users] augmented finite elements

[Tobias M. <tob...@un...>]

Dear all,

I hope that someone has experience or some good ideas to help me out: I 
am using a setup where the finite element basis is augmented by a 
function which contains the most rapidly changing features of the 
solution to allow for a considerably coarser grid.

However, when it comes to printing of the solution, I am stuck:
  - In best case, I'd like to use a refined mesh for printing. To 
represent the solution reasonably good, I have tried to copy the mesh 
into another one which I refine. My initial idea was to use a child 
class to MeshFunction to project the original solution with the 
augmentation function onto the finer grid; but since the 
"MeshFunction::operator()" is not virtual, one cannot overload it.
  - For the moment being, I am also happy to plot the solution on the 
coarser grid, hoping that the representation is not completely spoiled. 
However, I don't see a simple way to augment the solution-vector by a 
user-specified function?

I have tried to get this working manually, but at least in parallel mode 
always ended in a mess of functions and some contributions that were not 
available to the processor...

Many thanks in advance,
Tobias