ufl team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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.

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

--
Anders

Attachment: signature.asc
Description: Digital signature