On Fri, Sep 20, 2013 at 2:45 PM, Paul T. Bauman <[email protected]> wrote:
> On Fri, Sep 20, 2013 at 1:34 PM, Manav Bhatia <[email protected]>wrote: > >> >> Have there been any changes to System::project_solution() with >> FEMFunctionBase over the last 2-3 months? >> > > Current master shows FEMFunctionBase was touched 5 months ago. > > >> I had a piece of code that was working fine with it, but now is >> throwing >> an error at line 2339 in sytem_projection.C: >> > > Is this with the current git master? 2339 in system_projection.C is not a > libmesh_error()... What's the error you're getting? > 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. The stacktrace is (my internet machine is different from my coding machine, so I am typing this out): 0: 0 libMesh::print_trace 1: 1 libMesh::MacroFunctions::report_error 2: 2 QBase::n_points 3: 3 libMesh::ProjectFEMSolution::operator() 4: 4 libMesh::Threads::parallel_for<>... 5: 5 libMesh::System::project_vector 6: 6 libMesh::System::project_solution 7: 7 call from application code main Please let me know if you would need more information. > const unsigned int n_qp = qsiderule.n_points(); >> >> It seems like the FEMContext object (context) initialized at line 1843 >> in >> system_projection.C sets up its element_qrule and side_qrule with zero >> points, but edge_qrule has 2 points. >> > > This makes me a little nervous because we pushed API updates to > FEMContext, but it shouldn't have affected any of that (unless I screwed > something up obviously). > > Again, version of libMesh and the error your getting? stacktrace? > > Thanks, > > 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/20/13. http://pubads.g.doubleclick.net/gampad/clk?id=58041151&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
