[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