mesh.cell_orientations() not editable after intitialization

[Jørgen Dokken]

After the introduction of pybind, it is not possible to change cell orientations after initialization. There are two issues here:

• If this is behavior not wanted, the code below throws an error

• if this is wanted behavior, we restrict the possiblities of cell_orientations to cases that can be described by an Expression

from dolfin import *
mesh = UnitSquareMesh(3,3)
mesh.init_cell_orientations(Expression(("1","1"),degree=1))
print(mesh.cell_orientations())
mesh.cell_orientations()[0] = 1
print(mesh.cell_orientations())

[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

[Stephan Schmidt]

Any news on this? We are stuck with exactly the same problem. Interestingly, it used to work in older versions…