CBoundingBox CMeshBox::GenerateBoundingBox(void)
{
CBoundingBox testbox;
nova::nReal xmax = mMeshDef.nVertexList[0].x,
xmin = mMeshDef.nVertexList[0].x,
ymax = mMeshDef.nVertexList[0].y,
ymin = mMeshDef.nVertexList[0].y,
zmax = mMeshDef.nVertexList[0].z,
zmin = mMeshDef.nVertexList[0].z;
TVertexes::iterator it = mMeshDef.nVertexList.begin()+1;
for(; it != mMeshDef.nVertexList.end(); ++it)
{
xmax = std::max((*it).x, xmax);
xmin = std::min((*it).x, xmin);
ymax = std::max((*it).y, ymax);
ymin = std::min((*it).y, ymin);
zmax = std::max((*it).z, zmax);
zmin = std::min((*it).z, zmin);
}
nova::Vector3f maxpoint(xmax, ymax, zmax);
nova::Vector3f minpoint(xmin, ymin, zmin);
testbox.CalcSides(minpoint, maxpoint);
for(nova::nUInt32 i = 0; i < GetListenersCount(); i++)
{
CMeshBoxListener * lis =
dynamic_cast<CMeshBoxListener *>(GetListener(i));
lis->GenerateBoundingBoxListener(this);
}
return testbox;
}
int main(int argc, char** argv)
{
LibMeshInit init(argc, argv);
GetPot cl(argc, argv);
int dim = -1;
if (!cl.search("--dim"))
{
std::cerr << "No --dim argument found!" << std::endl;
usage_error(argv[0]);
}
dim = cl.next(dim);
Mesh mesh(dim);
if(!cl.search("--input"))
{
std::cerr << "No --input argument found!" << std::endl;
usage_error(argv[0]);
}
const char* meshname = cl.next("mesh.xda");
mesh.read(meshname);
std::cout << "Loaded mesh " << meshname << std::endl;
if(!cl.search("--newbcid"))
{
std::cerr << "No --bcid argument found!" << std::endl;
usage_error(argv[0]);
}
boundary_id_type bcid = 0;
bcid = cl.next(bcid);
Point minnormal(-std::numeric_limits<Real>::max(),
-std::numeric_limits<Real>::max(),
-std::numeric_limits<Real>::max());
Point maxnormal(std::numeric_limits<Real>::max(),
std::numeric_limits<Real>::max(),
std::numeric_limits<Real>::max());
Point minpoint(-std::numeric_limits<Real>::max(),
-std::numeric_limits<Real>::max(),
-std::numeric_limits<Real>::max());
Point maxpoint(std::numeric_limits<Real>::max(),
std::numeric_limits<Real>::max(),
std::numeric_limits<Real>::max());
if (cl.search("--minnormalx"))
minnormal(0) = cl.next(minnormal(0));
if (cl.search("--minnormalx"))
minnormal(0) = cl.next(minnormal(0));
if (cl.search("--maxnormalx"))
maxnormal(0) = cl.next(maxnormal(0));
if (cl.search("--minnormaly"))
minnormal(1) = cl.next(minnormal(1));
if (cl.search("--maxnormaly"))
maxnormal(1) = cl.next(maxnormal(1));
if (cl.search("--minnormalz"))
minnormal(2) = cl.next(minnormal(2));
if (cl.search("--maxnormalz"))
maxnormal(2) = cl.next(maxnormal(2));
if (cl.search("--minpointx"))
minpoint(0) = cl.next(minpoint(0));
if (cl.search("--maxpointx"))
maxpoint(0) = cl.next(maxpoint(0));
if (cl.search("--minpointy"))
minpoint(1) = cl.next(minpoint(1));
if (cl.search("--maxpointy"))
maxpoint(1) = cl.next(maxpoint(1));
if (cl.search("--minpointz"))
minpoint(2) = cl.next(minpoint(2));
if (cl.search("--maxpointz"))
maxpoint(2) = cl.next(maxpoint(2));
std::cout << "min point = " << minpoint << std::endl;
std::cout << "max point = " << maxpoint << std::endl;
std::cout << "min normal = " << minnormal << std::endl;
std::cout << "max normal = " << maxnormal << std::endl;
bool matcholdbcid = false;
boundary_id_type oldbcid = 0;
if (cl.search("--oldbcid"))
{
matcholdbcid = true;
oldbcid = cl.next(oldbcid);
if (oldbcid < 0)
oldbcid = BoundaryInfo::invalid_id;
}
AutoPtr<FEBase> fe = FEBase::build(dim, FEType(FIRST,LAGRANGE));
QGauss qface(dim-1, CONSTANT);
fe->attach_quadrature_rule(&qface);
const std::vector<Point> &face_points = fe->get_xyz();
const std::vector<Point> &face_normals = fe->get_normals();
MeshBase::element_iterator el = mesh.elements_begin();
const MeshBase::element_iterator end_el = mesh.elements_end();
for (; el != end_el; ++el)
{
Elem *elem = *el;
unsigned int n_sides = elem->n_sides();
for (unsigned int s=0; s != n_sides; ++s)
{
if (elem->neighbor(s))
continue;
fe->reinit(elem,s);
const Point &p = face_points[0];
const Point &n = face_normals[0];
//std::cout << "elem = " << elem->id() << std::endl;
//std::cout << "centroid = " << elem->centroid() << std::endl;
//std::cout << "p = " << p << std::endl;
//std::cout << "n = " << n << std::endl;
if (p(0) > minpoint(0) && p(0) < maxpoint(0) &&
p(1) > minpoint(1) && p(1) < maxpoint(1) &&
p(2) > minpoint(2) && p(2) < maxpoint(2) &&
n(0) > minnormal(0) && n(0) < maxnormal(0) &&
n(1) > minnormal(1) && n(1) < maxnormal(1) &&
n(2) > minnormal(2) && n(2) < maxnormal(2))
{
if (matcholdbcid &&
mesh.boundary_info->boundary_id(elem, s) != oldbcid)
continue;
mesh.boundary_info->remove_side(elem, s);
mesh.boundary_info->add_side(elem, s, bcid);
//std::cout << "Set element " << elem->id() << " side " << s <<
// " to boundary " << bcid << std::endl;
}
}
}
std::string outputname;
if(cl.search("--output"))
{
outputname = cl.next("mesh.xda");
}
else
{
outputname = "new.";
outputname += meshname;
}
mesh.write(outputname.c_str());
std::cout << "Wrote mesh " << outputname << std::endl;
return 0;
}
