On Sun, 21 Dec 2014, Manav Bhatia wrote: > It seems like line 445 in patch_recovery_error_estimator.C, which currently > reads > > > for (unsigned int i=0; i<psi.size(); i++) > F(i) = JxW[qp]*u_h*psi[i]; > > should be > > for (unsigned int i=0; i<psi.size(); i++) > F(i) += JxW[qp]*u_h*psi[i];
It looks like you're correct. I guess we just never noticed because a grossly misintegrated error indicator is still a half-decent error indicator? This looks like my fault too; it seems to be in the refactoring I did while adding support for additional integration norms. I'll commit the fix now; if anyone sees something we're missing they can revert later. Please tell me that you discovered this bug while you were leisurely reading through the code and not while you were pulling your hair out trying to figure out broken results? --- Roy ------------------------------------------------------------------------------ Download BIRT iHub F-Type - The Free Enterprise-Grade BIRT Server from Actuate! Instantly Supercharge Your Business Reports and Dashboards with Interactivity, Sharing, Native Excel Exports, App Integration & more Get technology previously reserved for billion-dollar corporations, FREE http://pubads.g.doubleclick.net/gampad/clk?id=164703151&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
