#208 SGrid works with dimworld > dim
Metadata
Property | Value |
---|---|
Reported by | Stefan Vater (vater@math.fu-berlin.de) |
Reported at | Nov 27, 2006 10:27 |
Type | Bug Report |
Version | Git (pre2.4) [autotools] |
Operating System | Linux |
Last edited by | Oliver Sander (oliver.sander@tu-dresden.de) |
Last edited at | Sep 4, 2007 07:50 |
Closed by | Oliver Sander (oliver.sander@tu-dresden.de) |
Closed at | Sep 4, 2007 07:50 |
Closed in version | Unknown |
Resolution | Implemented |
Comment | I added some additional fixes to SGrid and the test suite so SGrid now actually passes the test. The embedding is such that all dimworld-dim additional coordinates are always zero. I don't really see any direct use for that, but it is a nice test for the test. |
Description
This bug report might be related to #82 (closed).
In contrary to the grid-howto (Quote: "The SGrid class does only support the case where dimension and world dimension are equal.") it is possible to initialize an SGrid object with dimworld > dim.
For example, when I change line 12 in gettingstarted.cc to
typedef Dune::SGrid<3,4> GridType;
everything works and I get the output:
=> SGrid(dim=3,dimworld=4) level 0 codim[0]=18 codim[1]=75 codim[2]=104 codim[3]=48 leaf codim[0]=18 codim[1]=75 codim[2]=104 codim[3]=48 leaf dim=3 geomTypes=((cube, 3)[0]=18,(cube, 2)[1]=75,(cube, 1)[2]=104,(cube, 0)[3]=48)