[dune-fem] Bug in jacobianAll
Agnese, Marco
m.agnese13 at imperial.ac.uk
Tue Mar 3 20:19:05 CET 2015
I have further localize the bug.
In the method
template< class Point, class JacobianRangeArray >
void jacobianAll ( const Point &x, JacobianRangeArray &jacobians ) const
in the header file
/dune/fem/space/basisfunctionset/default.hh
there is the instruction
const GeometryType &geo = geometry();
which gives a segmentation fault.
The method geometry() simply returns entity_->geometry() and entity_!=nullptr.
Nevertheless it doesn't like entity_->geometry() .
I haven't understand yet what is going wrong. Any ideas?
Cheers,
Marco.
More information about the dune-fem
mailing list