ffc team mailing list archive

[Vuk Milisic <Vuk.Milisic@xxxxxxx>]

Hello,


If I have a mixed element and I want a scalar weight
for instance x[0], is the following piece of code correct ?

P2 = FiniteElement("Lagrange", "triangle", 2)

P1 = FiniteElement("Lagrange", "triangle", 1)

Vh = (P2 + P2 +P2) 

Th = Vh + P1

(v, q) = TestFunctions(Th)

(u, p) = TrialFunctions(Th)

(phi,pkmu) = Functions (Th)

Py = Function(P2)

a=(phi[i]*u[i]*v[i])*Py+pkmu*p*q



I can't compile that, something's wrong ?


VM



--



-----------------------------------
Vuk Milisic
Chargé de Recherche
LJK-IMAG UMR 5523
51, rue des Mathematiques -  B. P. 53 
38041 Grenoble Cedex 9  France 
Office   n° 66
Tel:  +33 4 76 63 57 38
Fax:  +33 4 76 63 12 63
http://ljk.imag.fr/membres/Vuk.Milisic/
-----------------------------------