fenics team mailing list archive

[Chong Luo <luo.chong@xxxxxxxxx>]

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. 


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







      ____________________________________________________________________________________
Be a better friend, newshound, and 
know-it-all with Yahoo! Mobile.  Try it now.  http://mobile.yahoo.com/;_ylt=Ahu06i62sR8HDtDypao8Wcj9tAcJ