#1334 A question regarding the levelView of a grid.
Metadata
| Property | Value |
|---|---|
| Reported by | Adrian Ngo (adrian.ngo@iwr.uni-heidelberg.de) |
| Reported at | Aug 21, 2013 17:21 |
| Type | Unknown |
| Version | 2.2 |
| Operating System | Unspecified / All |
| Last edited by | Oliver Sander (oliver.sander@tu-dresden.de) |
| Last edited at | Aug 23, 2013 08:17 |
| Closed by | Oliver Sander (oliver.sander@tu-dresden.de) |
| Closed at | Aug 23, 2013 08:17 |
| Closed in version | Unknown |
| Resolution | Fixed |
| Comment |
Description
Dear all,
refining a 10x10 elements YASP grid globally, say twice, gives me level-views with
grid.levelView(0).size(0) = 100 grid.levelView(1).size(0) = 400 grid.levelView(2).size(0) = 1600
So far so good. This is what I would expect from a level grid view. Each level of the level view contains the complete gridview.
No so with a locally refined UG(2d) cube grid (or a 3D ALUGRID): Again, we start with a 10x10 grid:
grid.levelView(0).size(0) = 100
Step1: After marking 10 elements for refinement, I get this:
grid.maxLevel() = 1 grid.levelView( grid.maxLevel() ).size(0) = 40 grid.leafView().size(0) = 130
Step2: After marking another 14 elements for refinement, I get this:
grid.maxLevel() = 2 grid.levelView( grid.maxLevel() ).size(0) = 16 grid.leafView().size(0) = 172
I cannot see what exactly a level view is representing here. Surely not the full view of a grid. Is there a special switch to obtain the full view on each level?
Background: What I would like to do is to create a series of adaptively refined solutions u_h {h = h0,...,hMAX} and then compute the L2 error of (u_h - u_hMAX) for each refinement level. u_hMAX is kind of a reference solution. Therefore, I need to store the solution and the corresponding gridview on each refinement level. I was hoping that levelView() would be helpful.
Any suggestions?
Best regards, Adrian