On Thu, 14 Feb 2008, Georg Wenig wrote:

> using libmesh 0.6.2 compiled with support for second derivatives the 
> following program crashes:

> The error is an invalid read attempt in fe_map.C, see below:
>
> error: attempt to subscript container with out-of-bounds index 0, but 
> container only holds 0 elements.
>
>
> #5  0x0000002a95abc368 in FEBase::compute_single_point_map (this=0x5697b0, 
> [EMAIL PROTECTED], elem=0x549d70, p=0)
>    at src/fe/fe_map.C:64

Thank you for the report.  On first glance, it looks like the problem
here is that FE<>::init_shape_functions is what resizes the storage
space for second derivatives, but that function gets overridden in
FEXYZ, which was never updated with second derivative support.  I'll
fix it in SVN this weekend.
---
Roy

-------------------------------------------------------------------------
This SF.net email is sponsored by: Microsoft
Defy all challenges. Microsoft(R) Visual Studio 2008.
http://clk.atdmt.com/MRT/go/vse0120000070mrt/direct/01/
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users

Reply via email to