Re: [Libmesh-users] negative Jacobian with system.point_value(...)

2011-08-08 Thread Roy Stogner
On Mon, 8 Aug 2011, John Peterson wrote: > On Mon, Aug 8, 2011 at 9:29 AM, robert wrote: > >> I found the problem. GMVIO changes the order of the tetrahedra >> vertices: > > There is no problem with GMVIO. You just can't combine a mesh written > in GMV format and a solution written in xda/r fo

Re: [Libmesh-users] negative Jacobian with system.point_value(...)

2011-08-08 Thread John Peterson
On Mon, Aug 8, 2011 at 9:29 AM, robert wrote: > Hello, > > I have seen that there has already been a topic like this one some years > ago. However, I couldn't find a solution for the problem. > >  At the end of a calculation I want to evaluate the results by using > system.point_value(int,Point).

Re: [Libmesh-users] negative Jacobian with system.point_value(...)

2011-08-08 Thread robert
I found the problem. GMVIO changes the order of the tetrahedra vertices: *.ele file: 192 4 1 1 102317 82 2 472415202 3 5128 9172 *.gmv file: cells 192 tet 4 10 17 23 8 tet 4 47 15 24 20 tet 4 51 9 28 17 In contras

[Libmesh-users] negative Jacobian with system.point_value(...)

2011-08-08 Thread robert
Hello, I have seen that there has already been a topic like this one some years ago. However, I couldn't find a solution for the problem. At the end of a calculation I want to evaluate the results by using system.point_value(int,Point). Everything works fine if I do it immediately. However, if