fenics team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

On Sat, Mar 22, 2008 at 04:50:13PM -0700, Chong Luo wrote:
> Ah, it worked after defining (u,p) as a system, like this:
> 
> deg=1
> vector = VectorElement("Lagrange", "tetrahedron", deg)
> scalar = FiniteElement("Lagrange", "tetrahedron", deg)
> system = vector + scalar
> 
> (v, q) = TestFunctions(system)
> (u, p) = TrialFunctions(system)
> 
> Sorry for creating so much garbage in the mailing list.
> Maybe I should think more before I post to the mailing list.

No problem. Good to see that it works. For further questions about
implementing forms, try posting to ffc-dev@xxxxxxxxxx instead.

Thanks
-- 
Anders


> 
> Best,
> Chong Luo
> ----- Original Message ----
> From: Chong Luo <luo.chong@xxxxxxxxx>
> To: FEniCS-dev@xxxxxxxxxx
> Sent: Saturday, March 22, 2008 5:54:07 PM
> Subject: [FEniCS-dev] IndexError: index (4) out of range (0<=index<=4) in
> dimension 1
> 
> Hi,
> 
> After careful thinking, I almost got the form compiled, except one error which
> I don't understand:
>   File "/home/aharon/local/lib/python2.4/site-packages/ffc/compiler/
> codegeneration/tensor/tensorgenerator.py", line 296, in
> __generate_element_tensor
>     a0 = A0.A0[tuple(i + a)]
> IndexError: index (4) out of range (0<=index<=4) in dimension 1
> 
> 
> And the form file is in the attachment, which I wrote to solve incompressible
> nonlinear elasticity problem.
> 
> Thank you!
> 
> Best,
> Chong Luo
> 
> 
> ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
> Be a better friend, newshound, and know-it-all with Yahoo! Mobile. Try it now.
> 
> 
> -----Inline Attachment Follows-----
> 
> _______________________________________________
> FEniCS-dev mailing list
> FEniCS-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/fenics-dev
> 
> 
> ━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━
> Looking for last minute shopping deals? Find them fast with Yahoo! Search.

> _______________________________________________
> FEniCS-dev mailing list
> FEniCS-dev@xxxxxxxxxx
> http://www.fenics.org/mailman/listinfo/fenics-dev