[Libmesh-devel] Renumbering nodes in copy_nodes_and_elements()

2012-05-14 Thread David Andrs
Hi! I bumped into this problem in UnstructuredMesh: at the end of copy_nodes_and elements(), we call this->prepare_for_use(false) and we may end up with different node numbering. Since this is called from constructors and copy-constructors, the newly created object may have different node numb

[Libmesh-devel] Invalid inital values for scalar variables of higher order

2012-03-28 Thread David Andrs
rrect, are we not using the system properly? Thanks, -- David Andrs 0001-Fixing-initial-condition-setup-for-scalar-variables-.patch Description: Binary data -- This SF email is sponsosred by: Try Windows Azure free f

Re: [Libmesh-devel] Side effect of DofMap::constrain_element_matrix()

2011-09-29 Thread David Andrs
Roy Stogner wrote on 09/28/2011 05:02:21 PM: > > On Wed, 28 Sep 2011, David Andrs wrote: > > > We ran into an issue with our framework here at INL when forming > full jacobian matrix as a preconditioner with AMR. What we hit > > in our problem was the different size of

[Libmesh-devel] Side effect of DofMap::constrain_element_matrix()

2011-09-28 Thread David Andrs
call that constrain_element_matrix and we use this modified array for something later on? Thanks, -- David Andrs /* $Id: ex4.C 2501 2007-11-20 02:33:29Z benkirk $ */ /* The Next Great Finite Element Library. */ /* Copyright (C) 2003 Benjamin S. Kirk */ /* This library is free software; you can

Re: [Libmesh-devel] Configuring with Trilinos

2011-09-28 Thread David Andrs
t) and then if that's > not found checking for Makefile.export.{aztecoo,nox} (which really are > lower case, for Trilinos 8 and 9 support). We just must have > forgotten to put the latter check in a conditional - that should be > easy to fix. >

Re: [Libmesh-devel] Improving Trilinos support

2011-09-28 Thread David Andrs
Roy Stogner wrote on 09/23/2011 01:53:58 PM: > > > On Tue, 20 Sep 2011, Roy Stogner wrote: > > > There's a bug here: you forgot to change PetscNonlinearSolver to > > NoxNonlinearSolver in a couple START_LOG lines. > > One more serious bug, too: there's a use of libmesh_cast_ref which > breaks

Re: [Libmesh-devel] Support for Trilinos 10.6.4

2011-09-15 Thread David Andrs
Roy Stogner wrote on 09/13/2011 09:54:32 AM: > > On Mon, 12 Sep 2011, David Andrs wrote: > > > Attached is a patch that can configure libMesh against 9.0.3 and 10.6.4. > > Works great for me. Unless someone else reports breakage I'd say > commit it. > >

Re: [Libmesh-devel] Support for Trilinos 10.6.4

2011-09-12 Thread David Andrs
Roy Stogner wrote on 09/08/2011 03:40:08 PM: > > On Thu, 8 Sep 2011, David Andrs wrote: > > > I'll look into it. I was able to get 9.0.3 from sandia webpage. > > Thanks! If you can get things to configure against 9.0.3, that ought > to be sufficient. I tried onl

Re: [Libmesh-devel] Support for Trilinos 10.6.4

2011-09-08 Thread David Andrs
Roy Stogner wrote on 09/08/2011 03:04:03 PM: > > > On Thu, 8 Sep 2011, David Andrs wrote: > > > > > Roy Stogner wrote on 09/08/2011 02:29:14 PM: > > > > > > > > This isn't working at all for me yet - did you forget to add the new > &

Re: [Libmesh-devel] Support for Trilinos 10.6.4

2011-09-08 Thread David Andrs
Roy Stogner wrote on 09/08/2011 02:29:14 PM: > > This isn't working at all for me yet - did you forget to add the new > m4/trilinos.m4 to the patch? Unfortunately yes. I must have overlooked it among all of my other untracked files. Attached is a new patch with this file.

[Libmesh-devel] Support for Trilinos 10.6.4

2011-09-08 Thread David Andrs
get the right configure script - this patch is not modifying it. If you have questions/comment, let me know here. Thanks, -- David Andrs 0001-Adding-support-for-Trilinos-10.6.4.patch Description: Binary data -- Doing

Re: [Libmesh-devel] ex23 error in svn head

2011-09-07 Thread David Andrs
Can you try the attached patch? It should apply cleanly on r4840. It works for me, but I want to be sure that it fixes your issue as well - I'm not really familiar with ex23. If it works, I'll check that in. Thanks, -- David Andrs David Knezevic wrote on 09/07/2011 09:47:31 AM: &

Re: [Libmesh-devel] ex23 error in svn head

2011-09-07 Thread David Andrs
of(T), xdr_translator()); This seems to be catched by the safe guards of gcc 4.6.1. (which I have on that linux machine), but not by 4.2.1 which is on my Mac. I'll make a patch to fix this problem. Thanks for reporting, -- David Andrs David Knezevic wrote on 09/06/2011 08:13:17

Re: [Libmesh-devel] Restarting from an XDA file with subdomain restricted variables

2011-08-30 Thread David Andrs
Roy Stogner wrote on 08/26/2011 12:39:20 PM: > > On Fri, 26 Aug 2011, Kirk, Benjamin (JSC-EG311) wrote: > > > I'm happy to commit this - any objections? > > David Andrs has write permissions and can commit the patch himself. Thanks for the review, both patches are

Re: [Libmesh-devel] Restarting from an XDA file with subdomain restricted variables

2011-08-23 Thread David Andrs
LIBMESH_VERSION macro (above) can be easily used for version comparisons. Note: I was able to achieve this with less lines of code then the peek-into-stream solution and this seems a lot cleaner to me. Thanks, -- David Roy Stogner wrote on 08/22/2011 12:49:20 PM: > > On Mon, 22 Aug 2011,

[Libmesh-devel] Restarting from an XDA file with subdomain restricted variables

2011-08-22 Thread David Andrs
could probably peek the stream for opening parenthesis to see if this info is available there or not and take a proper action based on that. If you have any comments, please let me know. Meanwhile, I'll work on a patch that will enable this. Thanks, -- David Andrs ---

Re: [Libmesh-devel] XDA restart fix

2011-08-01 Thread David Andrs
We didn't have a test for that before - one of our post-docs ran into this. But now, we do have such a test. -- David Roy Stogner wrote on 08/01/2011 09:26:03 AM: > > On Mon, 25 Jul 2011, David Andrs wrote: > > > We bumped into an issue with XDA restart here at INL. We

[Libmesh-devel] XDA restart fix

2011-07-25 Thread David Andrs
the loop in ES::read() went over all systems there and we already had more than that. Attached patch is fixing this issue by storing a list of systems that were found in the XDA file and than using it instead of the member variable for setting the system vectors. Thanks, -- David Andrs 0001

[Libmesh-devel] Restart from XDA with new variables added

2011-07-21 Thread David Andrs
the right thing. Should this assert go away or be modified to n_assigned_vals >= vector_length? Thanks, -- David Andrs -- 5 Ways to Improve & Secure Unified Communications Unified Communications promises

Re: [Libmesh-devel] Ghosted solution vector not closed in update()

2011-06-10 Thread David Andrs
Roy Stogner wrote on 06/10/2011 10:11:16 AM: > > On Fri, 10 Jun 2011, Derek Gaston wrote: > > > On Jun 10, 2011, at 9:35 AM, Roy Stogner wrote: > > > >> Could you identify where in your failure case the solution vector is > >> getting left open? > > > > This happens during the call to update()

[Libmesh-devel] Ghosted solution vector not closed in update()

2011-06-10 Thread David Andrs
Hi! So we bumped into another issue with the recent libMesh update. It seems that System::update() does not close() the solution vector any more, however this call is crucial if the solution vector is ghosted. Attached is the patch that fixes it. Thanks, -- David 0001-Ghosted-solution-has-t

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 rep

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

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 t

[Libmesh-devel] Bug in TransientSystem::re_update

2011-06-08 Thread David Andrs
first_local_dof = 0 and end_local_dof = 0, thus it fails on assert in PetscVector::localize on line 873, size is 0 and last_local_dof is max integer value. Attached is the patch that reverts one hunk from r4500 to fix this issue for us. Thanks, -- David Andrs 0001-Fixing-TransientSystem-re_update.patch

[Libmesh-devel] PETSc vector fixes

2011-03-02 Thread David Andrs
. Thanks, -- David Andrs petsc-fix.patch Description: Binary data -- Free Software Download: Index, Search & Analyze Logs and other IT data in Real-Time with Splunk. Collect, index and harness all the fast moving IT

Re: [Libmesh-devel] Improving periodic boundaries

2010-11-01 Thread David Andrs
Roy Stogner wrote on 10/27/2010 02:50:30 PM: > > On Wed, 27 Oct 2010, David Andrs wrote: > > > Roy Stogner wrote on 10/27/2010 01:19:07 PM: > > > > > What kind of translations are they doing? > > > > Currently they do two types: > > > &

[Libmesh-devel] Mutex lock-up when using periodic boundary conditions in debug mode

2010-10-25 Thread David Andrs
increment_constructor_count is empty. Could anyone verify this? Thanks, -- David Andrs -- Nokia and AT&T present the 2010 Calling All Innovators-North America contest Create new apps & games for the Nokia N8 for c

[Libmesh-devel] Improving periodic boundaries

2010-10-25 Thread David Andrs
ange PeriodicBoundaries container to store pointers to PeriodBoundary instead of storing the class itself. There is a known memory leak in this patch, which I will fix once the changes are approved Any comments and suggestions are appreciated. Thanks, -- David Andrs periodic-bc.patch Description: Bi

Re: [Libmesh-devel] Time for 0.6.5 release?

2010-09-30 Thread David Andrs
Roy Stogner wrote on 09/30/2010 11:27:03 AM: > > On Wed, 29 Sep 2010, David Andrs wrote: > > > Cody and I prepared a patch for compiler.m4 file which should > address the build issues with native Apple compilers. It should not interfere > > with other build sys

Re: [Libmesh-devel] Time for 0.6.5 release?

2010-09-29 Thread David Andrs
Cody and I prepared a patch for compiler.m4 file which should address the build issues with native Apple compilers. It should not interfere with other build systems (see the attachment). -- David Andrs Roy Stogner wrote on 09/28/2010 03:58:17 PM: > > On Tue, 28 Sep 2010, Cody Permann