ufl team mailing list archive

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

Anders Logg wrote:
> On Thu, Dec 03, 2009 at 10:37:19AM +0000, Garth N. Wells wrote:
>>
>> 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?
> 
> Both.
> 
> It uses estimate_max to know which degree to choose for unspecified
> elements.
> 

It looks to me that FFC is now using too few quadrature points which is
why the FFC tests are failing.

Garth

> Then estimate_total is used to get the total degree for exact
> integration in the tensor representation. The quadrature scheme will
> use estimate_max here also.
> 
> Also note that estimate_max does not subtract for derivatives.
> 
> --
> Anders