Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread David Andrs
Roy Stogner wrote on 06/08/2011 01:28:03 PM: > > On Wed, 8 Jun 2011, David Andrs wrote: > > > This was actually a serial run with adaptivity and this bug poped out > > after first refinement when calling EquationSystems::reinit(). > > > > I'll try to replicate this behavior with pure libMesh ca

Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread Roy Stogner
On Wed, 8 Jun 2011, David Andrs wrote: > This was actually a serial run with adaptivity and this bug poped out > after first refinement when calling EquationSystems::reinit(). > > I'll try to replicate this behavior with pure libMesh calls and send it to > you, it might be worth to include it int

Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread David Andrs
Roy Stogner wrote on 06/08/2011 12:45:06 PM: > > On Wed, 8 Jun 2011, David Andrs wrote: > > > > Would it be sufficient to change that assert from > > > libmesh_assert(last_local_idx < this->size()); > > > to > > > libmesh_assert(last_local_idx+1 <= this->size()); > > > ? > > > > Yes, it fixes

Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread Roy Stogner
On Wed, 8 Jun 2011, David Andrs wrote: > > Would it be sufficient to change that assert from > > libmesh_assert(last_local_idx < this->size()); > > to > > libmesh_assert(last_local_idx+1 <= this->size()); > > ? > > Yes, it fixes the issue as well. Good to hear; that's in SVN now. Sorry about t

Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread David Andrs
Roy Stogner wrote on 06/08/2011 11:52:55 AM: > > On Wed, 8 Jun 2011, David Andrs wrote: > > > Some of our tests here at INL started to fail after we updated > libMesh. I tracked the problem down into TransientSystem::re_update > where r4500 > > commented out the if-statement on line 152. With

Re: [Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread Roy Stogner
On Wed, 8 Jun 2011, David Andrs wrote: > Some of our tests here at INL started to fail after we updated libMesh. I > tracked the problem down into TransientSystem::re_update where r4500 > commented out the if-statement on line 152. Without it the code execution > continues into localize of olde

[Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread David Andrs
Hi! Some of our tests here at INL started to fail after we updated libMesh. I tracked the problem down into TransientSystem::re_update where r4500 commented out the if-statement on line 152. Without it the code execution continues into localize of older_local_solution and we have first_local_d