Method createWindowSubsurface

src/model/FloorspaceReverseTranslator.cpp:590–706

Creates one or more window subsurfaces along an edge

    // Creates one or more window subsurfaces along an edge
    void FloorspaceReverseTranslator_Impl::createWindowSubsurface(const FSWindow& window, Surface& osSurface, const FSEdge& edge, double minZ,
                                                                  double maxZ) {

      auto windowDefinition = window.windowDefinition();
      if (!windowDefinition.has_value()) {
        return;
      }

      const auto& vertex1 = edge.firstVertex();
      const auto& vertex2 = edge.secondVertex();
      Vector3d edgeVector = edge.edgeVector();

      edgeVector.setLength(1.0);
      Vector3d upVector(0, 0, 1);
      Vector3d crossVector = upVector.cross(edgeVector);

      double sillHeight = windowDefinition->sillHeight();
      double height = 0.0;
      double width = 0.0;
      if ((windowDefinition->windowDefinitionMode() == "Single Window") || (windowDefinition->windowDefinitionMode() == "Repeating Windows")) {
        height = windowDefinition->height();
        width = windowDefinition->width();
      } else if (windowDefinition->windowDefinitionMode() == "Window to Wall Ratio") {
        double length = edge.edgeVector().length();
        width = length - 0.0508;  // Allow for 1" either end of the window
        // Area of the wall * wwr gives area of the window divided by width of the window gives height of the window
        height = (maxZ - minZ) * length * windowDefinition->wwr() / width;
      }

      for (double alpha : window.alphas()) {
        // Calculate the window center from the alpha and sill height
        double x = (1.0 - alpha) * vertex1.x() + alpha * vertex2.x();
        double y = (1.0 - alpha) * vertex1.y() + alpha * vertex2.y();
        Point3d windowCenter(x, y, 0);
        edgeVector.setLength(0.5 * width);
        Point3d window1 = windowCenter + edgeVector;
        edgeVector.setLength(-0.5 * width);
        Point3d window2 = windowCenter + edgeVector;

        Point3dVector windowVertices;
        windowVertices.push_back(Point3d(window1.x(), window1.y(), minZ + sillHeight + height));
        windowVertices.push_back(Point3d(window1.x(), window1.y(), minZ + sillHeight));
        windowVertices.push_back(Point3d(window2.x(), window2.y(), minZ + sillHeight));
        windowVertices.push_back(Point3d(window2.x(), window2.y(), minZ + sillHeight + height));

        // Make sure the sub-surface orientation matches the surface
        Plane plane1(osSurface.vertices());
        Plane plane2(windowVertices);
        if (plane1.reverseEqual(plane2)) {
          std::reverse(windowVertices.begin(), windowVertices.end());
          Point3d temp = window1;
          window1 = window2;
          window2 = temp;
          edgeVector = edgeVector.reverseVector();
          crossVector = crossVector.reverseVector();
        }

        SubSurface subSurface(windowVertices, m_model);