[Dune] [#467] UGGrid test fails after calling a method of BasicGeometry

Dune dune at dune-project.org
Thu Nov 20 09:33:08 CET 2008


THIS IS AN AUTOMATED MESSAGE, DO NOT REPLY.

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

User who did this - Oliver Sander (sander) 

Attached to Project - Dune
Summary - UGGrid test fails after calling a method of BasicGeometry
Task Type - Bug Report
Category - Grid
Status - New
Assigned To - 
Operating System - Unspecified / All
Severity - Low
Priority - Normal
Reported Version - SVN
Due in Version - Undecided
Due Date - Undecided
Details - So far UGGrid's intersection geometries were handwritten 'generic geometry'-type classes.  In order to replace them with the new GenericGeometries I added proper BasicGeometries as base classes.  These base classes are set up, but not actually used yet.  One exception is the method jacobianInverseTransposed(), which had not been implemented at all so far.  Here the corresponding method in the BasicGeometry base class is called now.

The grid test crashes after calling this method, though.  It doesn't seem to be a bug in jacobianInverseTransposed(); I rather suspect something strange in the BasicGeometry memory handling (of course I may be wrong).  Here's what makes me think so:

Add the lines

std::cout << "i: " << i << "  local: " << coord_[i] << "   base: " << Base::operator[](i) << std::endl;
assert((coord_[i]-Base::operator[](i)).infinity_norm() < 1e-6);

to the operator[] (at uggridgeometry.hh:380), compile, and call the grid test.  You will notice that the corner coordinates of intersectionGlobal seem to change without intervention from the outside.

This is as far as I got. Some help from the makers of BasicGeometry would be appreciated.  Apologies in advance if it turns out that I screwed it up myself.




More information can be found at the following URL:
http://www.dune-project.org/flyspray/index.php?do=details&task_id=467

You are receiving this message because you have requested it from the Flyspray bugtracking system.  If you did not expect this message or don't want to receive mails in future, you can change your notification settings at the URL shown above.



More information about the Dune mailing list