Hi ,
I knew that that was the way things were implemented. The main problem seems to
be that I am using the snes-vi 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" <[email protected]>
Cc: [email protected]
Sent: Friday, February 15, 2013 11:54:28 AM
Subject: Re: [Libmesh-users] Multiple equatiosn with different orders of
interpolation
On 02/15/2013 11:48 AM, Kirk, Benjamin (JSC-EG311) wrote:
>
> On Feb 15, 2013, at 10:41 AM, "David Knezevic" <[email protected]>
> 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 Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
------------------------------------------------------------------------------
Free Next-Gen Firewall Hardware Offer
Buy your Sophos next-gen firewall before the end March 2013
and get the hardware for free! Learn more.
http://p.sf.net/sfu/sophos-d2d-feb
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
