SyFi

  1. Bug #997739
  2. Comment #0

Comment 0 for bug 997739

Revision history for this message
Johan Hake (johan-hake) wrote :

from dolfin import *
mesh = UnitCube(5,5,5)
VV = VectorFunctionSpace(mesh,"Lagrange",1)
u = TrialFunction(VV)
v = TestFunction(VV)
Aij = diag(Constant((1,1,1)))
Delta = Identity( mesh.ufl_cell().geometric_dimension())
# Does not work (but works with sfc)
assemble(inner(Aij*(grad(u)+Delta), grad(v))*dx)
# Works
assemble(inner(Aij*grad(u), grad(v))*dx)