I'm trying to make a new quadrature rule... that is different in that the points it integrates at are different on every element in the mesh.
This would seem to be an easy task but I am wondering about FE::reinit(). Reading through it, it seems as if it is missing the case where the quadrature points might change... ie it assumes that the quadrature points will stay constant as long as the element type (and order) doesn't change and if the shape doesn't change it is even more aggressive about not recomputing things. It does have a call to shapes_need_reinit()... but that is so that shape functions can override this caching behavior. It seems like we need a similar call to the quadrature rule object to know if the quadrature point locations have changed. Would anyone be opposed to me adding a qrule->has_constant_positions() (that would return true for all of our current rules and false for the new one I'm making)... or something like it (I'm open to suggestions!) Thanks! Derek ------------------------------------------------------------------------------ Centralized Desktop Delivery: Dell and VMware Reference Architecture Simplifying enterprise desktop deployment and management using Dell EqualLogic storage and VMware View: A highly scalable, end-to-end client virtualization framework. Read more! http://p.sf.net/sfu/dell-eql-dev2dev _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
