[dune-fem] Bug in jacobianAll

[Martin Nolte]

Hi Marco,

this seems to be related to dune-grid. Can you reproduce this behavior with a 
small code snipped based on dune-grid only? If so, would please consider filing 
a bug report.

In order to remove the entity pointer in a compatible manner, some complicated 
changes have recently been made. Maybe one of them broke our way of obtaining 
the geometry.

Best,

Martin

On 03/03/2015 08:19 PM, Agnese, Marco wrote:
> 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.
> _______________________________________________
> dune-fem mailing list
> dune-fem at dune-project.org
> http://lists.dune-project.org/mailman/listinfo/dune-fem
>

-- 
Dr. Martin Nolte <nolte at mathematik.uni-freiburg.de>

Universität Freiburg                                   phone: +49-761-203-5630
Abteilung für angewandte Mathematik                    fax:   +49-761-203-5632
Hermann-Herder-Straße 10
79104 Freiburg, Germany