Screenshot instructions:
Windows
Mac
Red Hat Linux
Ubuntu
Click URL instructions:
Rightclick on ad, choose "Copy Link", then paste here →
(This may not be possible with some types of ads)
From: Kirk, Benjamin (JSCEG311) <benjamin.kirk1@na...>  20130215 16:49:02

On Feb 15, 2013, at 10:41 AM, "David Knezevic" <dknezevic@...> wrote: > > On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: >> Hi, >> I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. > > Sure, that should be no problem. > > David The trick is the difference between geometric elements and finite elements. Your mesh will need Quad9 geometric elements, but you can add bilinear, biquadratic, hierarchic, whatever finite element approximations. Ben 
From: Subramanya Gautam Sadasiva <ssadasiv@pu...>  20130215 16:39:49

Hi, I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. Thanks, Subramanya Sadasiva "But memory eventually fades. Turbulences damp out, internal strains yield to plastic flow, concentration inhomogeneities diffuse to uniformity. Systems tend to subside to very simple states,independent of their specific history." Herbert Callen . Thermodynamics and an Introduction to Thermostatics. 
From: David Knezevic <dknezevic@se...>  20130215 16:41:31

On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: > Hi, > I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. Sure, that should be no problem. David 
From: Kirk, Benjamin (JSCEG311) <benjamin.kirk1@na...>  20130215 16:49:02

On Feb 15, 2013, at 10:41 AM, "David Knezevic" <dknezevic@...> wrote: > > On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: >> Hi, >> I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. > > Sure, that should be no problem. > > David The trick is the difference between geometric elements and finite elements. Your mesh will need Quad9 geometric elements, but you can add bilinear, biquadratic, hierarchic, whatever finite element approximations. Ben 
From: David Knezevic <dknezevic@se...>  20130215 16:54:37

On 02/15/2013 11:48 AM, Kirk, Benjamin (JSCEG311) wrote: > > On Feb 15, 2013, at 10:41 AM, "David Knezevic" <dknezevic@...> wrote: > >> On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: >>> Hi, >>> I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. >> Sure, that should be no problem. >> >> David > The trick is the difference between geometric elements and finite elements. Your mesh will need Quad9 geometric elements, but you can add bilinear, biquadratic, hierarchic, whatever finite element approximations. Yes, sorry my post didn't actually point out the answer! 
From: Subramanya Gautam Sadasiva <ssadasiv@pu...>  20130215 17:01:57

Hi , I knew that that was the way things were implemented. The main problem seems to be that I am using the snesvi solver and to apply the bounds , I do MeshBase::const_node_iterator node_it = _solbub._mesh>local_nodes_begin(); const MeshBase::const_node_iterator node_it_end = _solbub._mesh>local_nodes_end(); for (;node_it!=node_it_end;node_it++) { Node* node = *node_it; const unsigned int dofnum_phi = node>dof_number(sys_num,phi_var,0); const unsigned int dofnum_mu = node>dof_number(sys_num,mu_var,0); } This causes a segmentation fault. I guess the node iterator goes over ALL the nodes that it sees. How do I check if dof_number is returning something meaningful? Sorry my first email was a bit hasty. I sent it before I localized the error. Thanks, Subramanya Sadasiva "But memory eventually fades. Turbulences damp out, internal strains yield to plastic flow, concentration inhomogeneities diffuse to uniformity. Systems tend to subside to very simple states,independent of their specific history." Herbert Callen . Thermodynamics and an Introduction to Thermostatics.  Original Message  From: "David Knezevic" <dknezevic@...> Cc: libmeshusers@... Sent: Friday, February 15, 2013 11:54:28 AM Subject: Re: [Libmeshusers] Multiple equatiosn with different orders of interpolation On 02/15/2013 11:48 AM, Kirk, Benjamin (JSCEG311) wrote: > > On Feb 15, 2013, at 10:41 AM, "David Knezevic" <dknezevic@...> wrote: > >> On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: >>> Hi, >>> I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. >> Sure, that should be no problem. >> >> David > The trick is the difference between geometric elements and finite elements. Your mesh will need Quad9 geometric elements, but you can add bilinear, biquadratic, hierarchic, whatever finite element approximations. Yes, sorry my post didn't actually point out the answer!  Free NextGen Firewall Hardware Offer Buy your Sophos nextgen firewall before the end March 2013 and get the hardware for free! Learn more. http://p.sf.net/sfu/sophosd2dfeb _______________________________________________ Libmeshusers mailing list Libmeshusers@... https://lists.sourceforge.net/lists/listinfo/libmeshusers 
From: Kirk, Benjamin (JSCEG311) <benjamin.kirk1@na...>  20130215 17:42:05

I believe you need to first check DofObject::n_comp() for the node to make sure it's not 0. On Feb 15, 2013, at 11:02 AM, "Subramanya Gautam Sadasiva" <ssadasiv@...> wrote: > Hi , > I knew that that was the way things were implemented. The main problem seems to be that I am using the snesvi solver and to apply the bounds , I do > > > > MeshBase::const_node_iterator node_it = _solbub._mesh>local_nodes_begin(); > const MeshBase::const_node_iterator node_it_end = _solbub._mesh>local_nodes_end(); > for (;node_it!=node_it_end;node_it++) > > { > > Node* node = *node_it; > > const unsigned int dofnum_phi = node>dof_number(sys_num,phi_var,0); > const unsigned int dofnum_mu = node>dof_number(sys_num,mu_var,0); > > > } > > This causes a segmentation fault. I guess the node iterator goes over ALL the nodes that it sees. How do I check if dof_number is returning something meaningful? Sorry my first email was a bit hasty. I sent it before I localized the error. > > Thanks, > > Subramanya Sadasiva > > "But memory eventually fades. Turbulences damp out, internal strains yield to plastic flow, concentration inhomogeneities diffuse to uniformity. Systems tend to subside to very simple states,independent of their specific history." > Herbert Callen . Thermodynamics and an Introduction to Thermostatics. > >  Original Message  > From: "David Knezevic" <dknezevic@...> > Cc: libmeshusers@... > Sent: Friday, February 15, 2013 11:54:28 AM > Subject: Re: [Libmeshusers] Multiple equatiosn with different orders of interpolation > > > On 02/15/2013 11:48 AM, Kirk, Benjamin (JSCEG311) wrote: >> >> On Feb 15, 2013, at 10:41 AM, "David Knezevic" <dknezevic@...> wrote: >> >>> On 02/15/2013 11:39 AM, Subramanya Gautam Sadasiva wrote: >>>> Hi, >>>> I am trying to solve a problem with multiple system objects as part of the same equation syztems object. One of them ( A navier stokes solver ) needs quad9 elements and the other one a cahn hilliard solver, I want to solve using QUAD4 elements. Is this possible on a single mesh?. The code does not work when I try to do this. >>> Sure, that should be no problem. >>> >>> David >> The trick is the difference between geometric elements and finite elements. Your mesh will need Quad9 geometric elements, but you can add bilinear, biquadratic, hierarchic, whatever finite element approximations. > > Yes, sorry my post didn't actually point out the answer! > > > > >  > Free NextGen Firewall Hardware Offer > Buy your Sophos nextgen firewall before the end March 2013 > and get the hardware for free! Learn more. > http://p.sf.net/sfu/sophosd2dfeb > _______________________________________________ > Libmeshusers mailing list > Libmeshusers@... > https://lists.sourceforge.net/lists/listinfo/libmeshusers > >  > Free NextGen Firewall Hardware Offer > Buy your Sophos nextgen firewall before the end March 2013 > and get the hardware for free! Learn more. > http://p.sf.net/sfu/sophosd2dfeb > _______________________________________________ > Libmeshusers mailing list > Libmeshusers@... > https://lists.sourceforge.net/lists/listinfo/libmeshusers 
From: John Peterson <jwpeterson@gm...>  20130215 17:47:49

On Fri, Feb 15, 2013 at 10:41 AM, Kirk, Benjamin (JSCEG311) < benjamin.kirk1@...> wrote: > I believe you need to first check DofObject::n_comp() for the node to make > sure it's not 0. > Running in devel/debug mode will help you catch things like this too. We have libmesh_assert_less (comp, this>n_comp(s,var)); in DofObject::dof_number()  John 
Sign up for the SourceForge newsletter:
No, thanks