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,