[Dune] Alugrid leveliterator assignment

Carsten Graeser graeser at math.fu-berlin.de
Fri Dec 29 11:49:53 CET 2006


Hello Dune,
is there a reason why assignment of leveliterators for different levels
is allowed for codim=0 and not allowed for codim=2 in ALUSimplexGrid<2,2> ?

For example

	GridType::Codim<codim>::LevelIterator it = grid.lbegin<codim>(0);
	it = grid.lbegin<codim>(1);

works fine for codim=0. For codim=2 this gives

Dune::ALU2dGridLevelIterator <2, pitype, GridImp>&
Dune::ALU2dGridLevelIterator<2, pitype,GridImp>::operator=
(const Dune::ALU2dGridLevelIterator<2, pitype, GridImp>&)
[with Dune::PartitionIteratorType pitype =  All_Partition, GridImp = const Dune::ALU2dGrid<2, 2>]: 
Assertion `&marker_ == &org.marker_' failed.

in dune-grid/dune/grid/alugrid/2d/iterator_imp.cc:1151.
This assertion is not present in the codim=0 specialization.
Shouldn't such assigments by allowed since the level
is only a property of the object and not of the type ?

Carsten


----------------------------------------------------------------------
  Carsten Gräser

  Freie Universität Berlin
  FB Mathematik und Informatik (WE II)
  Arnimallee 2-6
  14195 Berlin, Germany

  office: +49-30 / 838-75349
  fax:    +49-30 / 838-54977
  email: graeser at math.fu-berlin.de




More information about the Dune mailing list