int main()
{
// Structure sub domain
class Structure : public SubDomain
{
bool inside(const Array<double>& x, bool on_boundary) const
{
return x[0] > 1.4 - DOLFIN_EPS and x[0] < 1.6 + DOLFIN_EPS
and x[1] < 0.6 + DOLFIN_EPS;
}
};
// Create mesh
auto mesh = std::make_shared<RectangleMesh>(Point(0.0, 0.0),
Point(3.0, 1.0), 60, 20);
// Create sub domain markers and mark everything as 0
MeshFunction<std::size_t> sub_domains(mesh, mesh->topology().dim());
sub_domains = 0;
// Mark structure domain as 1
Structure structure;
structure.mark(sub_domains, 1);
// Extract sub meshes
auto fluid_mesh = std::make_shared<SubMesh>(*mesh, sub_domains, 0);
auto structure_mesh = std::make_shared<SubMesh>(*mesh, sub_domains, 1);
// Move structure mesh
MeshGeometry& geometry = structure_mesh->geometry();
for (VertexIterator v(*structure_mesh); !v.end(); ++v)
{
std::array<double, 2> x = {{v->x()[0], v->x()[1]}};
x[0] += 0.1*x[0]*x[1];
geometry.set(v->index(), x.data());
}
// Move fluid mesh according to structure mesh
ALE::move(fluid_mesh, *structure_mesh);
fluid_mesh->smooth();
// Plot meshes
plot(fluid_mesh);
plot(structure_mesh);
interactive();
return 0;
}
int main()
{
// Structure sub domain
class Structure : public SubDomain
{
bool inside(const Array<double>& x, bool on_boundary) const
{
return x[0] > 1.4 - DOLFIN_EPS and x[0] < 1.6 + DOLFIN_EPS and x[1] < 0.6 + DOLFIN_EPS;
}
};
// Create mesh
RectangleMesh mesh(0.0, 0.0, 3.0, 1.0, 60, 20);
// Create sub domain markers and mark everything as 0
MeshFunction<std::size_t> sub_domains(mesh, mesh.topology().dim());
sub_domains = 0;
// Mark structure domain as 1
Structure structure;
structure.mark(sub_domains, 1);
// Extract sub meshes
SubMesh fluid_mesh(mesh, sub_domains, 0);
SubMesh structure_mesh(mesh, sub_domains, 1);
// Move structure mesh
MeshGeometry& geometry = structure_mesh.geometry();
for (VertexIterator v(structure_mesh); !v.end(); ++v)
{
const double* x = v->x();
geometry.x(v->index())[0] += 0.1*x[0]*x[1];
}
// Move fluid mesh according to structure mesh
fluid_mesh.move(structure_mesh);
fluid_mesh.smooth();
// Plot meshes
plot(fluid_mesh);
plot(structure_mesh);
interactive();
return 0;
}