I have no issues. Also, storing the Jacobian is preferable to wasted flops if you need it too. 

-Ben



On Jul 25, 2012, at 6:24 PM, "Paul T. Bauman" <ptbauman@gmail.com> wrote:

On Wed, Jul 18, 2012 at 2:34 PM, Paul T. Bauman <ptbauman@gmail.com> wrote:
1. What does everyone think about deprecating the get_dxidx() etc. functions in favor of getting the FEMap object and getting those data from there (const accessor to the FEMap object is contained in this patch)? Not really a big deal, just struck me as conceptually more correct; you may disagree, of course.

2. Currently, only JxW are stored in FEMap. However, I'd like to also store, and provide access to, the determinant of the Jacobian of the map (i.e. J, without xW). This quantity is needed in the curl calculation for HCurl conforming elements and in the shape function evaluations in HDiv and L2 conforming elements (among others). Since typically few FE objects are instantiated (and therefore FEMap objects), this didn't strike me as a particularly heavy memory overhead, but would be a big FLOP saver. I haven't added this yet, but would like to once the final version of this patch goes in.

Are there any thoughts on the points above or on the patch itself? I would like to commit this patch (updated as needed for current trunk). Thanks,

Paul
------------------------------------------------------------------------------
Live Security Virtual Conference
Exclusive live event will cover all the ways today's security and
threat landscape has changed and how IT managers can respond. Discussions
will include endpoint security, mobile security and the latest in malware
threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel