ufl team mailing list archive

[Martin Sandve Alnæs <martinal@xxxxxxxxx>]

Further simplification that still reproduces the problem:


cell = triangle
d = cell.d
I = Identity(d)

vector = VectorElement("Lagrange", cell, 2)
scalar = FiniteElement("Lagrange", cell, 1)
element = vector + scalar

v, q = TestFunctions(element)
u, p = TrialFunctions(element)

F = inner(grad(v), grad(u) - p*I)*dx

#a = lhs(F)
L = rhs(F)


Note that you didn't use a mixed element.

Martin



On Tue, Sep 22, 2009 at 9:16 PM, Harish Narayanan
<harish.mlists@xxxxxxxxx> wrote:
> Martin Sandve Alnæs wrote:
>> I agree that the error is completely incomprehensible,
>> but in sigma(u, p0) you're adding terms with a trialfunction
>> only in one of the terms. Can you try with a valid form?
>
> I was overzealous in reducing the size of the form. As you noted in your
> next e-mail, changing p0 to p results in the same error.
>
> Harish
>