mesh.cell_orientations() not editable after intitialization
Issue #1043
new
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 anExpression
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]
Any news on this? We are stuck with exactly the same problem. Interestingly, it used to work in older versions…