On Fri, 14 Jan 2011, yunfei zhu wrote: > I get warning like this: > WARNING: Second derivatives are not currently correctly calculated on > non-affine elements! > > I noted that this warning come from function FEBase::compute_map() in > fe_map.C file. > > 00497 #ifdef LIBMESH_ENABLE_SECOND_DERIVATIVES00498 if > (calculate_d2phi > <http://libmesh.sourceforge.net/doxygen/classlibMesh_1_1FEBase.php#a4170551fa0860822e21475e3c098bc13>)00499 > {00500 libmesh_do_once(00501 libMesh::err > <http://libmesh.sourceforge.net/doxygen/namespacelibMesh.php#ab0817b8542b150bfa3fb23393df99f28> > << "WARNING: Second derivatives are not currently "00502 > << "correctly calculated on non-affine elements!"00503 > << std::endl;);00504 }00505 #endif > > > I don't need to calculate the second derivatives for my problem and I did > not enable the second derivatives when I installed libmesh.
Yes, you did - they're enabled by default unless you add a configure option to explicitly disable them. If you're not using second derivatives, though, you don't need to worry about the warning. You can try prerequesting only the data you do need from the FE object - if you reinit() the object before requesting data then it'll precompute everything that you might request at a later time, including second derivatives, whether you need it all or not. > But why I get this warning. Because you had an element in your mesh that wasn't detected to be an "affine element" with linear mappings from master space to physical space. That actually shouldn't happen in a mesh that was created from a Tet4->Tet10 conversion unless you moved nodes afterward... if you want to hunt down the nodal coordinates of the offending element and send them to us it would be appreciated. --- Roy ------------------------------------------------------------------------------ Protect Your Site and Customers from Malware Attacks Learn about various malware tactics and how to avoid them. Understand malware threats, the impact they can have on your business, and how you can protect your company and customers by using code signing. http://p.sf.net/sfu/oracle-sfdevnl _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
