int main(int argc, char* argv[])
{
plbInit(&argc, &argv);
global::directories().setOutputDir("./");
global::IOpolicy().activateParallelIO(false);
if (argc != 9) {
pcout << "Wrong parameters; the syntax is:" << std::endl
<< (std::string)global::argv(0)
<< "[FLT | DBL | LDBL | INF] inputSTL.stl cut-plane-point-x cut-plane-point-y cut-plane-point-z" << std::endl
<< "cut-plane-normal-x cut-plane-normal-y cut-plane-normal-z" << std::endl;
exit(-1);
}
std::string precisionStr;
std::string inStlFileName;
Plane<T> cutPlane;
try {
global::argv(1).read(precisionStr);
global::argv(2).read(inStlFileName);
global::argv(3).read(cutPlane.point[0]);
global::argv(4).read(cutPlane.point[1]);
global::argv(5).read(cutPlane.point[2]);
global::argv(6).read(cutPlane.normal[0]);
global::argv(7).read(cutPlane.normal[1]);
global::argv(8).read(cutPlane.normal[2]);
}
catch (PlbIOException& exception) {
pcout << "Wrong parameters." << std::endl;
exit(-1);
}
Precision precision;
if (precisionStr == "FLT") {
precision = FLT;
} else if (precisionStr == "DBL") {
precision = DBL;
} else if (precisionStr == "LDBL") {
precision = LDBL;
} else if (precisionStr == "INF") {
precision = INF;
} else {
pcout << "Wrong precision command-line argument." << std::endl;
exit(-1);
}
TriangleSet<T>* set = 0;
try {
set = new TriangleSet<T>(inStlFileName, precision);
}
catch (PlbIOException& exception) {
pcout << "ERROR, could not read STL file " << inStlFileName
<< ": " << exception.what() << std::endl;
exit(-1);
}
// Cut the geometry in two parts.
TriangleSet<T> normalNegativePart;
int rv = 0;
rv = set->cutWithPlane(cutPlane, normalNegativePart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
normalNegativePart.writeBinarySTL("normalNegativePart.stl");
TriangleSet<T> normalPositivePart;
cutPlane.normal = -cutPlane.normal;
rv = 0;
rv = set->cutWithPlane(cutPlane, normalPositivePart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
normalPositivePart.writeBinarySTL("normalPositivePart.stl");
delete set;
return 0;
}
int main(int argc, char* argv[])
{
plbInit(&argc, &argv);
global::directories().setOutputDir("./");
plint n = 80; // Resolution (cylinder diameter).
T radius = (T) n / 2.; // Cylinder radius.
Array<T,3> center(0.0, 0.0, 0.0); // Cylinder center.
T length = (T) (4 * n); // Cylinder length.
plint nAxial = n / 2; // Two parameters needed for the creation of the triangularized cylinder surface.
plint nCirc = 3 * n / 2;
// Create a cylinder surface as a set of triangles.
TriangleSet<T> cylinder;
cylinder = constructCylinder<T>(center, radius, radius, length, nAxial, nCirc);
cylinder.writeAsciiSTL("cylinder.stl");
// Define the cut plane by a point and a normal vector.
Array<T,3> cutPlanePoint(length / 3.3, 0.0, 0.0);
Array<T,3> cutPlaneNormal(-1.0, -0.1, -0.5);
Plane<T> cutPlane;
cutPlane.point = cutPlanePoint;
cutPlane.normal = cutPlaneNormal;
// Cut the cylinder in two parts.
TriangleSet<T> rightCylinderPart;
int rv = 0;
rv = cylinder.cutWithPlane(cutPlane, rightCylinderPart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
rightCylinderPart.writeAsciiSTL("rightCylinderPart.stl");
TriangleSet<T> leftCylinderPart;
cutPlane.normal = -cutPlane.normal;
rv = 0;
rv = cylinder.cutWithPlane(cutPlane, leftCylinderPart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
leftCylinderPart.writeAsciiSTL("leftCylinderPart.stl");
// Use a cuboid to cut the cylinder (correct usage).
Array<T,3> lowerLeftCorner = cutPlanePoint - Array<T,3>(length, 5.0 * radius / 4.0, 5.0 * radius / 4.0);
Array<T,3> upperRightCorner = cutPlanePoint + Array<T,3>(length / 4.0, 5.0 * radius / 4.0, 5.0 * radius / 4.0);
Cuboid<T> cutCuboid;
cutCuboid.lowerLeftCorner = lowerLeftCorner;
cutCuboid.upperRightCorner = upperRightCorner;
TriangleSet<T> newRightCylinderPart;
cutPlane.normal = -cutPlane.normal;
rv = 0;
rv = cylinder.cutWithPlane(cutPlane, cutCuboid, newRightCylinderPart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
newRightCylinderPart.writeAsciiSTL("newRightCylinderPart.stl");
// Use a cuboid to cut the cylinder (wrong usage).
upperRightCorner = cutPlanePoint + Array<T,3>(0.0, radius, radius);
cutCuboid.upperRightCorner = upperRightCorner;
TriangleSet<T> newBrokenRightCylinderPart;
rv = 0;
rv = cylinder.cutWithPlane(cutPlane, cutCuboid, newBrokenRightCylinderPart);
if (rv != 1) {
pcout << "Problem with the surface cutting." << std::endl;
exit(1);
}
newBrokenRightCylinderPart.writeAsciiSTL("newBrokenRightCylinderPart.stl");
return 0;
}
