On Tue, 4 Mar 2008, Maxime Debon wrote:
> Looking for tips to check the correctness of Ke matrix and Fe vector
> in an handbook, I found this rule : A0 * u'' = (A0 * u')' - A0' * u'
>
> This latest transformation must be taken into account to adjust a
> non-constant coefficient to the diffusion t
Hi,
Looking for tips to check the correctness of Ke matrix and Fe vector
in an handbook, I found this rule : A0 * u'' = (A0 * u')' - A0' * u'
This latest transformation must be taken into account to adjust a
non-constant coefficient to the diffusion term. It also explains why
the result was
On Mon, 3 Mar 2008, Maxime Debon wrote:
> Taking the equation : x * u''(x) = 1
> between [0.5 ; 1.0]
> with u[0.5] = 0.5 and u[1.0] = 1.0,
>
> we would expect a value of 0.88045 at the point 0.900391 but we obtain
> 0.900391.
This problem is still giving you trouble? I recall seeing your
earlie
Hi,
To give a precision, the variable coefficient only generate an error
computing the diffusion term (u''(x)). In contrast, tests made with an
X dependent coordinate coefficient on convection term and u(x) term
gave a very powerful precision. Thus, the error can't be attributed to
previou
Hi,
Using spatial coordinates (q_point) in matrix Ke almost as it's
presented in the example 14 for matrix Fe, I may misunderstand
something important. Here is the code.
// -- //
const std::vector& qrule_points = fe->get_xyz();