Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-14 Thread Kirk, Benjamin (JSC-EG311)
On Dec 13, 2012, at 4:43 PM, Jens Lohne Eftang wrote: >> I'll second that request, but for a different reason: If there was an >> error in user code that could have been caught with the right library >> assertion, that might be nice to add. > It wasn't libmesh-related at all, just a stupid typo

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Jens Lohne Eftang
On 12/13/2012 05:39 PM, Roy Stogner wrote: > > On Thu, 13 Dec 2012, John Peterson wrote: > >> On Thu, Dec 13, 2012 at 3:23 PM, Jens Lohne Eftang >> wrote: >>> So this turned out to be a rather humiliating error on my part. Sorry >>> about the spam but thanks at least for pointing me in the right

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Roy Stogner
On Thu, 13 Dec 2012, John Peterson wrote: > On Thu, Dec 13, 2012 at 3:23 PM, Jens Lohne Eftang wrote: >> So this turned out to be a rather humiliating error on my part. Sorry >> about the spam but thanks at least for pointing me in the right direction! > > Care to share any details in case it mi

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread John Peterson
On Thu, Dec 13, 2012 at 3:23 PM, Jens Lohne Eftang wrote: > So this turned out to be a rather humiliating error on my part. Sorry > about the spam but thanks at least for pointing me in the right direction! Care to share any details in case it might be of use to others in future? -- John --

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Jens Lohne Eftang
So this turned out to be a rather humiliating error on my part. Sorry about the spam but thanks at least for pointing me in the right direction! Best, Jens On 12/13/2012 04:13 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Dec 13, 2012, at 1:43 PM, John Peterson wrote: > >> On Thu, Dec 13, 2012 at

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Kirk, Benjamin (JSC-EG311)
On Dec 13, 2012, at 1:43 PM, John Peterson wrote: > On Thu, Dec 13, 2012 at 12:40 PM, Roy Stogner > wrote: >> >> On Thu, 13 Dec 2012, Jens Lohne Eftang wrote: >> >>> 2) The output of EquationSystems::print_info() after initialization >>> (after adding the second system) is listed below. It lo

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread John Peterson
On Thu, Dec 13, 2012 at 12:40 PM, Roy Stogner wrote: > > On Thu, 13 Dec 2012, Jens Lohne Eftang wrote: > >> 2) The output of EquationSystems::print_info() after initialization >> (after adding the second system) is listed below. It looks sound to me. > > Wait one second: > >> EquationSystems >>

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Jens Lohne Eftang
On 12/13/2012 02:40 PM, Roy Stogner wrote: > > On Thu, 13 Dec 2012, Jens Lohne Eftang wrote: > >> 2) The output of EquationSystems::print_info() after initialization >> (after adding the second system) is listed below. It looks sound to me. > > Wait one second: > >> EquationSystems >> n_systems(

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Roy Stogner
On Thu, 13 Dec 2012, Jens Lohne Eftang wrote: > 2) The output of EquationSystems::print_info() after initialization > (after adding the second system) is listed below. It looks sound to me. Wait one second: > EquationSystems > n_systems()=2 >System #0, "io_system" ... > n_vectors()=2

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Jens Lohne Eftang
On 12/13/2012 12:11 PM, Kirk, Benjamin (JSC-EG311) wrote: > I created a unit test in the latest libMesh source tree, and for the buffers > you provide the code is behaving as expected. > https://github.com/libMesh/libmesh/pull/9 > > So the issue is not the new VariableGroups patch, which is consi

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-13 Thread Kirk, Benjamin (JSC-EG311)
On Dec 12, 2012, at 9:15 PM, Jens Lohne Eftang wrote: >> >> in your top-level directory and try that build - what happens? > I get the same behavior as with 6525. > > Thanks, > Jens Jens, I created a unit test in the latest libMesh source tree, and for the buffers you provide the code is beha

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 10:04 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Dec 12, 2012, at 6:35 PM, Jens Lohne Eftang wrote: > >>> and repeat your experiment, replying with the results again for the first >>> three nodes, also with the >>> node->debug_buffer() as before? >>> >>> I'd expect system 0, vars 0,

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
On Dec 12, 2012, at 6:35 PM, Jens Lohne Eftang wrote: >> >> and repeat your experiment, replying with the results again for the first >> three nodes, also with the >> node->debug_buffer() as before? >> >> I'd expect system 0, vars 0,1,2 to have ids 0, 96, and 192, respectively, >> from this b

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 06:17 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Dec 12, 2012, at 4:11 PM, Jens Lohne Eftang wrote: > >> With revision 6524 it's still broken. Yep, I'm pulling from svn. > With r6524, will you please uncomment the section > >// std::cout << "base_idx, var, vg, vig, ncg, comp=

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
One more thing - you may want to put node->debug_buffer(); libmesh_error(); So you can force an abort to track down the output - uncommenting those lines will create a lot of noise on the console! -Ben On Dec 12, 2012, at 5:17 PM, "Kirk, Benjamin (JSC-EG311)" wrote: > On Dec 12, 2012, a

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
On Dec 12, 2012, at 4:11 PM, Jens Lohne Eftang wrote: > With revision 6524 it's still broken. Yep, I'm pulling from svn. With r6524, will you please uncomment the section // std::cout << "base_idx, var, vg, vig, ncg, comp=" //<< base_idx << " " //

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 05:09 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Dec 12, 2012, at 3:52 PM, "Jens Lohne Eftang" wrote: > >> On 12/12/2012 04:31 PM, Kirk, Benjamin (JSC-EG311) wrote: >>> Sorry, the debug_buffer() is a new feature - from when I thought it was my >>> patch. ;-) >>> >>> Have you run thi

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
On Dec 12, 2012, at 3:52 PM, "Jens Lohne Eftang" wrote: > On 12/12/2012 04:31 PM, Kirk, Benjamin (JSC-EG311) wrote: >> Sorry, the debug_buffer() is a new feature - from when I thought it was my >> patch. ;-) >> >> Have you run this on debug mode? > Hmm, I must have messed something up when I ro

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 04:31 PM, Kirk, Benjamin (JSC-EG311) wrote: > Sorry, the debug_buffer() is a new feature - from when I thought it was my > patch. ;-) > > Have you run this on debug mode? Hmm, I must have messed something up when I rolled back to 6524 earlier. Now I did get debug_buffer() to work an

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
Sorry, the debug_buffer() is a new feature - from when I thought it was my patch. ;-) Have you run this on debug mode? Reading the system from disk initialized the equation systems, and then you add another system? Conceptually this should work, but... That's what the asserts are there for!

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 04:11 PM, Kirk, Benjamin (JSC-EG311) wrote: > On Dec 12, 2012, at 3:07 PM, Roy Stogner wrote: > >> On Wed, 12 Dec 2012, Jens Lohne Eftang wrote: >> >>> The system that is read from file has three variables. Now, when I >>> iterate of the nodes of mesh, and do >>> >>> index_0 = node->

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
On 12/12/2012 04:07 PM, Roy Stogner wrote: > On Wed, 12 Dec 2012, Jens Lohne Eftang wrote: > >> The system that is read from file has three variables. Now, when I >> iterate of the nodes of mesh, and do >> >> index_0 = node->dof_number(0, i, 0); >> index_1 = node->dof_number(0,0,0); >> >> I alwas g

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread John Peterson
On Wed, Dec 12, 2012 at 2:07 PM, Roy Stogner wrote: > > On Wed, 12 Dec 2012, Jens Lohne Eftang wrote: > >> The system that is read from file has three variables. Now, when I >> iterate of the nodes of mesh, and do >> >> index_0 = node->dof_number(0, i, 0); >> index_1 = node->dof_number(0,0,0); >>

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Kirk, Benjamin (JSC-EG311)
On Dec 12, 2012, at 3:07 PM, Roy Stogner wrote: > On Wed, 12 Dec 2012, Jens Lohne Eftang wrote: > >> The system that is read from file has three variables. Now, when I >> iterate of the nodes of mesh, and do >> >> index_0 = node->dof_number(0, i, 0); >> index_1 = node->dof_number(0,0,0); >> >

Re: [Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Roy Stogner
On Wed, 12 Dec 2012, Jens Lohne Eftang wrote: > The system that is read from file has three variables. Now, when I > iterate of the nodes of mesh, and do > > index_0 = node->dof_number(0, i, 0); > index_1 = node->dof_number(0,0,0); > > I alwas get index_0==index_1 also when i = 1,2. I assume you

[Libmesh-users] EquationsSystems with two systems; node iteration

2012-12-12 Thread Jens Lohne Eftang
Hi all, I have an EquationSystems which holds two systems: EquationSystems es(mesh); es.read(filename), READ, EquationSystems::READ_HEADER | EquationSystems::READ_DATA | EquationSystems::READ_ADDITIONAL_DATA); System & sys = es.add_