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
_______________________________________________
Libmesh-users mailing list
Libmesh-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-users
