Method translateSurface

src/sdd/MapGeometry.cpp:1242–1472

  }

  boost::optional<model::ModelObject> ReverseTranslator::translateSurface(const pugi::xml_node& element, openstudio::model::Space& space) {
    boost::optional<model::ModelObject> result;

    UnitSystem siSys(UnitSystem::SI);

    std::vector<openstudio::Point3d> vertices;

    pugi::xml_node polyLoopElement = element.child("PolyLp");
    if (!polyLoopElement) {
      LOG(Error, "Surface element 'PolyLp' is empty, cannot create Surface.");
      return boost::none;
    }

    std::vector<pugi::xml_node> cartesianPointElements = makeVectorOfChildren(polyLoopElement, "CartesianPt");
    for (std::vector<pugi::xml_node>::size_type i = 0; i < cartesianPointElements.size(); i++) {
      std::vector<pugi::xml_node> coordinateElements = makeVectorOfChildren(cartesianPointElements[i], "Coord");
      if (coordinateElements.size() != 3) {
        LOG(Error, "PolyLp element 'CartesianPt' does not have exactly 3 'Coord' elements, cannot create Surface.");
        return boost::none;
      }

      /* DLM: these unit conversions are taking about 75% of the total time to translate a large model

      // sdd units = ft, os units = m
      Quantity xIP(coordinateElements.at(0).toElement().text().as_double(), IPUnit(IPExpnt(0,1,0)));
      Quantity yIP(coordinateElements.at(1).toElement().text().as_double(), IPUnit(IPExpnt(0,1,0)));
      Quantity zIP(coordinateElements.at(2).toElement().text().as_double(), IPUnit(IPExpnt(0,1,0)));

      OptionalQuantity xSI = QuantityConverter::instance().convert(xIP, siSys);
      OS_ASSERT(xSI);
      OS_ASSERT(xSI->units() == SIUnit(SIExpnt(0,1,0)));

      OptionalQuantity ySI = QuantityConverter::instance().convert(yIP, siSys);
      OS_ASSERT(ySI);
      OS_ASSERT(ySI->units() == SIUnit(SIExpnt(0,1,0)));

      OptionalQuantity zSI = QuantityConverter::instance().convert(zIP, siSys);
      OS_ASSERT(zSI);
      OS_ASSERT(zSI->units() == SIUnit(SIExpnt(0,1,0)));

      vertices.push_back(openstudio::Point3d(xSI->value(), ySI->value(), zSI->value()));
      */

      double x = footToMeter * coordinateElements[0].text().as_double();
      double y = footToMeter * coordinateElements[1].text().as_double();
      double z = footToMeter * coordinateElements[2].text().as_double();
      vertices.push_back(openstudio::Point3d(x, y, z));
    }

    openstudio::model::Surface surface(vertices, space.model());
    surface.setSpace(space);

    pugi::xml_node nameElement = element.child("Name");
    std::string name;
    if (!nameElement) {
      LOG(Error, "Surface element 'Name' is empty.")
    } else {
      name = escapeName(nameElement.text().as_string());