dolfin team mailing list archive

[Johan Hake <johan.hake@xxxxxxxxx>]

On Thursday October 21 2010 01:34:08 Marie Rognes wrote:
> On 21. okt. 2010 00:34, Johan Hake wrote:
> > Btw:
> > 
> > Is the rank function in ufl borked:
> >   mesh = UnitSquare(1, 1)
> >   V = FunctionSpace(mesh, "CG", 1)
> >   v = TestFunction(V)
> >   u = Function(V)
> >   L = u*v*dx
> >   rank(L)
> >   => AttributeError: 'Form' object has no attribute 'shape'
> > 
> > or do I need to preprocess the form first?
> 
> rank operates on expr's (not on forms). 

Ahh...

> An input check might be
> appropriate.

Yes!

> Try this instead
> 
>     from ufl.algorithms import compute_form_arities
>     ranks = compute_form_arities(L)

It works on forms, but not on a linearized form

  du = TrialFunction(V)
  compute_form_arities(derivative(L, u, du))
  => UFLException: Found basis function in ... this is an invalid expression.

Johan

> --
> Marie
> 
> _______________________________________________
> Mailing list: https://launchpad.net/~dolfin
> Post to     : dolfin@xxxxxxxxxxxxxxxxxxx
> Unsubscribe : https://launchpad.net/~dolfin
> More help   : https://help.launchpad.net/ListHelp