Re: [Libmesh-users] Under Integration

2009-04-08 Thread Kirk, Benjamin (JSC-EG311)
Sent: Wed Apr 08 19:39:39 2009 Subject: Re: [Libmesh-users] Under Integration I guess I should explain what I'm trying to do. =) I have an 8 variable system of which one variable is pressure. When the pressure variable is added, I get a very poorly conditioned matrix. A suggestion that was

Re: [Libmesh-users] Under Integration

2009-04-08 Thread John Peterson
On Wed, Apr 8, 2009 at 6:27 PM, David Knezevic wrote: > You could do: > > AutoPtr qrule = > fe_type.default_quadrature_rule(dim,extra_quad_order); > > where extra_quad_order is negative. If you want to achieve "mass lumping" there are also Trapezoidal (QTrap) rules available. There's also Simpso

Re: [Libmesh-users] Under Integration

2009-04-08 Thread Andrea Hawkins
gt; > - Original Message - > From: John Peterson > To: David Knezevic > Cc: libmesh-users > Sent: Wed Apr 08 18:49:58 2009 > Subject: Re: [Libmesh-users] Under Integration > > On Wed, Apr 8, 2009 at 6:27 PM, David Knezevic wrote: > > You could do: > >

Re: [Libmesh-users] Under Integration

2009-04-08 Thread Kirk, Benjamin (JSC-EG311)
mass matrix? - Original Message - From: John Peterson To: David Knezevic Cc: libmesh-users Sent: Wed Apr 08 18:49:58 2009 Subject: Re: [Libmesh-users] Under Integration On Wed, Apr 8, 2009 at 6:27 PM, David Knezevic wrote: > You could do: > > Auto

Re: [Libmesh-users] Under Integration

2009-04-08 Thread David Knezevic
You could do: AutoPtr qrule = fe_type.default_quadrature_rule(dim,extra_quad_order); where extra_quad_order is negative. - Dave Andrea Hawkins wrote: > Hello- > > I was wondering if there is an easy implementation set up in libMesh for > under-integration. > > Thanks, > Andrea > --

[Libmesh-users] Under Integration

2009-04-08 Thread Andrea Hawkins
Hello- I was wondering if there is an easy implementation set up in libMesh for under-integration. Thanks, Andrea -- This SF.net email is sponsored by: High Quality Requirements in a Collaborative Environment. Download a