def solve(self, parameters=default_parameters()):
"Solve and return computed solution (u_F, p_F, U_S, P_S, U_M, P_M)"
# Create submeshes and mappings (only first time)
if self.Omega is None:
# Refine original mesh
mesh = self._original_mesh
for i in range(parameters["num_initial_refinements"]):
mesh = refine(mesh)
# Initialize meshes
self.init_meshes(mesh, parameters)
# Create solver
solver = FSISolver(self)
# Solve
return solver.solve(parameters)
