Ejemplo n.º 1
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;
}