diff --git a/src/mesh/mesh_base.C b/src/mesh/mesh_base.C index 9c4e8b88e36..8d5bbae2e10 100644 --- a/src/mesh/mesh_base.C +++ b/src/mesh/mesh_base.C @@ -1159,7 +1159,7 @@ std::string MeshBase::get_info(const unsigned int verbosity /* = 0 */, const boo if (verbosity > 1) { this->comm().min(info.bbox.min()); - this->comm().min(info.bbox.max()); + this->comm().max(info.bbox.max()); } } diff --git a/src/mesh/poly2tri_triangulator.C b/src/mesh/poly2tri_triangulator.C index 0120668815d..8c144392103 100644 --- a/src/mesh/poly2tri_triangulator.C +++ b/src/mesh/poly2tri_triangulator.C @@ -292,9 +292,12 @@ void Poly2TriTriangulator::triangulate_current_points() std::vector outer_boundary_points; std::vector> inner_hole_points(n_holes); + dof_id_type nn = _mesh.max_node_id(); + libmesh_error_msg_if + (!nn, "Poly2TriTriangulator cannot triangulate an empty mesh!"); + // Unless we're using an explicit segments list, we assume node ids // are contiguous here. - dof_id_type nn = _mesh.max_node_id(); if (this->segments.empty()) libmesh_error_msg_if (_mesh.n_nodes() != nn, diff --git a/src/mesh/unstructured_mesh.C b/src/mesh/unstructured_mesh.C index d7b2baddbd1..24de140b0f3 100644 --- a/src/mesh/unstructured_mesh.C +++ b/src/mesh/unstructured_mesh.C @@ -1773,6 +1773,11 @@ void UnstructuredMesh::stitching_helper (const MeshBase * other_mesh, bool enforce_all_nodes_match_on_boundaries, bool skip_find_neighbors) { +#ifdef DEBUG + // We rely on neighbor links here + MeshTools::libmesh_assert_valid_neighbors(*this); +#endif + // FIXME: make distributed mesh support efficient. // Yes, we currently suck. MeshSerializer serialize(*this); @@ -1935,9 +1940,13 @@ void UnstructuredMesh::stitching_helper (const MeshBase * other_mesh, { libMesh::out << "In UnstructuredMesh::stitch_meshes:\n" << "This mesh has " << this_boundary_node_ids.size() - << " nodes on boundary " << this_mesh_boundary_id << ".\n" + << " nodes on boundary `" + << this->get_boundary_info().get_sideset_name(this_mesh_boundary_id) + << "' (" << this_mesh_boundary_id << ").\n" << "Other mesh has " << other_boundary_node_ids.size() - << " nodes on boundary " << other_mesh_boundary_id << ".\n"; + << " nodes on boundary `" + << this->get_boundary_info().get_sideset_name(other_mesh_boundary_id) + << "' (" << other_mesh_boundary_id << ").\n"; if (h_min_updated) {