Re: [Libmesh-devel] Preconditioned adjoint solves with PETSc

2015-01-09 Thread Roy Stogner
On Thu, 8 Jan 2015, Jed Brown wrote: > Roy Stogner writes: > >> By "neither seems to do much"? >> >> I mean that if I set either "-ksp_pc_side right" or "-ksp_norm_type >> unpreconditioned", then instead of ending up falsely "converged" at >> 10 KSP preconditioned resid norm 1.580464771307e-13

[Libmesh-devel] PointLocatorTree Logic Issues?

2015-01-09 Thread Derek Gaston
I'm investigating an issue with PointLocators... and so I was digging through the logic in TreeNode... and came across some weirdness. There appears to be some logic mismatches between TreeNode::find_element() and TreeNode::find_element_in_children() I can see what the _intention_ was - but I don

Re: [Libmesh-devel] PointLocatorTree Logic Issues?

2015-01-09 Thread Derek Gaston
I seriously feel like I'm going crazy here. I commented out that "if" statement in TreeNode::find_element() but it still isn't helping. What's happening is that I have a particular setup that is _unstable_... i.e. I can run the same code with the same inputs and every once in a while it fails. A

[Libmesh-devel] Anyone using libMesh with MPICH?

2015-01-09 Thread Roy Stogner
I'm seeing the most astonishing bug: In devel mode, in introduction_ex4 of all places, miscommunications are causing parallel runs to hang in devel mode. Running through gdb I can actually watch as one process sends "0" and another receives "1". I haven't yet ruled out a bug on our part (this i

Re: [Libmesh-devel] Anyone using libMesh with MPICH?

2015-01-09 Thread Paul T. Bauman
I've been using mpich 3.1.x for months. Mostly the 0.9.4 branch. I'll give master a swing this weekend. > On Jan 9, 2015, at 6:06 PM, Roy Stogner wrote: > > > I'm seeing the most astonishing bug: > > In devel mode, in introduction_ex4 of all places, miscommunications > are causing parallel