ufl team mailing list archive

["Garth N. Wells" <gnw20@xxxxxxxxx>]

Anders Logg wrote:
> On Thu, Dec 03, 2009 at 10:20:06AM +0000, Garth N. Wells wrote:
>> The test for
>>
>>   estimate_max_polynomial_degree
>>
>> is failing. For
>>
>>         V2 = FiniteElement("CG", triangle, 2)
>>         VV = VectorElement("CG", triangle, 3)
>>         v2 = BasisFunction(V2)
>>         vv = BasisFunction(VV)
>>         print "****Degree ", estimate_max_polynomial_degree(v2*vv[0])
>>
>> UFL computes a max degree of 3, when it should obviously be 5.
>>
>> Anders: is this related the changes that you recently made?
> 
> Yes. The unit test should be changed.
> 
>   estimate_max_polynominal_degree
> 
> now actually returns the maximum polynomial degree which is 3.
> 
> I have added a new function
> 
>   estimate_total_polynomial_degree
> 
> which does what estimate_max_polynomial_degree used to do.
> That function will return 5.
>

OK. Which one do we want FFC to use?

Garth

> I'll update it later if you haven't fixed it.
> 
> --
> Anders