dolfin team mailing list archive

[Anders Logg <logg@xxxxxxxxx>]

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
>

-- 
Cannot assemble simple Expression (norm of Expression not working)
https://bugs.launchpad.net/bugs/497509
You received this bug notification because you are a member of DOLFIN
Team, which is subscribed to DOLFIN.

Status in DOLFIN: New

Bug description:
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.