There are a couple of points that I wanted to bring up (and actually meant
to originally... need to quit sending those emails at 3AM...) about this
patch.
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.
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
