Just found the the call to side_fe_reinit() in this for loop. It seems like the entire block from lines 2238 to 2247 should be moved to after the context.side_fe_reinit() ?
-Manav On Fri, Sep 20, 2013 at 4:21 PM, Manav Bhatia <bhatiama...@gmail.com> wrote: > I was comparing the code with AssemblyContributions in fem_system.C, and > see that context.side_fe_reinit() has not been called before the call to > get qsiderule.n_points(). Should that be added to the for loop on element > sides at line 2248 and move the get_side_qrule method call (line > 2238) after this call? > > -Manav > > On Fri, Sep 20, 2013 at 4:09 PM, Manav Bhatia <bhatiama...@gmail.com>wrote: > >> Yes, it is 2239. I am sorry about the typo and the resulting confusion. >> -Manav >> >> >> On Fri, Sep 20, 2013 at 4:07 PM, Roy Stogner <royst...@ices.utexas.edu>wrote: >> >>> >>> On Fri, 20 Sep 2013, Manav Bhatia wrote: >>> >>> Sorry I was not too specific. I am on the git master, and line 2339 is a >>>> call to qsiderule.n_points(). It has an assertion of >>>> libmesh_assert(!_points.empyt(**)), where the code throws an exception. >>>> >>> >>> I think Paul's (and my) confusion was a simple typo - you're referring >>> to line 2239, not 2339. >>> --- >>> Roy >>> >> >> > ------------------------------------------------------------------------------ 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/20/13. http://pubads.g.doubleclick.net/gampad/clk?id=58041151&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
