That works fine when you have one Testfunction like "v", How about when
you have the following set of Testfunction:
-->
P2 = VectorElement("Lagrange", "triangle", 1)
P1 = FiniteElement("Lagrange", "triangle", 1)
TH = P2 + P1(v,w) = TestFunctions(TH) --> analogous with what anders said, one should be able to do something like : (v,w)=( v + delta*A(v), w) However, this creates problem in ffc. An alternative is to multiply the whole equation sets by the added stabilization term "delta*A(v)" and add the new terms to the unstabilized variational form. But it is much easier if one can just change the weighing functions v to v+delta*A(v). is there anyway to do that?
--------------------------------------
Hatef Monajemi
Graduate Student
Dept of Civil and Environmental Eng
2032 MC Building, Carleton University
1125 Colonel By Drive, Ottawa
Ontario, Canada, K1S 5B6
_______________________________________________ DOLFIN-dev mailing list [email protected] http://www.fenics.org/mailman/listinfo/dolfin-dev
