Re: [Libmesh-users] XYZ spatial coordinates

2008-03-04 Thread Roy Stogner
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

Re: [Libmesh-users] XYZ spatial coordinates

2008-03-04 Thread Maxime Debon
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

Re: [Libmesh-users] XYZ spatial coordinates

2008-03-03 Thread Roy Stogner
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

Re: [Libmesh-users] XYZ spatial coordinates

2008-03-03 Thread Maxime Debon
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

[Libmesh-users] XYZ spatial coordinates

2008-02-27 Thread Maxime Debon
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();