- changed status to resolved
MeshCoordinates broken
Issue #760
resolved
This testcase segfaults for me:
mesh = UnitSquareMesh(10,10)
POS = MeshCoordinates(mesh)
V = VectorFunctionSpace(mesh, "CG", 1)
v = Function(V)
v.interpolate(POS)
This is because the constructor of MeshCoordinates uses _mesh
initilizing the constructor of the Expression
base class constructor before initializing it :
MeshCoordinates::MeshCoordinates(std::shared_ptr<const Mesh> mesh)
: Expression(_mesh->geometry().dim()), _mesh(mesh)
{
// Do nothing
}
Comments (1)
-
- Log in to comment
Fixed in ed87e96.