There are two possibilities here... You could declare two fe objects for the
velocity components, each using a different quadrature rule. Then for integrals
of the pressure space against the velocity space just choose the right fe.
The other possibility would be to use the optional argument to f
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
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 given to me was to under-integrate the pressure
variable. So, I was thinking to just
If you mean can the library properly under-integrate, say, the mass matrix for
bilinear lagrange using 1pt quadrature no, but it should be an easy add...
That is where the problems come in, right? For example you can solve the
transient heat eqn under-integrated if you do something with the ma
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
> --
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
On Wed, Apr 8, 2009 at 10:55 AM, Stefan wrote:
>
> Hi,
>
> it seems that the triangle module close the whole application when it
> runs out of memory for example.
> Also it seems that Shewchucks code is buggy for quality PSLG meshes with
> a triangle size lower than 0.01
> The problem is in line 7
Hi,
it seems that the triangle module close the whole application when it
runs out of memory for example.
Also it seems that Shewchucks code is buggy for quality PSLG meshes with
a triangle size lower than 0.01
The problem is in line 7369 of the triangle.c
if ((dist < 1.001 * dist2) && (dis
