[Dune] Notification from Dune

Flyspray dune at hal.iwr.uni-heidelberg.de
Fri Feb 2 09:35:54 CET 2007


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

A new Flyspray task has been opened.  Details are below.

User who did this: - Uli Sack (usack)

Attached to Project - Dune
Summary - ownersFather of ALUvertices does not find father
Task Type - Bug Report
Category - Alu3D
Status - Unconfirmed
Assigned To - 
Operating System - All
Severity - Low
Priority - Normal
Reported Version - SVN
Due in Version - 
Due Date - Undecided
Details - ownersFather and consequently positionInOwnersFather do not
work properly with ALUSimplexGrid<2,2>. Invariably there is no father
found. The following code reproduces the problem:

#ifdef HAVE_CONFIG_H
# include "config.h"     
#endif

#include <fstream>
#include <iostream>

#include <dune/grid/alugrid.hh>
#include <dune/grid/io/file/dgfparser/dgfalu.hh>

using namespace Dune;

int main (int argc, char *argv[]) try
{
	// The grid
	const int dim = 2;
	typedef ALUSimplexGrid<dim,dim> GridType;
	Dune::GridPtr<GridType> gridptr;
	gridptr = Dune::GridPtr<GridType>("data/square.dgf");
	GridType& grid = *gridptr;
	grid.globalRefine(3);


	typedef GridType::Codim<dim>::LevelIterator NodeLevelIterator;
	typedef GridType::Codim<dim>::LeafIterator NodeLeafIterator;
	typedef GridType::Codim<0>::Entity::EntityPointer ElementPointer;

	FieldVector<double,dim> posInOwnersFather;
	
	//~ NodeLeafIterator nodeLeafEnd = grid.leafend<dim>();
	//~ for (NodeLeafIterator nodeLeafIt = grid.leafbegin<dim>();
nodeLeafIt!=nodeLeafEnd; ++nodeLeafIt)
	//~ {
		//~ ElementPointer ownersFather = (*nodeLeafIt).ownersFather();
		//~ posInOwnersFather = (*nodeLeafIt).positionInOwnersFather();
	//~ }
	
	NodeLevelIterator nodeLevelEnd = grid.lend<dim>(3);
	for (NodeLevelIterator nodeLevelIt = grid.lbegin<dim>(3);
nodeLevelIt!=nodeLevelEnd; ++nodeLevelIt)
	{
		ElementPointer ownersFather = (*nodeLevelIt).ownersFather();
		posInOwnersFather = (*nodeLevelIt).positionInOwnersFather();
	}
			
	return 0;
}

catch (Exception e)
{
	std::cout << e << std::endl;
}


More information can be found at the following URL:
http://hal.iwr.uni-heidelberg.de/flyspray/?do=details&id=252

You are receiving this message because you have requested it from the
Flyspray bugtracking system.  You can be removed from future
notifications by visiting the URL shown above.





More information about the Dune mailing list