Skip to content
Snippets Groups Projects
Commit 5487807d authored by Martin Nolte's avatar Martin Nolte
Browse files

use ReferenceElement instead of GenericReferenceElement

[[Imported from SVN: r174]]
parent e682678b
No related branches found
No related tags found
No related merge requests found
......@@ -55,11 +55,11 @@ namespace Dune
GeometryType type = geometry.type();
const GenericReferenceElement< ctype, mydim > &refElement = GenericReferenceElements< ctype, mydim >::general(type);
const ReferenceElement< ctype, mydim > &refElement = ReferenceElements< ctype, mydim >::general(type);
// Test whether the return value of the method 'center' corresponds to the center of the
// reference element. That is the current definition of the method.
const FieldVector<ctype, coorddim> center = geometry.global(refElement.position(0,0));
const FieldVector< ctype, coorddim > center = geometry.global( refElement.position( 0, 0 ) );
if( std::abs( (geometry.center() - center).two_norm() ) > 1e-8 )
DUNE_THROW(Exception, "center() is not consistent with global(refElem.position(0,0)).");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment