// Preparations for optional smoothing and layer treatment
    if ((doBoundarySmoothing) or (doLayerTreatment))
    {
        // Resize and initialize storage variables
        prismIslands1.setSize(mesh.nPoints(), UNDEF_LABEL);
        prismIslands2.setSize(mesh.nPoints(), UNDEF_LABEL);
        prismIslands3.setSize(mesh.nPoints(), UNDEF_LABEL);
        pointHops1.setSize(mesh.nPoints(), UNDEF_LABEL);
        pointHops2.setSize(mesh.nPoints(), UNDEF_LABEL);
        pointHops2.setSize(mesh.nPoints(), UNDEF_LABEL);
        pointNormals1.setSize(mesh.nPoints(), Zero);
        pointNormals2.setSize(mesh.nPoints(), Zero);
        pointNormals3.setSize(mesh.nPoints(), Zero);
        innerPrismPoints1.setSize(mesh.nPoints(), UNDEF_VECTOR);
        innerPrismPoints2.setSize(mesh.nPoints(), UNDEF_VECTOR);
        innerPrismPoints3.setSize(mesh.nPoints(), UNDEF_VECTOR);
        outerPrismPoints1.setSize(mesh.nPoints(), UNDEF_VECTOR);
        outerPrismPoints2.setSize(mesh.nPoints(), UNDEF_VECTOR);
        outerPrismPoints3.setSize(mesh.nPoints(), UNDEF_VECTOR);
        innerPrismPointLabels1.setSize(mesh.nPoints(), UNDEF_LABEL);
        innerPrismPointLabels2.setSize(mesh.nPoints(), UNDEF_LABEL);
        innerPrismPointLabels3.setSize(mesh.nPoints(), UNDEF_LABEL);
        outerPrismPointLabels1.setSize(mesh.nPoints(), UNDEF_LABEL);
        outerPrismPointLabels2.setSize(mesh.nPoints(), UNDEF_LABEL);
        outerPrismPointLabels3.setSize(mesh.nPoints(), UNDEF_LABEL);