On Tue, May 5, 2015 at 2:04 PM, David Knezevic <david.kneze...@akselos.com>
wrote:
> On Tue, May 5, 2015 at 3:50 PM, Roy Stogner <royst...@ices.utexas.edu>
> wrote:
>
>>
>> On Tue, 5 May 2015, David Knezevic wrote:
>>
>> I think the problem is that calling "elem->print_info()" then calls
>>> "elem->volume()", which reinits an FE object, which triggers
>>> negative jacobian error again, and then you keep looping forever.
>>>
>>
>> Ha! We should have seen that coming.
>>
>> The best fix is probably not to revert the print_info(), though, but
>> to use a static bool to make sure we only hit it once. "failing=true;
>> print_info(); libmesh_error();" if it's false; zero out stuff and
>> return garbage (on the theory that *some* info is better than none) if
>> it's true.
>>
>
>
> Sounds fine to me... I can do this later, unless someone else wants to
> take care of it first.
>
> Also, going back to my negative jacobian issue: I realized that this
> element is actually not affine. So it seems like John's comment could be
> right, i.e. that the jacobian is negative at at least one of the quadrature
> points, and that Cubit doesn't pick this up because it's not testing the
> jacobian at the same points.
>
When you export second-order elements, I think Cubit automatically snaps
the mid-edge and mid-face nodes to the geometry.
--
John
------------------------------------------------------------------------------
One dashboard for servers and applications across Physical-Virtual-Cloud
Widest out-of-the-box monitoring support with 50+ applications
Performance metrics, stats and reports that give you Actionable Insights
Deep dive visibility with transaction tracing using APM Insight.
http://ad.doubleclick.net/ddm/clk/290420510;117567292;y
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel
