On Wed, Dec 12, 2012 at 2:25 PM, Jens Lohne Eftang <[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] > >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. 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] > >> 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] > > 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] > 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
