added functions to get outer edges

This commit is contained in:
jakob.schratter 2026-01-24 13:44:41 +01:00
commit db7bd92b2d

View file

@ -774,9 +774,9 @@ const std::vector<int> Mesh::OuterEdges() const
cout << "All boundary edges: " << _ebedges.size() << endl; cout << "All boundary edges: " << _ebedges.size() << endl;
cout << "Outer boundary edges: " << outerEdges.size() << endl; cout << "Outer boundary edges: " << outerEdges.size() << endl;
cout << _ebedges; // cout << _ebedges;
cout << endl; // cout << endl;
cout << outerEdges << endl; // cout << outerEdges << endl;
return outerEdges; return outerEdges;
} }
@ -1017,22 +1017,8 @@ Mesh::Mesh(std::string const &fname)
Mesh::Mesh(std::string const &filename, std::string const &subdomain_filename) : Mesh(filename) Mesh::Mesh(std::string const &filename, std::string const &subdomain_filename) : Mesh(filename)
{ {
ElementSubdomains = ReadElementSubdomains(subdomain_filename); ElementSubdomains = ReadElementSubdomains(subdomain_filename);
//EdgeSubdomains = ReadEdgeSubdomains(filename);
} }
// const vector<int> Mesh::ReadEdgeSubdomains(std::string const &filename) const
// {
// vector<int> edgeSubdomains(_nedge);
// ifstream ifs(filename);
// if (!(ifs.is_open() && ifs.good())) {
// cerr << "Mesh::ReadEdgeSubdomains: Error cannot open file " << filename << endl;
// assert(ifs.is_open());
// }
// return edgeSubdomains;
// }
const vector<int> Mesh::ReadElementSubdomains(string const &dname) const const vector<int> Mesh::ReadElementSubdomains(string const &dname) const
{ {