Re: [Libmesh-users] Example 11

2011-10-09 Thread Roy Stogner
On Sun, 9 Oct 2011, Christopher Petrovitch wrote: > Ah, it seems I miss understood the function n_vars(). I believe this is > whats needed. > > if (elem->node(n) == side->node(ns) && elem->get_node(n)->n_comp(0,p_var)!=0) This looks more right, yes. --- Roy --

Re: [Libmesh-users] Example 11

2011-10-09 Thread Christopher Petrovitch
Ah, it seems I miss understood the function n_vars(). I believe this is whats needed. if (elem->node(n) == side->node(ns) && elem->get_node(n)->n_comp(0,p_var)!=0) Thanks! Chris On Oct 9, 2011, at 4:33 PM, Kirk, Benjamin (JSC-EG311) wrote: > Since a node is a DofObject you can ask it if it

Re: [Libmesh-users] Example 11

2011-10-09 Thread Christopher Petrovitch
Ok, that makes sense. So something like: if (elem->node(n) == side->node(ns) && side->get_node(ns)->n_vars(p_var)!=0) to set the pressure variable in example 11? Christopher Petrovitch cpetr...@purdue.edu On Oct 9, 2011, at 4:33 PM, Kirk, Benjamin (JSC-EG311) wrote: > Since a node is a Dof

Re: [Libmesh-users] Example 11

2011-10-09 Thread Kirk, Benjamin (JSC-EG311)
Since a node is a DofObject you can ask it if it has any components for a particular variable of a particular system. If the number of components is 0 you obviously can't set anything. -Ben On Oct 9, 2011, at 12:06 PM, "Christopher Petrovitch" wrote: > Hello, > > I was wondering how you wo

[Libmesh-users] Example 11

2011-10-09 Thread Christopher Petrovitch
Hello, I was wondering how you would go about setting pressure boundary conditions in example 11. I followed the method used in the example and the program segfaults, which I suppose it should since the loop is over a Quad9, and pressure is setup to be a first order approximation. So, setting