OK, I've got it. PR forthcoming. Manav was on the right track moving the side_fe_reinit around. reinit needs to be get called to properly initialize the quadrature rule. However, the reason this didn't fail in 2D or on the edge in 3D was because QGauss (the default quadrature rule here) internally inits the rule if the dimension is 1D (sides in 2D, edges in 3D). Fix forthcoming.
Thanks again for reporting, Manav, and sorry for the trouble. Best, Paul On Sat, Sep 21, 2013 at 10:17 AM, Paul T. Bauman <[email protected]> wrote: > No problem! This is still good info - probably means this has been broken > since I originally wrote that code. > > > > On Sep 21, 2013, at 7:49 AM, Manav Bhatia <[email protected]> wrote: > > Hi Paul, > > Now that I think of it, I had never before tested this for 3d cases. > > Sorry that I made the wrong claim. > > Manav > > > > On Sep 21, 2013, at 2:23 AM, "Paul T. Bauman" <[email protected]> wrote: > > On Fri, Sep 20, 2013 at 10:26 PM, Paul T. Bauman <[email protected]>wrote: > >> >> Thanks for the test case! I'll play with this in a little while and see >> if I can reproduce. >> > > OK, I can reproduce. Will figure out the problem and get a PR in. You said > this worked previously - any chance you have a hash or version number > associated with that? > > Thanks for reporting. > > Best, > > Paul > > ------------------------------------------------------------------------------ LIMITED TIME SALE - Full Year of Microsoft Training For Just $49.99! 1,500+ hours of tutorials including VisualStudio 2012, Windows 8, SharePoint 2013, SQL 2012, MVC 4, more. BEST VALUE: New Multi-Library Power Pack includes Mobile, Cloud, Java, and UX Design. Lowest price ever! Ends 9/22/13. http://pubads.g.doubleclick.net/gampad/clk?id=64545871&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
