On Tue, May 5, 2015 at 6:34 PM, John Peterson wrote:
>
>
> On Tue, May 5, 2015 at 4:08 PM, David Knezevic > wrote:
>
>> On Tue, May 5, 2015 at 5:48 PM, John Peterson
>> wrote:
>>
>>>
>>>
>>> On Tue, May 5, 2015 at 2:04 PM, David Knezevic <
>>> david.kneze...@akselos.com> wrote:
>>>
Al
On Tue, May 5, 2015 at 4:08 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 5:48 PM, John Peterson
> wrote:
>
>>
>>
>> On Tue, May 5, 2015 at 2:04 PM, David Knezevic <
>> david.kneze...@akselos.com> wrote:
>>
>>>
>>> Also, going back to my negative jacobian issue: I realized that this
>>> ele
On Tue, May 5, 2015 at 5:48 PM, John Peterson wrote:
>
>
> On Tue, May 5, 2015 at 2:04 PM, David Knezevic > wrote:
>
>>
>> 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 jaco
On Tue, May 5, 2015 at 2:04 PM, David Knezevic
wrote:
>
> 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 tha
On Tue, May 5, 2015 at 2:27 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 4:24 PM, Roy Stogner
> wrote:
>
>>
>> On Tue, 5 May 2015, John Peterson wrote:
>>
>> static_bool is a reasonable fix, but I'd also be OK with just
>>> removing the call to elem->volume() from Elem::print_info(). In
On Tue, May 5, 2015 at 4:24 PM, Roy Stogner
wrote:
>
> On Tue, 5 May 2015, John Peterson wrote:
>
> static_bool is a reasonable fix, but I'd also be OK with just
>> removing the call to elem->volume() from Elem::print_info(). In
>> theory, you can always estimate that by knowing the node
>> loc
On Tue, 5 May 2015, John Peterson wrote:
static_bool is a reasonable fix, but I'd also be OK with just
removing the call to elem->volume() from Elem::print_info(). In
theory, you can always estimate that by knowing the node
locations...
One more alternative: try/catch in elem->volume(). If
On Tue, May 5, 2015 at 2:04 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 3:50 PM, Roy Stogner
> 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 trigge
On Tue, May 5, 2015 at 1:38 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 3:07 PM, David Knezevic > wrote:
>
>> On Tue, May 5, 2015 at 2:30 PM, John Peterson
>> wrote:
>>
>>>
>>>
>>> On Tue, May 5, 2015 at 12:20 PM, David Knezevic <
>>> david.kneze...@akselos.com> wrote:
>>>
On Tue, M
On Tue, May 5, 2015 at 3:50 PM, Roy Stogner
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 for
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 i
On Tue, May 5, 2015 at 3:07 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 2:30 PM, John Peterson
> wrote:
>
>>
>>
>> On Tue, May 5, 2015 at 12:20 PM, David Knezevic <
>> david.kneze...@akselos.com> wrote:
>>
>>> On Tue, May 5, 2015 at 10:27 AM, John Peterson
>>> wrote:
>>>
O
On Tue, May 5, 2015 at 2:30 PM, John Peterson wrote:
>
>
> On Tue, May 5, 2015 at 12:20 PM, David Knezevic <
> david.kneze...@akselos.com> wrote:
>
>> On Tue, May 5, 2015 at 10:27 AM, John Peterson
>> wrote:
>>
>>>
>>>
>>> On Tue, May 5, 2015 at 7:10 AM, David Knezevic <
>>> david.kneze...@aksel
On Tue, May 5, 2015 at 12:20 PM, David Knezevic
wrote:
> On Tue, May 5, 2015 at 10:27 AM, John Peterson
> wrote:
>
>>
>>
>> On Tue, May 5, 2015 at 7:10 AM, David Knezevic <
>> david.kneze...@akselos.com> wrote:
>>
>>> [0] /home/dknez/software/libmesh-src/src/fe/fe_map.C, line 618, compiled
>>> M
On Tue, May 5, 2015 at 10:27 AM, John Peterson wrote:
>
>
> On Tue, May 5, 2015 at 7:10 AM, David Knezevic > wrote:
>
>> [0] /home/dknez/software/libmesh-src/src/fe/fe_map.C, line 618, compiled
>> May 5 2015 at 08:08:28
>> ERROR: negative Jacobian: -1.79627e-07 in element 7051
>>
>
> Is your fe
On Tue, May 5, 2015 at 7:10 AM, David Knezevic
wrote:
> [0] /home/dknez/software/libmesh-src/src/fe/fe_map.C, line 618, compiled
> May 5 2015 at 08:08:28
> ERROR: negative Jacobian: -1.79627e-07 in element 7051
>
Is your fe_map.C not the version in master? Line 618 doesn't have a
libmesh_error
On Tue, 5 May 2015, David Knezevic wrote:
> The mesh looks fine to me (it's 1.3MB so I can't attach it here, but I'd be
> happy to email it to anyone who is interested).
Don't care to see the whole mesh, but if you "catch throw" in the
debugger and print out the list of point locations in the f
On Tue, May 5, 2015 at 10:10 AM, John Peterson wrote:
>
>
> On Tue, May 5, 2015 at 7:10 AM, David Knezevic > wrote:
>
>> I'm getting a negative jacobian error for a mesh created in CUBIT. CUBIT
>> says that the smallest jacobian in the mesh is around 1.e-4, but when I
>> load the mesh into libMe
On Tue, May 5, 2015 at 7:10 AM, David Knezevic
wrote:
> I'm getting a negative jacobian error for a mesh created in CUBIT. CUBIT
> says that the smallest jacobian in the mesh is around 1.e-4, but when I
> load the mesh into libMesh and call fe->reinit(elem) on each element, I get
> the error belo
I'm getting a negative jacobian error for a mesh created in CUBIT. CUBIT
says that the smallest jacobian in the mesh is around 1.e-4, but when I
load the mesh into libMesh and call fe->reinit(elem) on each element, I get
the error below.
The mesh looks fine to me (it's 1.3MB so I can't attach it h
20 matches
Mail list logo
