[Dune] [#1004] Make VTKFunctionPtr a const-pointer

[Jö Fahlke]

Am Tue, 24. Jan 2012, 15:23:58 +0100 schrieb Benedikt Oswald:
> although related to the VTK writer, perhaps my comment is somewhat off-topic.
> 
> I observe, when using a virtual resampling of the solution with the VTK subsampling functionality
> that I 'lose' certain parts of the solution, specifically when using subsampling level = 2.
> 
> I have attached an example, in much compressed form as jpg.
> 
> 
> 
> In the case of subsampling level = 1, there is no such problem.
> 
> It would be nice to have a look at this problem as well, in case you plan to modify
> the VTK writer functionality anyway.

Just to let everyone know the resolution: This was a numerical artifact.  The
cut-plane in paraview happened to be right on the border between grid
elements.  The SubsamplingVTKWriter uses nonconforming mode -- thus it can
happen that there is a gap between neighbouring grid elements with a width in
the order of the numerical precision (e.g. 1e-16).  In this case, moving the
cut-plane bei 1e-14 in either direction worked around the issue.

Bye,
Jö.

-- 
Jorrit (Jö) Fahlke, Interdisciplinary Center for Scientific Computing,
Heidelberg University, Im Neuenheimer Feld 368, D-69120 Heidelberg
Tel: +49 6221 54 8890 Fax: +49 6221 54 8884

This is the first age that's paid much attention to the future, which
is a little ironic since we may not have one.
-- Arthur C Clarke
-------------- next part --------------
A non-text attachment was scrubbed...
Name: signature.asc
Type: application/pgp-signature
Size: 828 bytes
Desc: Digital signature
URL: <https://lists.dune-project.org/pipermail/dune/attachments/20120201/8212f801/attachment.sig>