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)?

--
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

Reply via email to