On 12/12/2012 04:37 PM, Andrs, David wrote: > On Wed, Dec 12, 2012 at 2:25 PM, Jens Lohne Eftang <[email protected] > <mailto:[email protected]>> wrote: > > On 12/12/2012 04:19 PM, Andrs, David wrote: > > On Wed, Dec 12, 2012 at 1:52 PM, Jens Lohne Eftang > <[email protected] <mailto:[email protected]>>wrote: > > > >> 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_system<System>("sys"); > >> sys.add_variable("u"); > >> es.init(); > >> > >> 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); > >> > > Maybe it is a typo here, but the first argument to dof_number is > system > > index: so they should be something like 0 and 1? The second one > is variable > > number, so it should be 0 (as long as you have just one variable > in each > > system)? > > Yep, that's just a typo in my email. Thanks though. > > > So, with 1 equation_system class with 2 systems where each system has > 1 variable, there are going to be two solution vectors. Now, if DOFs > are associated in a deterministic manner (and I believe they are), it > is very likely that the same node is gonna have the same global index > in both systems, because that number is a global within the solution > vector. > Not for the same variable number. System 0 has three variables (so i = 0,1,2 in my code), and System 1 has one variable. Nevertheless the global returned indices are the same whatever i is ...
Jens > However, I am not an libMesh expert, so I might be wrong... > -- > David > > Best, > Jens > > > > -- > > David > > > > > >> I alwas get index_0==index_1 also when i = 1,2. > >> > >> Any idea what I might be missing here? Bot n_vars() and > n_dofs() on the > >> two systems returns what they're supposed to ... > >> > >> Best, > >> Jens > >> > >> > >> > > ------------------------------------------------------------------------------ > >> LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial > >> Remotely access PCs and mobile devices and provide instant support > >> Improve your efficiency, and focus on delivering more value-add > services > >> Discover what IT Professionals Know. Rescue delivers > >> http://p.sf.net/sfu/logmein_12329d2d > >> _______________________________________________ > >> Libmesh-users mailing list > >> [email protected] > <mailto:[email protected]> > >> https://lists.sourceforge.net/lists/listinfo/libmesh-users > >> > > > > ------------------------------------------------------------------------------ > > LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial > > Remotely access PCs and mobile devices and provide instant support > > Improve your efficiency, and focus on delivering more value-add > services > > Discover what IT Professionals Know. Rescue delivers > > http://p.sf.net/sfu/logmein_12329d2d > > _______________________________________________ > > Libmesh-users mailing list > > [email protected] > <mailto:[email protected]> > > https://lists.sourceforge.net/lists/listinfo/libmesh-users > > > > ------------------------------------------------------------------------------ > LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial > Remotely access PCs and mobile devices and provide instant support > Improve your efficiency, and focus on delivering more value-add > services > Discover what IT Professionals Know. Rescue delivers > http://p.sf.net/sfu/logmein_12329d2d > _______________________________________________ > Libmesh-users mailing list > [email protected] > <mailto:[email protected]> > https://lists.sourceforge.net/lists/listinfo/libmesh-users > > > > > -- > David Andrs ------------------------------------------------------------------------------ LogMeIn Rescue: Anywhere, Anytime Remote support for IT. Free Trial Remotely access PCs and mobile devices and provide instant support Improve your efficiency, and focus on delivering more value-add services Discover what IT Professionals Know. Rescue delivers http://p.sf.net/sfu/logmein_12329d2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
