On Jul 4, 2013, at 8:15 AM, David Knezevic <dkneze...@seas.harvard.edu> wrote:
>> The reason I ask is that we've identified what I believe to be a bug >> with Petsc block matrices where they are not honoring >> KEEP_NONZERO_PATTERN. This is causing a ton of mallocs for us in >> MOOSE. So, for now, we're just disabling block matrix support in >> libmesh until we get some clarity on the issue... > > OK, interesting. Disabling block matrix support (using Ben's hard-coding > suggestion) fixed the mallocs for me as well, so sounds like I might > have been running into the same issue. Do we have a consensus path forward here for 0.9.2? I am OK leaving the '-mat_new_nonzero_allocation_err false' option in the two RB examples that trigger the issue if everyone else is. As for Derek's concern that there is a bug in PETSc, if that is the case there is definitely a configure option to disable the blocked matrix support to work around that problem. I would like to get this out though, as the blocked matrix support provides pretty considerable speedup and memory savings. I'd also like to polish off 0.9.2 so we can get on with e.g. pthreads and other changes. -Ben ------------------------------------------------------------------------------ See everything from the browser to the database with AppDynamics Get end-to-end visibility with application monitoring from AppDynamics Isolate bottlenecks and diagnose root cause in seconds. Start your free trial of AppDynamics Pro today! http://pubads.g.doubleclick.net/gampad/clk?id=48808831&iu=/4140/ostg.clktrk _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
