Thanks everyone. Sorry about the silly error.
I did run into one more problem before I was able to get it working. The
mesh generator I was using uses a different handedness of tetrahedra
than libmesh was expecting. The result was that I was getting an error
message saying that I had a negative Jacobian. It took me a little while
to figure-out what the cause was. So, when documentation is written
about the format, it might be nice to include that point.
Roy Stogner wrote:
> On Thu, 29 Jan 2009, Brent Kraczek wrote:
>> ERROR: Finite element LAGRANGE on geometric element TET
>> only supports FEInterface::max_order = 1, not fe_type.order = 2
> John seems to have figured out the problem here... although the error
> message should be telling you TET4, not TET, which may be a regression
> from some of my recent code; I'll go double-check that now.
>> WARNING: Finite element LAGRANGE on geometric element TET
>> could not be p refined past FEInterface::max_order = 1
>> ERROR: unsigned char too small to hold Elem::_p_level!
>> Recompile with Elem:_p_level set to something bigger!
>> What I have found through digging is that when it goes to
>> reinitialize the mesh (DofMap::reinit) it finds that fe_type.order ==
>> 2; elem->p_level()==255
> This is probably a weird side-effect of the same problem (the p
> refinement code tries to bump down the p level to what the element can
> support, but because there's been no p refinement the level just goes
> from 0 to -1...); I'll add a check to fix that. But it could be an
> I/O issue too; let us know if it persists.