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.
Thanks again, Brent 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. > --- > Roy ------------------------------------------------------------------------------ This SF.net email is sponsored by: SourcForge Community SourceForge wants to tell your story. http://p.sf.net/sfu/sf-spreadtheword _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
