Cannot assemble simple Expression (norm of Expression not working)

[Johan Hake]

from dolfin import *
f = Expression("sin(3.0*x[0])*sin(3.0*x[1])*sin(3.0*x[2])")
mesh = UnitCube(2,2,2)
n = norm(f, mesh=mesh)

Raise:
UFLException: Unable to extract form data. Reason: Missing cell definition in form.

[Anders Logg (logg)]

On Wed, Dec 16, 2009 at 06:59:33PM -0000, Johan Hake wrote:
> Public bug reported:
>
> from dolfin import *
> f = Expression("sin(3.0*x[0])*sin(3.0*x[1])*sin(3.0*x[2])")
> mesh = UnitCube(2,2,2)
> n = norm(f, mesh=mesh)
>
> Raise:
> UFLException: Unable to extract form data. Reason: Missing cell definition in form.
>
> ** Affects: dolfin
> Importance: Undecided
> Status: New

I'm on it.

--
Anders

[Anders Logg (logg)]

I have added an optional argument cell to Expression which means you
can do this:

  from dolfin import *
  f = Expression("sin(3.0*x[0])*sin(3.0*x[1])*sin(3.0*x[2])", cell="tetrahedron")
  mesh = UnitCube(2,2,2)
  n = norm(f, mesh=mesh)

But we should ultimately figure out the cell in assemble. The mesh
argument is passed there and it could figure out the cell type and
send that along to the JIT compiler. But that will have to wait until
the new reorganized version of FFC is ready.

--
Anders

On Wed, Dec 16, 2009 at 06:59:33PM -0000, Johan Hake wrote:
> Public bug reported:
>
> from dolfin import *
> f = Expression("sin(3.0*x[0])*sin(3.0*x[1])*sin(3.0*x[2])")
> mesh = UnitCube(2,2,2)
> n = norm(f, mesh=mesh)
>
> Raise:
> UFLException: Unable to extract form data. Reason: Missing cell definition in form.
>
> ** Affects: dolfin
> Importance: Undecided
> Status: New
>

[Garth Wells (garth-wells)]

Is this bug resolved now?

[Anders Logg (logg)]

Should be fixed now.